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The recent interest in aspects common to quantum information and condensed matter has 
prompted a flory of activity at the border of these disciplines that were far distant untill few 
years ago. Numerous interesting questions have been addressed so far. Here we review an im- 
portant part of this field, the properties of the entanglement in many-body systems. We discuss 
the zero and finite temperature properties of entanglement in interacting spin, fermion and boson 
model systems. Both bipartite and multipartite entanglement will be considered. In equilibrium 
we show how entanglement is tightly connected to the characteristics of the phase diagram. The 
behavior of entanglement can be related, via certain witnesses, to thermodynamic quantities thus 
offering interesting possibilities for an experimental test. Out of equilibrium we discuss how to 
generate and manipulate entangled states by means of many-body Hamiltonians. 
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I. INTRODUCTION 

Entanglement expresses the "spooky" no n-locality in- 
herent to quantum mechanics ( Belli . Il987l ). Because of 
that, it gave rise to severe skepticisms since the early 
days of quantum mechanics. It was only after the semi- 
nal contribution of John Bell that the fundamental ques- 
tions related to the existence of entangled states could 
be tested experimentally. In fact, under fairly general 
assumptions, Bell derived a set of inequalities for corre- 
lated measurements of two physical observables that any 
local theory should obey. The overwhelming majority 
of experiments done so far are in favor of quantum me- 
chanics thus demo nstrating tha t quantum entanglement 
is physical reality ( Peresl Il9930 . 1 

Entanglement has gained renewed interest with 
the development of quantu m information sci- 
ence ([Nielsen and Chuanel . l2000f ). In its framework, 
quantum entanglement is viewed at as a precious 
resource in quantum information processing. It is e.g. 
believed to be the main ingredient of the quantum 
speed-up in quantum computation and communication. 
More over several quant um protocols, as teleporta- 
tion (Be nnett et al. I. I1993D just to mention an important 
example, can be realized exclusively with the help of 
entangled states. 

The role of entanglement as a resource in quan- 
tum information has stimulated intensive research try- 
ing to unveil both its qualitativ e and quant i tative 
aspect s dBengtsson and Zvczkowskl 20061: iBrufj, 12002 ; 
EisertJ. 120061: iHorodecki et all 120071: (pienio and Vedral . 
1998 [ iPlenio and Virmanil l2007t IVedrall 120021: IWootters! . 
20011 ). To this end, necessary criteria for any entangle- 



ment measure to be fulfilled have been elaborat ed and 
lead to the notion of an entanglement monotone ( Vidall . 



1 There are states that do not violate Bell inequal ities and never- 
theless arc entangled (Mcthot and Scarani, 2000). 



I2000D allowing to attach a precise number to the entan- 
glement encoded in a given state. There is a substantial 
bulk of work for bipartite systems, in particular for the 
case of qubits. Many criteria have been proposed to dis- 
tinguish separable from entangled pure states, as for ex- 
ample the Schmidt rank and the von Neumann entropy. 
The success in the bipartite case for qubits asked for ex- 
tensions to the multipartite case, but the situation proved 
to be far more complicated: different classes of entangle- 
ment occur, which are inequivalent not only under de- 
terministic local operations and classical c ommunication , 
but e ven under their stochastic analogue ( Bennett et all 
120011) . 

In the last few years it has become evident that quan- 
tum information may lead to further insight into other 
areas of phy sics as statist ical mechanics and quantum 
field theory dPresk ill. 2000). The attention of the quan- 
tum information community to systems intensively stud- 
ied in condensed matter has stimulated an exciting cross- 
fertilization between the two areas. Methods developed 
in quantum information have proved to be extremely use- 
ful in the analysis of the state of many-body systems. 
At T = many-body systems are most often described 
by complex ground state wave function which contain 
all the correlations that give rise to the various phases 
of matter (superconductivity, ferromagnetism, quantum 
hall systems, . . .). Traditionally many-body systems have 
been studied by looking for example at their response 
to external perturbations, various order parameters and 
excitation spectrum. The study of the ground state of 
many-body systems with methods developed in quantum 
information may unveil new properties. At the same time 
experience built up over the years in condensed matter 
is helping in finding new protocols for quantum compu- 
tation and communication: A quantum computer is a 
many-body system where, differently from "traditional 
ones" , the Hamiltonian can be controlled and manipu- 
lated. 

The amount of work at the interface between sta- 
tistical mechanics and quantum information has grown 
tremendously during the last few years, shining light 
on many different aspects of both subjects. In partic- 
ular, there has been an extensive analysis of entangle- 
ment in quantum cri t ical models (lOsborne and Nielsen! . 
120021 : lOsterloh et a/1 120021: IVidal et all I2003D . Tools 
from quantum information theory also provided impor- 
tant support for numerical methods, as the density ma- 
trix renormalization group or the design of new efficient 
simula tion strategies fo r many - body systems (see for ex- 
ample IVerstraete et all l2004d: I Vidall . I2OO3I |2004|). Spm 
netw orks have been proposed as quantum channels (jBosel . 
|2003[ ) by exploiting the collective dynamics of their low 
lying excitations for transporting quantum information. 

Despite being at its infancy, this new area of research 
has grown so fast that a description of the whole field 
is beyond the scope of a single review. Many interest- 
ing facets of this branch of research will therefore remain 
untouched here. In this review we will only discuss the 
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properties of entanglement in many-body systems. The 
models which will be considered include interacting spin 
networks, itinerant fermions, harmonic and bosonic sys- 
tems. All of them are of paramount interest in condensed 
matter physics. 

This review is organized as follows. In the next Sec- 
tion we give a brief overview on the concepts and mea- 
sures of entanglement, with particular attention to those 
measures that we will use later on. In Section IIIII we 
then proceed with a brief introduction to several mod- 
els of interacting many-body systems which will be sub- 
ject of the review. We will discuss various aspects of 
quantum correlations starting from the pairwise entan- 
glement, Section IIV1 we then proceed with the proper- 
ties of block entropy, Section [V] , and localizable entan- 
glement, Section EH In these three Sections it is espe- 
cially relevant the connection between entanglement and 
quantum phase transitions. The effect of a finite temper- 
ature is considered in Section fVIII The characterization 
of entanglement in many-body systems requires also the 
understanding of multipartite entanglement. This topic 
will be reviewed in Section IVIIII From the point of view 
of quantum information processing, dynamical proper- 
ties of entanglement are important as well. They will be 
addressed in Section IIXI The conclusions, the outlook 
and a very short panorama of what is left out from this 
review are presented in the concluding Section. 



II. MEASURES OF ENTANGLEMENT 

The problem of measuring entanglement is a vast 
and lively field of research in its own. Numerous dif- 
ferent methods have been proposed for its quantifica- 
tion. In this Section we do not attempt to give an 
exhaustive review of the field. Rather we want to in- 
troduce those measures that are largely being used to 
quantify entanglement in many-body systems. Com- 
prehensiv e overviews of entanglement mea s ures can be 
found i n dBengtsson and Zvczkowskil [20061: iBrufj 2002 ; 
Eisertl. 120061: iHorodecki et all l2007t iPlenio and Vedral 



ten as 



1998; Plenio and Virma ni. 2007; Vedral, 2002: iWootters 



20011 ) In this context, we also outline a method of detect- 
ing entanglement, based on entanglement witnesses. 



A. Bipartite entanglement in pure states 

Bipartite entanglement of pure states is conceptu- 
ally well understood, although quantifying it for lo- 
cal dimensi ons higher than two still bears theoretical 
chall enges ( Horodecki et all 120071 IVirmani and Plenkl 
l2000h . A pure bipartite state is not entangled if and only 
if it can be written as a tensor product of pure states of 
the parts. Moreover for every pure bipartite state | i^ab) 
(with the two parts, A and B), two orthonormal bases 
{| ipA,i}} and {| (j)B,j}} exist such that | i/jab) can be writ- 



AB 



^2 &i I i>j 



(i) 



where on are positive coefficients. This decomposition 
is called the Schmidt decomposition and the particular 
basis coincide with the eigenbasis of the corresponding 
reduced density operators Pb/a — tr a/bQ^ab)) = 
Y,i a f | 4>B/A,i) (ipB/A,i | The density operators p A and 
Pb have common spectrum, in particular they are equally 
mixed. Since only product states lead to pure reduced 
density matrices, a measure for their mixedness points 
a way towards quantifying entanglement in this case. 
Given the state \iPab)i we can thus take its Schmidt 
decomposition, Eq.|T]), and use a suitable function of the 
cti to quantify the entanglement. 

An entanglement measure E is fixed uniquely after im- 
posing the following conditions: 1) E is invariant under 
local unitary operations (=> E is indeed a function of the 
«i's only); 2) E is continuous (in a certain sense also 
in the asympt otic limit of infinite copie s of the state; 
see e.g. Ref. (jPlenio and Virmanil . 120071 )): 3) E is ad- 
ditive, when several copies of the system are present: 
E{\i>AB)®\<t>AB)) =E{\i,AB))+E{\(j)AB)). The unique 
measure of entanglement satisfying all the above condi- 
tions is the von Neumann entropy of the reduced density 
matrices 



S(p A ) = S(p B ) 



(2) 



this is just the Shannon entropy of the moduli squared 
of the Schmidt coefficients. In other words: under the 
above regularity conditions, the answer on the question 
of how entangled a bipartite pure state is, is given by the 
von Neumann entropy of (either of) the reduced density 
matrices. The amount of entanglement is generally diffi- 
cult to define once we are away from bipartite states, but 
in several cases we can still gain some insight into many- 
party entanglement if one considers different bipartitions 
of a multipartite system. 

It is worth to notice that a variety of purity measures 
are admissible when the third condition on additivity is 
omitted. In principle, there are infinitely many measures 
for the mixedness of a density matrix; two of them will 
typically lead to a different ordering when the Hilbcrt 
space of the parts has a dimension larger than two. 
In contrast, if we trace out one of two qubits in a pure 
state, the corresponding reduced density matrix pa con- 
tains only a single independent and unitarily invariant 
parameter: its eigenvalue < 1/2. This implies that each 
monotonic function [0, 1/2] i — ^ [0, 1] of this eigenvalue can 
be used as an entanglement measure. Though, also here 
an infinity of different mixedness measures exists, all lead 
to the same ordering of states with respect to their entan- 
glement, and in this sense all are equivalent . A re levant 
example is the (one-) tangle ( Coffman et all 120001) 



Ti [p A ] = 4detpA 



(3) 
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By expressing p A in terms of spin expectation val- 
ues, it follows that Tx[p A ] = 1 - 4((S X ) 2 + (S y ) 2 + 
(S z ) 2 ) where (S a ) = tr A {p A S a ) and S a = \a c \ a a 
{a = x,y,z} being the Pauli matrices, For a pure 
state | tpAs) of two qubits the relation t\ = \ (ip* \ a\ ® 

T2 a pplies , where C is the 
concurrence (|Hill and Wootterl Il997t IWootterd . Il998h 



Wl 2 =: CMar)] 



for pure states of two qubits, a measure of pairwise en- 
tanglement (see next Section), and * indicates the com- 
plex conjugation in the eigenbasis of o~ z . The von Neu- 
mann entropy can be expressed as a function of the 

(one-) tangle S[p A ] = h (j^ (l + \/l - n[/>A]^ where 

h(x) =: —x log 2 x 
tropy. 



(1 — x) log 2 (l — x) is the binary en- 



B. Pairwise qubit entanglement in mixed states 

Subsystems of a many-body (pure) state will gener- 
ally be in a mixed state. In this case different way 
of quantifying entanglement can be introduced. Three 
important representatives are the entanglement cost 
Ec, the distillable entang lement Ed (both defined in 



Ref. (IBennett et a/.l.ll996afr) and th e entanglement of for- 



mation Ep ( Bennett et all , Il996bl ). Whereas Ed and 
Eq are asymptotic limits of multi-copy extraction prob- 
abilities of Bell states and creation from such states, 
the entanglement of formation is the amount of pure 
state entanglement needed to create a single copy of 
the mixe d state. Although recent progress have been 
achieved ( Paz-Silva and Reinal l2007t ), the full additivity 



of the Ep fo r bipartite systems has not been established 
yet (see e.g. (jVidal et aZj . l2002h ). 



The conceptual difficulty behind the calculation of Ep 
lies in the infinite number of possible decompositions of a 
density matrix. Therefore, even knowing how to quantify 
bipartite entanglement in pure states, we cannot simply 
apply this knowledge to mixed states in terms of an av- 
erage over the mixtures of pure state entanglement. The 
problem is that two decompositions of the same density 
matrix usually lead to a different average entanglement. 
Which one do we choose? It turns out that we must take 
the minimum over all possible decompositions, simply 
because if there is a decomposition where the average 
is zero, then this state can be created locally without 
need of any entangled pure state, and therefore Ep = 0. 
The same conclusion can be drawn from the require- 
ment that entanglement must not increase on average 
by means of local operations including classical commu- 
nication (LOCC). 

The entanglement of formation of a state p is therefore 
defined as 



E F (p) := min^pjSXpAj) 



(4) 



mann entropy of the reduced density matrix p A ^ := 
tr b \ il>j) {il>j |- Eq.® is the so called convex roof (also 
the expression convex hull is found in the literature) of 
the entanglement of formation for pure states, and a de- 
composition leading to this convex roof value is called an 
optimal decomposition. 

For systems of two qubits, an analytic expression for 
Ep does exist and it is given by 



cr=± 



o-CHp) 



(5) 



wher e C(p) is the so called concurrence (IWootterd . fl998l 
l200lf ), the convex roof of the pure state concurrence, 
which has been defined in the previous section. Its con- 
vex roof extension is encoded in the positive Hermitean 
matrix R = Jpp^fp = Jp(o y ® a y )p* [a y ® <j y )^Tp , with 
eigenvalues A 2 > • • ■ > A 2 in the following way 



C = maxjAi — X 2 — A3 — A 4 , 0} 



(6) 



As the entanglement of formation is a monotonous func- 
tion of the concurrence, also C itself or its square r 2 - 
called also the 2-tangle - can be used as entanglement 
measures. This is possible due to a curious peculiarity 
of two-qubit systems: namely th at a continuous variety 
of optimal decompositions exist ( Woottersl [T998). The 
concurrence C and the tangle T\ both range from (no 
entanglement) to 1. 

By virtue of ©, the concurrence in a spin- 1/2 chain 
can be computed in terms of up to two-point spin corre- 
lation functions. As an example we consider a case where 
the model has a parity symmetry, it is translational in- 
variant and the Hamiltonian is real; the concurrence in 
this case reads 



dj = 2max{0,C/ i ,C//} 



(7) 



where C], 

= 



9if ~ 9f V] 



.<//; + 9ij 1 - \l (1/4 + 9if 

, „ , , 3 ff - 1/4, with g%* = (S?^ 
M z = (S z ). A state with dominant fidelity of parallel and 
anti-parallel Bell states is characterized by dominant C 1 
and C 11 , respectively. This was shown in ( Fubini et all 
I2006D . where the concurrence was expres sed in terms of 
the fu lly entangled fraction as defined in I Bennett et al\ , 
1996b). A systematic analysis of the relation between 
the concurrence (together with the 3-tangle, see section 
IILFJ and the correlat ion functions has been presented 
in (|Glaser et ^.1 . 120031 ). 

The importance of the tangle and the concur- 
ren ce is due to the m onogamy inequality derived 
in ( Coffman et all I2000T ) for three qubits. This in- 
equali ty has been proved to hold al so for n-qubits sys- 
tem ( Osborne and Verstraetej . [2006). In the case of 
many-qubits (the tangle may depend on the site i that is 
considered) it reads 



Ml and 
"Sf) and 



where the minimum is taken over all realizations of the 
state p A B = Y^i ;Pj\ 1 Pj)(' t Pj\i an, 4 S{pA,j) is the von Neu- 



< n. 



(8) 
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The so called residual tangle n,< — Sj^jCy' 
sure for multipartite entanglement not stored in pairs of 
qubits only. We finally mention that the antilinear form 
of the pure state concurrence was the key for the first ex- 
plicit construction of a convex roof, and hence its exten- 
sion to mixed states ( Hill and Woottersl . \l<M lUhlmannl . 
I2000t IWootterl 119981) . 

Another measure of entangleme nt we mention is th e 
relative entropy of entanglement ( Vedral et all Il997l) . 
It can be applied to any number of qubits in princi- 
ple (or any dimension of the local Hilbert space). It 
is formally defined as E(a) := mia p ^T> 5"(cr| , where 
S^crllp) = tr a [In a — In p] is the quantum relative en- 
tropy. This relative entropy of entanglement quantifies 
the entanglement in a by its distance from the set T> of 
separable states (since T> is compact, the minimum is as- 
sumed always). The main difficulty in computing this 
measure is to find the disentangled state closest to p. 
This is in general an open problem, even for two qubits. 
In the presence of certain symmetries - which is the case 
for e.g. eigenstates of certain models - an analytical ac- 
cess is possible. In these cases, the relative entropy of 
entanglement becomes a very useful tool. The relative en- 
tropy reduces to the entanglement entropy in the case of 
pure bi-partite states; this also means that its convex roof 
extension coincides with the entanglement of formation , 
and i s readily deduced from the concurrence ( Woottersl . 
Il998h . 

We close this summary on the pairwise entanglement 
by commenting on the notion on the quantum mutual in- 
formation. Groisman et al quantified the work necessary 
to er ase the total correlation s existing in a bipartite sys- 
tem ()Groisman et ali , [2005) . The entanglement can be 
erased by a suitable random ensemble of unitary transfor- 
mations acting on one of the parts, but a certain amount 
of classical correlation among the two partners may sur- 
vive. The work necessary to erase all correlations is given 
by the quantum mutual information 



x-direction on one of the sites is transformed into a Bell 
state. 

The conce pt of localizable entanglement has been in - 
troduced in ijpppp et~al\ . 120051 : IVerstraete et al. I. l2004ah . 
It is defined as the maximal amount of entanglement that 
can be localized, on average, by doing local measurements 
in the rest of the system. In the case of N parties, the 
possible outcomes of the measurements on the remaining 
N — 2 particles are pure states \ip s ) with corresponding- 
probabilities p s . The localizable entanglement Ei oc on 
the sites i and j is defined as the maximum of the aver- 
age entanglement over all possible outcome states | ips)^ 



Eioc(i,j) 



(10) 



where £ is the set of all possible outcomes (p s , |V>s)) of 
the measurements, and E represents the chosen measure 
of entanglement of a pure state of two qubits (e.g. the 
concurrence). Although very difficult to compute, lower 
and upper bounds have been found which allow to deduce 
a number of consequences for this quantity. 

An upper bound to the localizable entanglement is 
given by the entanglement of assistance ( Laustsen et ali , 
2003) obtained from localizable entanglement when also 
global and joint measurements were allowed on the N — 2 
spins . A lower bound of the localiza ble entanglement 
comes from the following theorem ( Verstraete et ali 
2004a): Given a (pure or mixed) state of N qubits with 
reduced correlations Qlf = (SfSj) — (S")(Sj) between 
the spins i and j and directions a and (3 then there al- 
ways exist directions in which one can measure the other 
spins such that this correlation do not decrease, on av- 
erage. It then follows that a lower bound to localizable 
entanglement is fixed by the maximal correlation func- 
tion between the two parties (one of the various spin-spin 
correlation functions Q"^) 2 . 



Iab = S(pa) + S(p B ) - S(pab) 



(9) 



D. Entanglement witnesses 



C. Localizable entanglement 

A different approach to entanglement in many-body 
systems arises from the quest to swap or transmute dif- 
ferent types of multipartite entanglement into pairwise 
entanglement between two parties by means of general- 
ized measures on the rest of the system. In a system 
of interacting spins on a lattice one could then try to 
maximize the entanglement between two spins (at posi- 
tions i and j) by performing measurements on all the 
others. The system is then partitioned in three regions: 
the sites i, j and the rest of the lattice. This concen- 
trated pairwise entanglement can then be used e.g. for 
quantum information processing. A standard example is 
that the three qubit Greenberger-Horne-Zeilinger (GHZ) 
state (l/\/2)(| 000) + | 111)) after a projective measure in 



It is important to realize that not just the quantifica- 
tion of many-party entanglement is a difficult task; it is 
an open problem to tell in general, whether a state of n 
parties is separable or not. It is therefore of great value 
to have a tool that is able to merely certify if a certain 
state is entangled. An entanglement witness W is a Her- 
mitean operator which is able to detect entanglement in 
a state. The basic idea is that the expectation value of 
the witness W for the state p under consideration exceeds 



2 It ha s been argued recently {Gour, 2006; Gour and Sp ekkensl . 
2006) that in order to extend the entanglement of assistance and 
the l ocalizable en tanglement to being an entanglement mono- 
tone llVidalll2000l) one should admit also local operations includ- 
ing classical communication on the extracted two spins, this was 
named entanglement of collaboration. 
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certain bounds only when p is entangled. An expectation 
value of W within this bound however does not guarantee 
that the state is separable. Nonetheless, this is a very ap- 
pealing method also from an experimental point of view, 
since it is sometimes possible to relate the presence of the 
entanglement to the measurement of few observables. 

Simple geometric ideas help to explain the witness 
operator W at work. Let T be the set of all density 
matrices and let £ and V be the subsets of entangled 
and separable states, respectively. The convexity of T> 
is a key property for witnessing entanglement The en- 
tanglement witness is then an operator defining a hyper- 
plane which separates a given entangled state from the 
set of separable states. The main scope of this geo- 
metri c approach is then to o ptimize the witness oper- 
ator ( Lewenstein et all l2000ft or to replace the hyper- 
plane by a c urved manifol d, tangent to the set of sepa- 
rable states ( Guhnel . l 2004ft ffor other geometric aspects 
of entanglement see (Bengtsson and Zvczkowskil 120061: 
lKlvachkoL[200llLeinaas et a/.l . l2006h ). We have the free- 
dom to choose W such that tr (pdW) < for all disen- 
tangled states po £ T>. Then, tr (pW) > implies that 
p is entangled. A caveat is that the concept of a wit- 
ness i^sjiaMirvajiairtjj^ operations (see 
e.g. ICavalcanti and Terra-Cunhal 12005ft . 

Entanglement witnesses are a special case of a more 
general concept, namely that of positive maps. These 
are injective superoperators on the subset of positive op- 
erators. When we now think of superoperators that act 
non-trivially only on part of the system (on operators 
that act non trivially only on a sub-Hilbcrt space), then 
we may ask the question whether a positive map on the 
subspace is also positive when acting on the whole space. 
Maps that remain positive also on the extended space 
are called completely positive maps. The Hermitean time 
evolution of a density matrix is an example for a com- 
pletely positive map. Positive but not completely positive 
maps are important for entanglement theory. There is a 
remarkable isomo rphism between posi tive maps and Her- 
mitean operators (iJamiolkowskil 1 1972ft. This can be used 
to prove a key theorem (Horodc cki et all Il996f ): A state 
Pab is entangled if and only if a positive map A exists 
(not completely positive) such that (1a <B> ^-b)pab < 0. 
For a two dimensional local Hilbert space the situation 
simplifies considerably in that any positive map P can 
be written as P = CP X + CP 2 T B , where CPx and CP 2 
are completely positive maps and Tb is a transposition 
operation on subspace B. This decomposition tells that 
for a system of two qubits the lack of complete positivity 
in a positive map is due to a partial transposition. This 
partial transposition clearly leads to a positive operator 
if the state is a tensor product of the parts. In fact, also 
the opposite is true: a state of two qubits pab is separa- 
ble if and only if p^B — that is, its partial transposition 
is positive. This is very simple to test and it is known 
as the Peres -Horodecki criterion ( Horodecki et aZ.I . Il996t 
iPerest Il996ft . The properties of entangled states under 
partial transposition lead to a measure of entanglement 



known as the negativity. The negativity Nab of a bipar- 
tite state is defined as the absolute value of the sum of the 
negative eigenvalues of Pab- The logarithmic negativity 
is then defined as 



E N =log 2 2(2A^ B + l). 



(11) 



For bipartite states of two qubits, p^% ha s at most 
one negative eigenvalue ( Sanpera et all 1l998). For gen- 
eral multipartite and higher local dimension this is only 
a sufficient condition for the presence of entanglement. 
There exist entangled states with a posi tive partial trans- 
pose known as bound e ntangled states ( Acin et aZl . 120011 : 
iHorodecki et aUI 1998ft . 



E. Multipartite entanglement measures 

Both the classification of entanglement and its 
quantification are at a preliminary st age even for 



Briand et al 



distinguishable partic les (se e however 

2003112004 iDiir et all l2000l iLuaue and Thibonl . [ 2005 : 
Mandilara et all 120061: iMivake and Wadatll IgOOa 



Osterloh and SiewertU2005l l200rj IVerstraete et qL| . 120021 
and references therein). We restrict ourselves to those 
approaches which have been applied so far for the study 
of condensed matter systems discussed in this review. 
It has already been mentioned that several quantities 
are useful as indicators for multipartite entanglement 
when the whole system is in a pure state. The entropy 
of entanglement is an example for such a quantity 
and several works use multipartite measures con- 
structed from and rela ted to it 
2005 



2003: Coffman et al 



(see e.g. 



Barnum et 



de Oliveira et al 



Love et all 120061: iMever and Wallachl 1200a IScott 



2006a; 



20041 ). These measures are of 'collective' nature - in 
contrast to 'selective' measures - in the sense that they 
give indication on a global correlation without discerning 
among the different entanglement classes encoded in the 
state of the system. 

The geometric measure of entanglement quantifies 
the entanglement of a pure state through the mini- 
mal di stance of the state from the set of pure pro duct 
states (|Vedral et all 1 19971: IWei and Goldbartl . 12003ft 



E 9 m 



log 2 max 



(12) 



where the maximu m is on all product states $. As dis- 
cussed in detail in ( We i and Goldbartl . I2003T ) , the previ- 
ous definition is an entanglement monotone if the convex- 
roof extension to mixed states is taken. It is zero for sep- 
arable states and rises up to unity for e.g. the maximally 
entangled n-particle GHZ states. The difficult task in its 
evaluation is the maximization over all possible separable 
states and of course the convex roof extension to mixed 
states. Despite these complications, a clever use of the 
symmetries of the problem renders this task accessible 
by substantially reducing the number of parameters (see 
Section En]}. 
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Another example for the collective measures of mul- 
tipartite entanglement as mentioned in the beginning 
of this section are the measures i ntrodu ced by Meyer 
and Wallach dMever and Wallacbl l2002h and by Bar- 
num et al rtBarnum et all I2004L l2003h . In the case 



of qubit system the Q-measure of Meyer and Wallach 



in dCoffman et all 2000f)) of the state ( 


Barnum et all 


2004: Brennen, 


20031: Mever and Wallach, 2002) 



TV 



E, 



gl 



N 



E^ 2 



3=1 



The notion of generalized entanglement introduced 
in ( Barnum et all , [20041 l2003h relaxes the typically cho- 
sen partition into local subsystems in real space. The 
generalized entanglement measure used by Barnum et al. 
is the purity relative to a distinguished Lie algebra A of 
observables. For the state lib) it is defined as 



PA = r ri{[p A wm 2 } 



(14) 



where Va is the projection map p — ► Va(p)- If 
the set of observables is defined by the operator basis 
{Ai, A 2 , . . . , A L } then P A = Yh=\{ A i) 2 fr° m which the 
reduction to Eq. (Ti"3"]) in the case of all local observables is 
evident. This conceptually corresponds to a redefinition 
of locality as induced by the distinguished observable set, 
beyond the archetype of partition in the real space. It 
defines an observer dependent concept of entanglement 
adapted to e.g. experimentally accessible or physically 
relevant observables. In this case, the generalized entan- 
glement coincides with the global entanglement of Meyer 
and Wallach. 

Another approach pursued is the generalization of 
the concurrence. For the quantification of pairwise en- 
tanglement in higher dimensional local Hilbert spaces, 
the c oncept of concur r ence vectors has been fo rmu- 
lated ([Audenaert et aU feOOll iBadziag et all l2002f) be- 
sides the I-concurrence ( Rungta et all 2001 ). A concur- 
rence vect or was also proposed for multipartite systems 
of qubits ( Akhtarshenasl [20051 ). It consists in applying 
the pure state concurrence formula to a mixed two-site 
reduced density matrix. It coincides with the true con- 
currence if and only if the eigenbasis of the density ma- 
trices provide optimal decompositions. 

The n-tangle is a straightforward extension of the con- 
currence to multipartit e states as the overlap o f the s tate 
with its time-reversed (|Wong and Christensenl . 1200 ll ) . It 
vanishes identically for an odd number of qubits, but 
an entanglement monotone is obtained for an even num- 
ber of qubits. It detects products of even-site entangled 
states in addition to certain genuine multipartite entan- 
gled states: it detects the multipartite GHZ or cat state, 
but not for example the four qubit cluster state. 

Three classes of states inequivalent un- 
der SLOCC (Stochastic LOCC) exist for four 



qubits (|Osterloh and Siewertl . 120051 l2006h . Repre- 
sentatives are the GHZ state, the celebrated cluster 
state and a third state, which is also measured by the 
4-qubit Hyperdeterminant. Class selective measures 
are constructed from two basic elements, namely the 
operator a y employed for the concurrence, and the 



operator <Tu ■ a* 1 



1 • 1 



o~ z ■ g z where 



the • is a tensor product indicating that the two 
operators are acting on different copies of the same 
qubit. Both are invariant under sl(2, C) operations 
on the qubit. The 3-tangle is then expressed as 



(13) T 3 [lp] = (^1* | • (tp* | O fj 



The multilinearity, however, makes it problematic to 
employ th e procedure of convex roof const ruction pre- 
sented in (jUhlmanrJ . l200d : IWoottersL Il998| ) for general 
mixtures. 

F inally we ment ion the approach pursued 
in (iGiihne et al. |, I2005D (see also [Sh arma and Sharmal . 
I200fih where different bounds on the average energy 
of a given system were obtained for different types of 
n-particle quantum correlated states. A violation of 
these bounds then implies the presence of multipar- 
tite entanglement in the system. The starting point 
of Giihne et al. is the notion of n- separability and 
k-producibility which admit to discriminate particular 
types of n-particle correlations present in the system. 
A pure state | tp) of a quantum systems of N parties is 
said to be n-separable if it is possible to find a partition 
of the system for which | ip) = |0i)|02) ■ ■ ■ \4>n)- A pure 
state | ip) can be produced by k-party entanglement (i.e. 
it is k-producible) if we can write | V) = |0i)|$2) • 1 1 \<j>m) 
where the \<j>i) are states of maximally k parties; by 
definition m > N/k. It implies that it is sufficient to 
generate specific k-party entanglement to construct the 
desired state. Both these indicators for multipartite 
entanglement are collective, since they are based on the 
factorizability of a given many particle state into smaller 
parts, k-separability and -producibility both do not 
distinguish between different k-particle entanglement 
classes (as e.g. the k -particle W-state s and different 
k-particle graph states ( Hein et~cd\ , I2004T ). like the GHZ 
state). 



F. Indistinguishable particles 

For indistinguishable particles the wave function is 
(anti-) symmetrized and therefore the definition of entan- 
gled states as given in the previous Section does not ap- 
ply. In particular, it does not make sense to consider each 
individual particle as parts of the partition of the sys- 
tem. Having agreed upon a definition of entanglement, 
concepts as entanglemen t cost or distillation rema i n per- 
fectly valid. Followi ng ( Ghirardi and Marinattol . 120031 : 
iGhirardi et all 120021) one can address the problem of 
defining entanglement in an ensemble of indistinguish- 
able particles by seeing if one can attribute to each of the 
subsystems a complete set of measurable properties, e.g. 
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momenta for free pointless particles. Quantum states 
satisfying the above requirement represent the separable 
states for indistinguishable particles. 

There is another crucial difference between the en- 
tanglement of (indistinguishable) spin-1/2 particles and 
that of qubits. Let us therefore consider two fermions 
on two sites. Whereas the Hilbert space H s of a two- 
site spin lattice has dimension dim7i s = 4, the Hilbert 
space Ttf for two fermions on the same lattice has di- 
mension dim7iy = 6. This is due to the possibility that 
both fermions, with opposite spins, can be located at the 
same lattice site. When choosing the following number- 
ing of the states 1 1) = }\ | 0) =: c\ 1 1 0), | 2) = f\ | 0) =: 



J 



0), I 3) = f\ | 0) =: 4, t | 0), 1 4) = fl | 0) =: | 0) 

and the definition | i, j) = fjfj | 0), there are Bell states 
analogous to those occurring for distinguishable parti- 
cles (| 1,3) ± | 2, 4))/V2 and (| 1,4) ± |2,3))/\/2. There 
are however new entangled states, as (| 1, 2) ± | 3, 4))/v / 2, 
where both fermions take the same position. The local 
Hilbert space is made of four states labelled by the oc- 
cupation number and the spin, if singly occupied. The 
site- entanglement of indistinguishable particles is then 
defined as the entanglement of the corresponding Fock 
states. It can be measured e.g. by the local von Neu- 
mann entropy. This quantity is the analogue to the one- 
tangle for qubits, but the local Hilbert space dimension 
is 4 due to the possibility of having empty and dou- 
bly o ccupied sites. Also th e quantum mutual informa- 
tion ( Groisman et all 120051 ). see Eq.©, can be defined 
in this way, quantifying the total amount (classical and 
quantum) of correlations stored in a given state of a sec- 
ond quantized system. 

Although from a mathematical point of view the en- 
tanglement of indistinguishable particles can be quanti- 
fied, the major part of the literature on second quan- 
tized systems that we discuss in this review considers the 
site-entanglcment described above or the entanglement 
of degrees of freedom, singled out from a suitable set of 
local quantum numbers (e.g. the spin of the particle at 
site i). In both cases, entanglement measures for distin- 
guishable particles (see Sections IIV.C.1I and IV. F[) can be 
used. With this respect, this Section has a different scope 
than the others on the quantification of entanglement; al- 
though most of the discussion which follows will not be 
used later on, we believe that it will be of interest for 
further studies of entanglement in itinerant many-body 
systems. 



1. Two Fermion entanglement 

Due to the antisymmetry under particle exchange, 
there is no Schmidt decomposition for Fermions. Nev- 
ertheless, a Fermionic analogue to the Schmidt rank, 
which classifies entanglement in bipartite systems of dis- 
tinguishable particles does exist: the so called Slater 
rank. A generic state of two-electrons on two lattice sites 



can be written as | lo) := J2 i . =1 LOi.j | i, j) where w is a 
4x4 matrix which can be assumed antisymmetric and 
normalized as tr lo^lo = i. Since here the local entities 
whose entanglement shall be studied, are the particles, 
unitary transformations act on the 4-dimensional single 
particle Hilbert space. Due to the indistinguishability 
of the particles, the transformation must be the same 
for each of the particles. Given a unitary transforma- 
tion U G SU(4) such that fj := Ujkf k , the transformed 
state is given by | lo') where lo' :— UuU T . The above 
unitary transformation preserves the antisymmetry of lo 
and can transform every pure state of two spin-1/2 par- 
ticles on two sites into a state corresponding to the nor- 
mal form of lo. In fact, every two-particle state within a 
D-dimensional single particle Hilbert space can be trans- 
formed into the normal form lo s = diag{Zi, . . . , Z r , Zq} 



where Z 



j 



j(j y and {Zo)ij = for i, j e {1, . . . , D 



2r}. In the previous expression r is then called the 
Slater rank of the pure Fe rmion state ( Eckert et all [20021: 
ISchliemann et all l2001al fbl) . A pure Fermion state is en- 
tangled if and only if its Slater rank is larger than 1. It is 
important to notice that the above concept of entangle- 
ment only depends on the dimension of the Hilbert space 
accessible to each of the particles (this includes indistin- 
guishable particles on a single D-level system). 

For electrons on an L-site lattice the "local" Hilbert 
space dimension is 2L, and the question, whether a pure 
state living in a 2L-dimensional single particle Hilbert 
space has full Slater r ank, can be answered by consider- 
ing t he Pfaffian of lo ( Caianello and Fubinil . Il952t iMuirl . 
1960) 



7reS< 2=1 



l),w(2j) 



(15) 



which is non-zero only if lo has full Slater rank L. In 
the above definition denotes those elements ir of the 
symmetric group S2L with ordered pairs, i.e. 7r(2m— 1) < 
7t(2to) for all m < L and ir(2k — 1) < 7r(2m — 1) for 
k < m. Notice that relaxing the restriction to just 
leads to a combinatorial factor of 2 L L\ by virtue of the 
antisymmetry of lo and hence can we write 



P fH 



1 



2L 

j'lj— ,22£=1 



■3X, — ,l2L, 



. LOj 



where e n, '" ,nL is the fully antisymmetric tensor with 



_1,2,...,2L 



= 1. There is a simple relation between the 
Pfaffian and the determinant of an antisymmetric even- 
dimensional matrix: pf[w] 2 = det[w]. 

For the simplest case of two spin-1/2 Fermions on 
two lattice sites the Pfaffian reads pf[w] = ^1,2^3,4 — 
wi,3o;2,4 + LJi,4t02,3- Normalized in order to range in the 
interval [0, 1] thi s has been called t he Fermionic con- 
currencc C[\lo)} (jEckert et all 120021 : ISchliemann et all 
l2001allbh 



c[| w )] = |<a| w )| = 8|pf[ w ]| 



(17) 
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1 r-i jkl , ,* 



has been termed the dual to ui. 



2 1 - "k,l 

u>) — : P | a;) is t he analog to the conjugated state 



where u 
Then, 

in dHill and Woottersl Il997t lUhlmannl . l2000l : IWootterl . 
Il998h leading to the concurrence for qubits. It is im- 
portant to notice that the Pfaffian in Eq. ([T5|) is invari- 
ant under the complexification of su(2L), since it is the 
expectation value of an antilinear operator, namely the 
conjugation T> for the state \uj). Since this invariant 
is a bilinear expression in th e state coefficien ts, its con- 
vex roof is readily obtained (|Uhlmannl . [2000) by means 
of the positive eigenvalues Af of the 6x6 matrix R — 
^fpDpD^/p. The conjugation 2?, expressed in the basis 
{| 1,2), 1 1,3), 1 1,4), I 2, 3), I 2, 4), I 3, 4)} takes the form 
T>o C, where C is the complex conjugation and the only 
non-zero elements of Vq are (£>o)i6 = 1> (^0)25 = — 1 3 
(2? ) 3 4 = 1, (^0)43 = 1, (^0)52 = -1, and (2?o) 8 i = 1- 
Notice that the center part of this matrix is precisely 
(jy ® <Jy and indeed corresponds to the Hilbert space of 
two qubits. The remaining part of the Hilbert space gives 
rise to an entanglement of different values for the occu- 
pation number. This type o f entanglement has been re- 
ferred to as the fluffy bunny (jVerstraete and Cirad . l2003t 
IWiseman et all 120031 ) in the literature. 

For a single particle Hilbert space with dimension 
larger than 4 one encounters similar complications as for 
two distinguishable particles on a bipartite lattice and 
local Hilbert space dimension larger than 2, i.e. for two 
qudits. This is because different classes of entanglement 
occur, which are characterized by different Slater rank as 
opposed to their classification by different Schmidt rank 
for distinguishable particles. The Sla ter rank can be ob- 
tained by looking at Pfaffian minors <|Muirl . ll960l) : if the 
Slater rank is r, all Pfaffian minors of dimension larger 
than 2r are identically zero. 



2. Multipartite Entanglement for Fermions 

For indistinguishable particles the only classification 
available up to now is to check whether or not a pure 
state has Slater ra nk one. Eckert et a l. formulated two 
recursive lemmata ( Eckert et all [20021 ) that can be sum- 
marized as follows: let an iV-electron state be contracted 
with N — 2 arbitrary single electron states encoded in 



the vectors a-? as ai/J | 0) (j — 1, . . . , N — 2 and sum 
convention) to a two-electron state. Then the Pfaffian of 
the two-electron state is zero if and only if the original 
state (and hence all intermediate states in a successive 
contraction) has Slater rank one. This means that all 
4-dimensional Pfaffian minors of to are zero. 

Instead of the Pfaffian of u>, also the single-particle 
reduced density matrix can be considered, and its 
von Neumann entropy as a measure for the quan- 
tum entanglement h a s been analyzed in ( Li et all l200lt 
iPaskauskas and Youl . l200ll ). It is important to remind 
that for distinguishable particles the local reduced den- 
sity matrix has rank one if and only if the original state 
were a product. This is no longer true for indistinguish- 



able particles. For an iV-particle pure state with Slater 
rank one the rank of the single-particle reduced density 
matrix coincides with the number of particles, N. A sub- 
tlety is that a measure of entanglement is obtained after 
subtraction of the constant value of the von Neumann 
entropy of a disentangled state. This must be taken into 
account also for the extension of the measure to mixed 
states. 



3. "Entanglement of particles" 

Entanglement in the presence of super selection 
rules (SSR) induced b y particle conservatio n has 
been discussed in Refs. (iBartlett and Wiseman! . 120031 : 
ISchuch et all I2003L [200l IWiseman and Vaccarol [2003). 
The mai n difference in the concept of entanglement of 
particles ( Wiseman and Vaccarol l2003h from the entan- 
glement of indistinguishable particles as described in the 
preceding section (but also to that obtained from the re- 
duced density matrix of e.g. spin degrees of freedom of 
indistinguishable particles) consists in the projection of 
the Hilbert space onto a subspace of fixed particle num- 
bers for either part of a bipartition of the system. The 
bipartition is typically chosen to be space-like, as moti- 
vated from experimentalists or detectors sitting at dis- 
tinct positions. E.g. two experimentalists, in order to 
detect the entanglement between two indistinguishable 
particles, must have one particle each in their laboratory. 

This difference induced by particle number supersclcc- 
tion is very subtle and shows up if multiple occupan- 
cies occur at single sites for Fermions with some inner 
degrees of freedom, as the spin. Their contribution is 
finite for finite discrete lattices and will generally scale 
to zero in the thermodynamic limit with vanishing lat- 
tice spacing. Therefore both concepts of spin entan- 
glement of two distant particles coincide in this limit. 
Significant differences are to be expected only for finite 
non-dilute systems. It must be noted that the same re- 
strictions imposed by SSR which change considerably the 
concept of entanglement quantitatively and qualitatively, 
on the other hand enable otherwise impossible protocols 
of qu antum information processing ( Schuch et all 120031 
2004) which are based on variances about the observable 
fixed by superselection. 

Wiseman and Vaccaro project an TV-particle state 
I V>jv) onto all possible subspaces, where the two parties 
have a well defined numb er (714, = N — ua) of par - 
ticles in their laboratory ( Wiseman and Vaccarol [20031 . 



Let j ^[n^]) be the respective projection, and let p nA be 
the weight ( if)[nA\ \ ip[nA] ) / ( ipN | ^jv ) of this projection. 
Then the entanglement of particles E p is defined as 



E p [\^ n )]=Y,PnE M [^[ 



TlA\ 



(18) 



where Em is some measure of entanglement for distin- 
guishable particles. Although this certainly represents 
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a definition of entanglement appealing for experimen- 
tal issues, it is sensitive only to situations, where e.g. 
the two initially indistinguishable particles eventually are 
separated and can be examined one-by-one by Alice and 
Bo b. Consequently, "local operat ions" have been defined 
in ( Wiseman and Vaccarol 12003) as those performed by 
Alice and Bob in their laboratory after having measured 
the number of particles 3 . 

Verstraete and Cirac pointed out that the presence of 
SSR gives rise to a new resource which has to be quan- 
tified. They have proposed to replace the quantity E p 
with the SSR- entanglement of formation. This is defined 



E 



(SSR), 
f I 



mm 



E 



PnE M [i>r, 



where the minimization is performed over all those de- 
composition of the density matrix where the \^) n are 
eigenstates of the total number of particles ( Schuch et all 
2003. l2004r i. 



4. Entanglement for Bosons 

The quantification and classification of boson entan- 
glement is very close in spirit to that of Fermions as de- 
scribed in Section III. F. II In the bosonic case the matrix 
ui introduced in the previous section is symmetric un- 
der permutations of the particle numbers. Consequently, 
for any two-particle state of indistinguishable bosons, w 
can be diagonalized by means of unitary transformations 
of the single particle basi s. This leads to th e Schmidt 
decomposition for bosons (jEckert et al. L l2002h . An curi- 
ous feature distinguishing this case from the entangle- 
ment measures of distinguishable particles is that the 
Schmidt decomposition is not unique. In fact, any two 
equal Schmidt coefficients admit for a unitary transfor- 
mation of the two corresponding basis states, such that 
the superposition of the two doubly occupied states can 



3 As a potential difference between the entanglement of photons 
as opposed to that of massive bosonic particles, it has been 
claimed that certain superselection rules may hold for massive 
particles only. One such claim is that we would in practice 
not be able to build coherent superpositions of states contain- 
ing a di fferent number of ma ssive particles (for a recent discus- 
sion see iBartlett et oil . 1200^ 1. This superselection rule would, 
for instance, prohibit creating a superposition of a Hydrogen 
atom and a Hydrogen molecule. However, the origin and va- 
lidity of any superselection rule remains a very much debated 
subject. The arguments pro superselection rules usually involve 
some symmetry considerations, or some decoherence mechanism. 
On the other hand, it turns out that if we allow most general 
operations in quantum mechanics, we no longer encounter any 
superselection restriction s. Recent work [Dowling et all 2006a ; 
iTerra Cunha et al.. 2006 shows that it should be possible to co- 
herently superpose massiv e particles and to observe a violation 
of certain Bell inequalities Terra Cunha et al., 2006 also for this 



be writ ten as a symmetrized s t ate o f two ort h ogona l 
states (jGhirardi and MarinattoL 120051 : iLi et all l200lh . 
This is the reason why it is not directly the Schmidt rank, 
but rather the reduced Schmidt rank - obtained after 
having removed all double degeneracies of the Schmidt 
decomposition - that determines whether or not a state 
is entangled. This non-uniqueness of the Schmidt rank is 
also responsible for the ambiguity of the von Neumann 
entropy or other purity measures of the single particle 
reduced density matrix as an entangle ment measure for 
Bosons ( Ghirardi and MarinattoL l2005h . 

With Zi being the Schmidt coefficients with degener- 
acy <7i, the reduced Schmidt rank is at most # + 2 {3-}, 
where {.} denotes the non- integer part. As a conse- 
quence, a Schmidt rank larger than two implies the pres- 
ence of entanglement. Schmidt rank 2 with degenerate 
Schmidt coefficients can be written as a symmetrized 
product of orthogonal states and consequently is disen- 
tangled (|Ghirardi and MarinattoL l2005f ) . This features is 
also present in the TV-boson case, where in presence of 
up to AT- fold degenerate Schmidt coefficients the corre- 
sponding state can be rewritten as a symmetrization of 
a product. 

For bipartite systems uj has full Schmidt rank if det lu =/= 
0. A Schmidt rank 1 can be verified by the same con- 
traction technique described for the Fermion case in the 
previous section, where the Pfaffian must be replaced by 
the determinant. This applies to both th e bipartite and 
the multipartite case (Eck ert et aLl . l2002D . 



G. Entanglement in harmonic systems 

In this section we concentrate on the entan- 
glement between d istinct modes of har monic 
oscillators (see dAdesso and Illuminatll . l2007t 
iBraunstein and van LoockL 2005 ) for recent reviews 
on the subject). The entanglement in this case is 
termed as continuous variable entanglement in the 
literature (to be distinguished from the entanglement of 
indistinguishable bosonic particles; see Section |II.F|) . 

Dealing with higher dimensional space of the local de- 
grees of freedom generally involves complications which 
are not tamable within the current knowledge about en- 
tanglement. The Peres-Horodecki criterion, just to men- 
tion an important example, is not sufficient already for 
two three-level systems, 3x3. The situation simplifies 
if only so called Gaussian states of the harmonic oscil- 
lator modes are considered. This restriction makes the 
infinite dimensional case even conceptually simpler than 
the finite dimensional counterparts. In order to explain 
what Gaussian states are, we introduce th e Wigner dis- 
tribution function W(p, q) (|Wignerl . Il932h . For a single 
degree of freedom it is defined from the density operator 
p as 

W{r, p);=— dr' (r + r'\p\r- r') e% pr ' , (19) 

7™ 
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where r and p are conjugate position and momentum 
variables of the degree of freedom. The connection be- 
tween bosonic operators a, a* and phase space operators 
f, p is a = (f + ip)/y/2, = (ay = (r — ip)/y/2. More 
degrees of freedom are taken into account in a straight 
forward manner. A (mixed) state p is then called Gaus- 
sian when its Wigner distribution function is Gaussian. 
Examples for such states are coherent pure states | a) , 
a | a) = a \ a) with a g C, and arbitrary mixtures of co- 
herent states p = J d 2 a P(a) | a) (a |, determined by the 
so called P- distribution P(z). Such states are termed 
classical if the Wigne r funct i on an d the P-distribution 
are non- negative (see ( SimonlfeoOOl )). 

The key quality of Gaussian states is that they are 
completely classified by second moments, which are en- 
coded in the symmetric so called (co-)variance matrix 
with the uncertainties of the phase space coordinates as 
entries. For two bosonic modes the phase space is four- 
dimensional and the covariancc matrix V is defined as 



V af . 



{A| Q) A&} 



d 4 C A^ a A^W({^}) , (20) 



where the curly brackets on the left hand side indicate the 
anti-commutator. The components of £ Q , a — 1,...,4 
are (fi,fii, rz,p%) and A£ Q := £ Q — (£ a ); the averages 
(.) is taken with respect to the given two-mode density 
matrix p, or equivalently, using the Wigner distribution 
of p. Then, the canonical commutation relations assume 
the compact form [£ Q ,£/3] = iftQa,/3 with fl — ia y 1. 
When expressed in terms of V, the Heisenberg uncer- 
tainty relation can be invoked in invariant form with re- 
spect to canonical t ransformations as det V ~>\ (see e.g. 
ISimon et a/.l . ll994h . The set of the real linear canonical 
transformation generates the symplectic group Sp(2n,\R) 
that plays an important role in the theory. Being a sym- 
plectic matrix, V can be brought in its diagonal form V n 
by means of symplectic transformations. The elements 
on the diagonal are then called the symplectic eigenvalues 
of V. An analysis of V n tt has unveiled an even more pow- 
erful invariant form of the Heisenberg uncertainty prin- 
ciple, V + ^il > 0, where the positive semi-definiteness 
means that all symplectic eigenvalues are non-negative. 
The uncertainty relation can hence be cast directly in 
terms of the symplectic eigenvalues of the covariance ma- 
trix V, which are the absolute values of the eigenvalues 
of -iQV. 

Some of the aspects of the harmonic systems can be 
disclosed by recognizing that the Gaussian structure of 
the bosonic states can be thought as a certain limit of 
the algebraic structure of the qubits in the sense that 
Sp(2,IR) ~ SL(2, C). The latter is the invariance group 
relevant f or qubit en t anglement classification and quan- 
tification four et all |2000t lOsterloh and Siewertl l2005t 
IVerstraete et all 120031 ) ~ 

We now introduce the notion of bipartite entangle- 
ment for Gaussian states. In complete analogy to the 
finite-dimensional case, a state is termed separable if it 
is a mixture of product states. In particular, all classical 



states, i.e. 

P = J d 2 Zl d 2 z 2 P( Zl ,z 2 ) | z x ) (z x \®\ z 2 ) (z 2 | (21) 

with positive P (z\,Z2) are se parable. 

It was Simon ( Simonl . l2000[) that first proved the Peres- 
Horodecki Positive Partial Transpose criterion being nec- 
essary and sufficient for entanglement of two harmonic 
oscillator modes, again in complete analogy to a system 
of two qubits. The effect of the transposition of the den- 
sity matrix is a sign change in the momentum variables of 
the Wigner function (fH?)) . Consequently, a partial trans- 
position induces a sign change of those momenta in the 
phase space vector, where the transposition should act 
on. For an entangled state, the partial transposition V 
of its covariance matrix V might then have symplectic 
eigenvalues smaller than ti/2. This can be detected by 
the logarithmic negativity as defined from the symplec- 
tic (doubly degenerate) eigenvalues {ci ; i = 1 . . . n} of 
V/h (IVidal and Werner! . I2002D 



En(V) = -^log 2 (2c i ) 



(22) 



;.=i 



These important results paved the way towards a system- 
atic analysis of multipartite systems of distinguishable 
bosonic modes. 



III. MODEL SYSTEMS 

This section is devoted to the basic properties of the 
model systems that will be analyzed in the rest of the re- 
view (in several cases we concentrate on one-dimensional 
systems). 



A. Spin models 



Interacting s pin models ( Auerbachl Il998t 
ISchollwock et all |2004 ) provide a paradigm to de- 
scribe a wide range of many-body systems. They 
account for the effective interactions in a variety of very 
different physica l contexts rang i ng from high energy t o 
nuclear physics (jBelitskv et all 12004 iPolvakovl 1 19771 ). 
In condensed matter beside d escribing the pro p erties 
of magnetic compounds (see iMatsumoto et al. , 2004 
for a recent survey), they capture several aspects 
of high-temperature superconductors, quantum Hall 
systems, heavy fermions, just to mention few impor- 
tant examples. Hamiltonians for interacting spins 
can b e realized artificially in Josep hson junctions ar- 
rays (Fa zio and van der Zantl l200ll ) or with neutral 
atoms loa d ed in optical lattices (iDuan et all 120031 : 
iJane et all 120031 : iPorras and Cirad . 1200J). Interacting 
spins are paradigm systems fo r qua ntum information 
processing (|Nielsen and Chuand . [2000h . 
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1. Spin-1/2 models with short range interactions 

A model Hamiltonian for a set of localized spins in- 
teracting with nearest neighbor exchange coupling on a 
rf-dimcnsional lattice is 



n(j,A,h z /jy- 



(i,3) 



E^ 2 



(23) 



In the previous expression i,j are lattice points, (•) con- 
straints the sum over nearest neighbors and S" (a — 
x,y, z) are spin-1/2 operators. A positive (negative) ex- 
change coupling J favors antiferromagnetic (ferromag- 
netic) ordering in the xy— plane. The parameters 7 and 
A account for the anisotropy in the exchange coupling in 
z direction, and h z is the transverse magnetic field. There 
are only very few exact results concerning ^(7, A, h z / J) 
in dimension d > 1. The ground state of Eq. (|2"3"|) is 
in general entangled. It exists however, for any value 
of the coupling constants 7 and A, J > a point in 
d= 1,2 (f or bipartite lattices) where the ground state i s 
factorized (jKurmann et al l ll982tlRoscilde et aZ.I . [2005bh . 
It occurs at the so called fa ctorizing field hf given by 
hf = I J-y/ (1 + A) 2 — (7/2) 2 where z is the coordination 
number. 

In d = 1 the model is exactly solvable in several im- 
portant cases. In the next two paragraphs we discuss the 
anisotropic quantum XY model ( A = and < 7 < 1) 
and the XXZ model (7 = 0). Also the XYZ- model in 
zero field, 7 7^ 0, A 7^ c an be solved exac tly but it will 
not be discussed here (see iTakatiashi. [19991 for a review). 



a. A = 0: Quantum XY -model The quantum Ising 
model corresponds to 7 = 1 while the (isotropic) XX- 
model is obtained for 7 = 0. In the isotropic case 
the model possesses an additional symmetry resulting 
in the conservation of the magnetization along the z- 
axis. For any v alue of the anisotropy the model can 
be sq l yed exac t ly dBarouch and McCovlll97J : iLieb et all . 
119611: iPfeutvl . Il97dh . By first applying the Jordan- 
Wigner transformation c k — e ITr °j 1 "j <j~ (with 
g ± = (l/2)(a x ±ia y )) the XY model can be transformed 
onto a free fermion Hamiltonian 



z\A. 



1,3^3 



Cj + 



~(cK iC t+h.c.) 



2 ^ 



A, 



(24) 

In the previous equations Ci,c\ are the annihilation 
and creation operators for the spinless Jordan- Wigner 
fcrmions. The two matrices A, B are defined as Aj^ — 
-J(h,j+x + Sj,k+i) - h z 5j^. and B^ k = -jj(5 k j+i - 
fij,k+i)- For the case of periodic boundary conditions 
on the spins, an extra boundary term appears in the 
fermionic Hamiltonian which depends on the parity of 
the total number of fermions Np- Notice that although 



Np does not commute with the Hamiltonian the parity of 
Np is conserved. A generic quadratic form, like Eq. (JMJ, 
can be diagonalized in terms of the normal-mode spinless 
Fermi operators by first going to the Fourier space and 
then performing a Bogoliubov transformation. 

The properties of the Hamiltonian are governed by the 
dimensionless coupling constant A = J/2h. In the in- 
terval < 7 < 1 the system undergoes a second order 
quantum phase transition at the critical value A c = 1. 
The order parameter is the magnetization in x-direction, 
(S x ), which is different from zero for A > 1. In the or- 
dered phase the ground state has a two-fold degereracy 
reflecting a global phase flip symmetry of the system. 
The magnetization along the z-direction, (S z ), is differ- 
ent from zero for any value of A, but presents a singular 
behavior of its first derivative at the transition. In the 
whole interval < 7 < 1 the transition belongs to the 
Ising universality class. For 7 = the quantum phase 
transition is of the Berezinskii-Kosterlitz-Thouless type. 

As it was discussed in Section Hi. Al and III.BI one- and 
two site- entanglement measures can be related to various 
equal-time spin correlation functions (in some important 
cases also the block entropy can be reduced to the evalu- 
ation of two-point correlators) M"(i) = (ip\Sf(t)\ip) and 
g ^(t) = (ip\Sf{t)S^(t)\ip). These correlators have been 
calculated for th is class of models in the case of ther- 
mal equilibrium dBarouch and McCovl . ll9Tlb ILieb et all , 
119611 : IPfeutvl Il970h . These can be recast in the form 
of Pfaffians that for stationary states reduce to Toeplitz 
determinants (i.e. determinants, whose entries depend 
only on the difference of their row and column num- 
ber), it can be demonstrated that the equal time cor- 
relat ion functions can be expre ssed as a sum of Pfaffi- 
ans (|Amico and OsterlohL |2004 . 



b. 7 = 0: XXZ model The two isotropic points 
A = 1 and A = — 1 describe the antiferromagnetic 
and ferromagnetic chains respectively. In one dimen- 
sion the XXZ Heisenberg model can be solved ex- 
actly by the Bethe Ansa tz technique ( Bethel . Il93lh 
(see e.g. iTakahashil . Il999l) and the correlation func- 
tions icaii_be_e2cp£essed in ter ms of certain determinants 
fsee iBogoliubov et all [l993 for a review). Correlation 
functions, especially for intermediate distances, are in 
general difficult to evaluate, alt hough important steps in 
this direction have been m ade ( Gohmann and Korepinl . 
2000: Kitanine et al. L 1 1999ft . 

The zero temperature phase diagram of the XXZ model 
in zero magnetic field shows a gapless phase in the inter- 
val — 1 < A < 1. Outside this interval the excitations are 
gapped. The two phases are separated by a Berezinskii- 
Kosterlitz-Thouless phase transition at A = 1 while at 
A = — 1 the transition is of the first order. In the pres- 
ence of the external magnetic field a finite energy gap 
appears in the spectrum. The universality class of the 
transition is not affected, as a result of the cons ervation 
of the total spin-z component (Takahaslul. Il999h . 
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When one moves away from one dimension, exact re- 
sults are rare. Nevertheless is it now established that the 
ground state of a two-dime nsional antiferromagnet pos- 
sesse s Neel long range order ( Dagottol . Il994t iManousakisl 
Il99lh . 



model (lAmico et all l200lt |Pi Lorenzo et all l2002t 
lOrtiz et alL\200$ ). 



3. Frustrated spin-1/2 models 



2. Spin-1/2 models with infinite range interaction 

In this case each spin interacts with all the other 
spins in the system with the same coupling strength 
n= [SfSf + ^SV] - E,h 4 • S, . For site- 

independent magnetic field hf — h a \/i, this model 
was orig inally proposed by Lipkin, Meshkoy and Glick 
(LMG) (iLipkin et all Il965t iMeshkov et all . Il965allbh to 
describe a collective motion in nuclei. In this case the 
dynamics of the system can be described in terms of a 
collective spin S a = £j ■ The previous Hamiltonian 
reduces to 



H = -^[(S x ) 2 + 1 (S y ) 2 ]-h-S 



(25) 



Since the Hamiltonian commutes with the Casimir op- 
erator S 2 the eigenstates can be labeled by the rep- 
resentation S of the collective spin algebra, at most 
linear in the number N of spins; this reduces (from 
2 N to N/2) the complexity of the problem. A further 
simplification is achieved at the supersymmetric point 
corresponding to J 2 j — 4h z , where the Hamiltonian 
can be factorized in two terms linear in the collective 
spin ( Unanvan and Fleischhauerll2003h ; then the ground 
state can be obtained explicitly. For a ferromagnetic cou- 
pling ( J > 0) and h x = h y = the system undergoes a 
second order quantum phase transiti on at A r = 1, charac - 
terized by mean field critical indices (jBottet et aLl . ll982h . 
The average magnetization (for any 7) m z — (S z )/N sat- 
urates for A < A c while it is suppressed for A > A c . For 



= the mode l exhib its a first 
(|Vidal et all l2006t ) while for 



h y = 0, h z < 1 and 7 
order transition at h x = 
an antiferromagnetic coupling and h y — a first order 
phase transition at h z = occurs, where the magnetiza- 
tion saturates abruptly at the same value m z — 1/2 for 
any 7's. 

The model Hamiltonian introduced at the beginning 
of this section embraces an important class of interact- 
ing fermion systems as well. By interpreting the non 
homogeneous magnetic field as a set of energy levels 
(h z )i = — £j, for h x = hy = and 7 = 1, it expresses 
the BCS model. This can be realized by noticing that 
the operators Sj := CjfCji, Sj~ := (Sjy 5| := 

( c ]| c jT + c jx c il — -0/2 span the su(2) algebra in the repre- 
sentation 1/2. In the fermion language the Hamiltonian 



t J 



reads U B cs = Ej,«r={t,i} € A° c i° ~ i 

Both the LMG and the BCS type models can 
be s olve d exactly by Bethe A nsatz dRichardsonl 



_ Richardson and Sherman! . 11964 ) (see 
also iDukelskv et all ~ 20041 for a, review) as they 
classical descendants of the six vertex 



Frustration arises in systems where certain local con- 
straints prevent the system from reaching a local energy 
minimum. The constraints can be of geometric nature 
(for example the topology of the underlying lattice) or of 
dynamical nature (two terms in the Hamiltonian tending 
to favor incompatible configurations). A classical exam- 
ple of the first type is that of an antiferromagnet in a 
triangular lattice with Ising interaction. At a quantum 
mechanical level this phenomenon can result in the ap- 
pearance of ground state degeneracies. The equilibrium 
and dynamical properties of frustrate d systems h ave been 
extensively studied in the literature (Dicp, 20Q5j) both in 
classical and quantum systems. 

A prototype of frustrated models in one dimension is 
the antiferromagnetic Heisenberg model with nearest and 
next-nearest neighbor interactions. This class of models 
were d iscussed originally to st udy the spin-Peierls tran- 



sition ( Schollwock et al. , |2004 ) . The Hamiltonian reads 



N 



H a = J 2j(Sj • Si+i + cvSj • S 



i+2) 



(26) 



are quasi 



Analytical calcul ations dHaldand. 19821). corrob orated by 
numerical result (jOkamoto and Nomural . 1 1992h indicate 
that at a ~ 1/4 there is a quantum phase transition to a 
dimerized 2-fold degenerate ground state, where singlets 
are arranged on doubled lattice constant distances. Such 
a phase is characterized by a finite gap in the low lying 
excitation spectrum. 

The Majumdar-Ghosh m odel (jMaiumdarl Il97fi 
iMaiumdar and Ghoshlll969al lbh is obtained from Eq.pfi]) 
for a — 1/2. The exact ground state can be solved by 
means of matrix product states (see next Section) and it 
is shown to be disordered. It is a doubly degenerate va- 
lence bond state made of nearest neighbor spin singlets. 
Although all two-point correlation function vanish, a fi- 
nite four-spin correlation function does reflect an ordered 
dimerization. 



4. Spin-1 models 

Spin-1 systems where originally considered to study 
the quantum dynamics of magnetic solitons in antif erro- 
magnets with single ion anisotropy (Mikeska, 1995). In 
one dimension, half-integ er and integer spi n chains have 
very different properties ( Haldane! . ll983al lbl). Long range 
order which is establis hed in the ground state of systems 
with half-integer spin ( Lieb et a/Ill96i[ ), may be washed 
out for integer spins. In this latter case, the system has 
a gap in the excitation spectrum. A paradigm model of 
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FIG. 1 A cartoon of the nearest neighbour valence bond state, 
exact ground state of the spin-1 model in Eq. (|27[) for /3 = 1/3 
(AKLT- model) . The ground state is constructed regarding 
every S = 1 in the lattice sites as made of a pair of S = 1/2, 
and projecting out the singlet state. The singlets are then 
formed taking pairs of S — 1/2 in nearest neighbor sites. 



particular state considered, if the state is only slightly 
entangled then the dimension o f A is bounded by some 
Dmps- MPS, first discussed in ( Fannes et aZ.l . ll992l ). ap- 
pear naturally in the Density Matrix Renorm alization 
Group (DMRG) (jOstlund and Rommerl . [l995h . In one- 
dimensional non-critical systems they describe faithfully 
the ground state. Infact, as shown by Vidal, matrix prod- 
uct states constit ute an efficien t representation of slightly 
entangled states (|Vidall . |200§) . 



interacting spin-1 systems is 



Si+i) 2 



(27) 



i=0 



The resulting gapped phase arises because of the pres- 
ence of zero as an eigenvalue of Sf; the corresponding 
eigenstates represent a spin excitation that can move 
freely in the chain , ultimately disordering the groun d 
state of the system (jGomez -Santosl Il99ll iMikeskal Il995l) . 
The so called string order parameter was proposed to 
capture the resulting 'floating' Neel order, made of al- 
ternati ng spins | 1), \ j) with strin gs of |0)'s in be- 
tween ( den Niis and Rommelsd . Il989l ) 



i+R-l 

Ostrtng = D lim (S"( TT ' 
k=i+l 



+RI 



(28) 



The ground state of physical systems described by Hamil- 
tonians of the form of Eq . (l27l) h as been studied in great 
details (jSchollwock et all |2004[) . Various phase transi- 



tions have been found between antiferromagnetic phases, 
Haldane phases, and a phase characterized by a large 
density of vanishing weights (Sf = 0) along the chain. 



a. The Affleck-Kennedy-Lieb-Tasaki (AKLT) Model. Some 
features of the phenomenology leading to the destruc- 
tion of the antiferromagnetic order can be put on a 
firm ground for (3 = 1/3 (AKLT model), where the 
groun d state of the Hamilt onian in Eq. (f2"Tj) is known ex- 
actly ( Affleck et all Il988l ). In this case it was proved 



that the ground state is constituted by a sea of nearest 
neighbour valence bond states, separated from the first 
excitation by a finite gap with exponentially decaying 
correlation functions. Such a state is sketched in FigfT] 
In fact it is a Matrix Product State (MPS), i.e. it belongs 
to the class of states which can be expressed in the form 



D 



i>MPS) = 



TrAf 



A SN 
■ N 



Si, ■ 



,sn) (29) 



si,...,sjv 



where the matrices (A s k z )i m parametrize the state; |sj > 
denotes a local basis of the D-dimensional Hilbert space; 
the trace contracts the indices l,m labelling bond states 
of the auxiliary system (namely the spin 1/2 for the 
AKLT model). The dimensions of A depends on the 



B. Strongly correlated fermionic models 

The prototype model of i nteracting fermion s on a lat- 
tice is the Hubbard model (Es sler et all. I2004TI 



« = -*Ei 

(ij) 



C i,cr C j,<? 



h.c.} + U^2ni^n ltl - /j,N (30) 



where Cj )(T , c\ a are fermionic operators: {c^, = 
Sij5 aa '. The coupling constant U describes the on-site 
repulsion, t is the hopping amplitude. 

The Hubbard model possess an u(l) © su(2) sym- 
metry expressing the charge conservation: u(X) = 
spanjA^ = ^2j a fija-} and the invariance under spin 
rotation: su(2) = spanlS 12 = X)j( n iT — n jl)^ + = 

J2j c l x c j,l' S~ = (S + Y}- Such a symmetry allows one to 
employ the total charge and magnetization as good quan- 
tum numbers. At half filling n = N/L = 1 (jj, = U/2) 
the symmetry is enlarged to so(4) = su(2) © su(2) by the 
generator rj = X^i(~) J ' c .?',T c .7M together with its hermitean 
conjugate (|Yang and Zhand . ll990l ). It was demonstrated 
that |W) = (ri) N \gs) are eigenstates of the Hubbard 
model (in any dimension), characterized by off-diagonal- 
long-ran ge-order via the mechanism of the so called r\- 
pairing (jYand . 1 19891 ). 

In one dimension the Hubbard model undergoes a 
Mott transition at U — of the Berezinskii-Kosterlitz- 
Thoul ess type. B y mea ns of the Bethe Ansatz solu- 
tion ( Lieb and Wul . Il968f) it can be demonstrated how 
the bare electrons decay in charge and spin excitations. 
The phenomenon of spin-charge separation occurs at low 
energies away from half filling. For repulsive interaction 
the half- filled band is gapped in the charge sector; while 
the spin excitations remain gaplcss. The mirror-inverted 
situation occurs for attractive interac tion where the gap 
is in the spin excitations instead fsee lEssler et al. 1. 120041 
for a recent review). 

The Hubbard model in a magnetic field was proved to 
exhibit two quantum critical points at hff 1 = 4(|Z7| ± 1) 
and half filling for U < 0, while there is one at h r c — 
4(Vt 2 + U 2 - U) for U > (lYang et all l2000h . 

If a nearest neighbor Coulomb repulsion 



,t 



TvjaTvj+la' 



is taken into account in Eq. ([3H)l . 



a spin density wave and a charge density wave phase 
appear. A transition to a phase separat ion of high 
densi ty and low density regions (see e.g. Ref. IClav et all . 
1999) is also present. 
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The bond charge extended Hubbard model, origi- 
nally proposed in t he context of high T c superconduc- 
tivity (jHirschl . H989V include further correlations in the 
hopping process already involved in (|30|) . The Hamilto- 
nian reads 



H 



^ni yl -t[l-x(ni- a +n i+1 .- a )]c\ a Ci +lya +h.c. 

(31) 

(for x — the Eq. (|3Tj) coincides with the Hubbard model 
([3H)l ). For x 7^ the hopping amplitudes are modu- 
lated by the occupancy of the sites involved in the pro- 
cesses of tunneling. Because of the particle-hole symme- 
try, x can be restricted in [0, 1] without loss of general- 
ity. For x — 1 the correlated hopping term commutes 
with the interaction. In this case the exact ground state 
was shown to exhibit a variety of quantum phase tran- 
sitions between insulators and superconducting regimes, 
controlled by the Coulomb repulsion parameter U. For 
x = 1 the phase diagram is shown in Section HVl FigJTS] 
At U/t = 4 and n = 1, a superconductor- insulator 
quantum phase transition occurs; for —4 < U/t < 4 
the ground state is characterized by off-diagonal long- 
range order; the low lying excitations are gapless. For 
U/t = —4 a further quantum critical point projects the 
ground state into the Hilbert subspace spanned by singly 



ground state into tne MUDert subspace spanned by singly 
and doubly occupied states ( Arrachea and Aligial . Il994t 
ISchadschneideii Il995l ) . For intermediate x the model has 
not been solved exactly. Numerical calculations indicate 
a superconducting-insulator transition controlled by U 
and parametrized by x. Specifically, for < x < 1/2 the 
phase is gapped at any non vanishing U ; for 1/2 < x < 1 
the onset to a superconduct i ng pha se was evidenced at 
some finite U ( Anfossi et aUl2005al) . 



the spin-bosons interaction and Ao is a ultraviolet cut- 
off frequency. The power s characterizes the bath. For 
s = 1 the bath is called Ohmic, in this case the model 
has a second order quantum phase transition at a = 1 
from under-damped to over-damped oscillations (where 
the value of the spin is frozen). The value a = 1/2 
identifies a cross over regime where the two-level sys- 
tem is driven from coherent to incoherent oscillations. 
If the bath is super-Ohmic (s > 1), the quantum crit- 
ical point is washed out, while a cross over occurs at 
a ~ log(A /<5). For sub-Ohmic baths (s < 1), sev- 
eral studies indicate the existe n ce of a quantum crit- 
ical point (jSpohn and Dumckel . Il985l). The question 



ical point ipponn and uumcm iiijaol). ine q u estion 
however is not co mplete l y settled dBulla et all 20031 : 
iKehrein and Mielkd . Il996t IStauber and Mielkd . l2002fh 

An interesting case is also that of a spin interacting 
with a single bosonic mode, A n = A<5 n) o- 

H = - S -S x + wo(aJao + \) + y S z (aJ + a ) (33) 

Such model describes for example an atom in- 
teracti ng with a monochr o matic electromagnetic 
field (ICohen-Tannoudii et all 19921) via a dipole 
force (jJavnes and Cummingd . Il963l ). Recently, the 
dynamics corresponding to (1331) wa s inten sely studied in 
relation to io n traps (ICirac et all Il992t) and quantum 
computation ( Hughes et all Il998f) . The model defined 
in Eq. (|33|) with 5 = 1/2 (Jaynes-Cummings model) 
was generalized and solved exac tly to consider generic 
spin (|Tavis and Cummingsl . Il969t ) in order to discuss the 
super-radiance phenomenon in cavity-QED. 



D. Harmonic lattices 



C. Spin-boson models 

A prototype model in this class is that of a quan- 
tum system coup led to a bath of harmonic oscillators 
(see IWeisd . 1999 for a review of open quantum mechan- 
ical systems) known also as the Caldeira-Leggett model. 
In this case the quantum system is a two level system. 
This class of models was intensely investigated to study 
the quantum-to-classical transi tion and the correspond- 
ing loss of quantum coherence ( Zurekl . l2003f h 

The spin-boson Hamiltonian has the form 

H sb = -^S x + ^„(4a„ + i) + is 2 ^A„(4 + a„) , 

n n 

(32) 

it can be demo nstrated to be equivalent to the aniso tropic 
Kondo model (jAnderson et all 11970 : lGuineal . ll985l ). The 
coupling constants {A„} fix the spectral density of the 
bath: J(cj) = (7r/2) A^<5(w — uj n )/co n . At low energy 

n 

the spectral function can be represented as a power-law: 
J(uj) cx 2au} s A^ s where a is the parameter controlling 



The Hamiltonian for a lattice of coupled har- 
monic oscillators (short: harmonic lattice) can be ex- 
pressed in terms of the phase space vector £ T = 
(<7i, . . .,q n ;Pi, ■ ■ ■ ,p n ) as 



h = r 











(34) 



where U is the n x n interaction matrix for the co- 
ordinates. If the system is translational invariant the 
matrix U is a Toeplitz matrix wi th periodic bound- 
aries , also called circulant matrix (jHorn and Johnson . 
1994). In the case of finite range interaction of the form 
J2 r Sfe=i Krilk+r ~~ Ik) 2 and assuming periodic bound- 



ary conditions, its entries are U 



l + 2£r 



and 



Ujj+ r = — a r with a r — 2K r /mu) 



2 . Since the Hamilto- 
is quadratic in the canonical variables its dy- 
namical algebra is sp(2n,WL). Then the diagonalization 
can be achieved by RHR^ 1 where R = ®£=i exp (i6 a G a ) 
where G a is the generic Hermitean element of sp(2,IR). 

As we discussed in Section III.GI the key quantity 
that characterizes the properties of harmonic sys- 
tems is the covariance matrix defined in Eq. (|20[) . 
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For the resulting decoupled harmonic oscillators it is 
diag{(m^/^"a;)" 1 , . . . , (m v /^w)~ 1 ; to v /77i"w, . . . , m^w), 
where r\j are the eigenvalues of U. Employing the virial 
theorem for harmonic oscillators, the covariance matrix 
for a thermal state with inverse temperature (3 = 1/ksT 
can be calculated as well 



V 



(mwVU)- 1 ^ 








(mwVU)N, 3 



(35) 



where = 1„ + 2/(exp (/3wVU) - 1„) The range of the 
position or momentum correlation functions is related to 
the low lying spectrum of the Hamiltonian. For a gapped 
systems the correlations decay exponentially. The ab- 
sence of a gap (some eigenvalues of U tend to zero for an 
infinite system) leads to critical behavior of the system 
and characteristic long ranged correlations. A rigorous 
and detailed discussion of the relations between the gap 
in the energy spectrum and the properties of the corre- 
lations can be found in (|Cramer and Eisertl . 120061 ) . 



IV. PAIRWISE ENTANGLEMENT 

At T — many-body systems are most often described 
by complex ground state wave functions which contain 
all the correlations that give rise to the various phases 
of matter (superconductivity, ferromagnetism, quantum 
hall systems, . . .). Traditionally many-body systems have 
been studied by looking for example at their response 
to external perturbations, various order parameters and 
excitation spectrum. The study of the ground state of 
many-body systems with methods developed in quantum 
information may unveil new properties. In this Section 
we classify the properties of the ground state of a many- 
body system according to its entanglement. We concen- 
trate on spin systems. Spin variables constitute a good 
example of distinguishable objects, for which the prob- 
lem of entanglement quantification is most developed. 
We will discuss various aspects starting from the pair- 
wise entanglement, we then proceed with the properties 
of block entropy and localizable entanglement. Most of 
the calculations are for one-dimensional systems where 
exact results are available. Section IIV.BI will overview 
the status in the d-dimensional case. Multipartite en- 
tanglement in the ground state will be discussed later in 
Section LYIIIl 



A. Pairwise entanglement in spin chains 

1. Concurrence and magnetic order 

The study of entanglement in interacting spin sys- 
tems was i nitiated with the works on isotropic Heisen- 
berg rings (lArnesen et all l200lb iGunlvcke etaH . \200ll 
lO'Connors and Woottersl . l200ll) . O'Connors and Woot- 
ters aimed at finding the maximum pairwise entangle- 
ment that can be realized in a chain of N qubits with peri- 



odic boundary conditions. Starting from the assumption 
that the state maximizing the nearest neighbor concur- 
rence C(l) were an eigenstate of the z component of the 
total spin (llshizaka and Hiroshimal l2000l : IMunro et all 
120011 : IVerstraete et all l200ll ) the problem was recast to 
an optimization pr ocedure simi lar in spirit to the coordi- 
nate Bethe ansatz ( Bethel . fl93lf) : 4 : the search for the op- 
timal state was restricted to those cases which excluded 
the possibility to find two nearest neighbor up spins. For 
fixed number of spins N and p spins up, the state can be 



written as 



= Ei 



<ii< - <i p <N 



,) (b are 



the coefficients and ij are the positions of the up-spins) 
therefore mapping the spin state onto a particle state 
such that the positions of the p particles correspond to 
those of the up spins. The maximum concurrence within 
this class of states could be related to the ground state 
of this gas of free spinless particles with the result 



C(l) = ~E gs 



N—p 

sin -jrf~ — 

N — p 



(36) 



Eq. (|36p gives a lower bound for the maximal attain- 
able concurrence. The isotropic antiferromagnetic chain 
was considered as the physical system closest to a per- 
fectly dimerized system (classically, with alternating up 
and down spins). It was noticed however that the con- 
currence of the ground state of the antiferromagnetic 
chain is actually smaller than the value of the ferromag- 
netic cha in, indicating that the situation is more com- 
plicated (jO'Connors and Woottersl l200lh . In order to 
clarify this point, a couple of simple examples are use- 
ful. For a system of N = 2 spins the ground state 
is a singlet. However for general N (we assume an 
even number of sites) the ground state is not made of 
nearest-neighbor singlets (Resonant Valence Bond (RVB) 
state). For example the N = 4 the ground state is 
\gs) = (l/v / 6)[2|0100)+2|1000)-|1001)-|0110)-|0011)- 
1 1100)] , different from the product of two singlets. It can 
be seen that the effect of the last two components of 
the state is to reduce the concurrence with respect to 
its maximum attainable value. Given the simple rela- 
tion Eq. (|36|) between the nearest neighbor concurrence 
and the ground state energy, the deviation from the 
RVB state can be quantified by looking at the difference 
from the exact ground state energy corresponding to the 
maximum concurrence. This maximum value is reached 
within the set of eigenstates with zero total magnetiza- 
tion (the 'balanced' states in lO'Connors and Woottersl 
2001), indicating that the concurrence is maximized only 
on the restricted Hilbert space of z-rotationally invariant 
states. Indications on how to optimize the concurrence 
were discussed in ( Hiesmavr et all 120061 : iMever et all 



4 Such a method relies in the existence of a 'non interacting picture' 
where the wave function of the system can be written as a finite 
sum of plane waves; the ansatz is successful for a very special 
form of the scattering among such non interacting pictures. 
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|2004| ). The solu tion to the problem for N — > oo was 
recently given in (jPoulsen et al. , 2006). It turns out that 
the state s with nearest-neighbors aligned spins (not in- 
cluded in lO'Connors and Woottersl 12001). correspond to 
a 'density-density' interaction in the gas of the spinlcss 
particles considered above, that hence are important for 
the analysis, (in the analogy of the coordinate Bethe 
ansatz method, this p rovides the 'inter acting picture'). 
Following the ideas of (jWolf et all 120031 ) , the problem to 
find the optimum concurrence was shown to be equivalent 
to that of finding the ground state energy of an effective 
spin Hamiltonian, namely the XXZ model in an exter- 
nal magnetic field. The optimal concurrence is found in 
the gapless regime of the spin model with a magnetiza- 
tion M z = 1 — 2p/ N. It was further demonst r ated that 
states considered in (lO'Connors and Woottersl l200lh ac- 
tually maximize the concurrence for M z > 1/3 (for 
< M z < 1/3 the states contain nearest neighbor up- 
spins). 

The concurrence, beyond nearest-neighbors, in 
isotropic Heisenberg antiferroma gnets in an external 
magnetic fie l d wa s discu s sed in ( Arnesen et all l200lt 
iFubini et all 120061 : IWand . l2002al) . The combined effect 
of the magnetic field a nd the anisotropy in He isenberg 
magnets was studied in ( Jin and KorepinLl2004ah making 



use of the exact results existing for the one-dimensional 
XXZ model. It turns o ut that the concurrence increase s 
with the anisotropy A ( Kartsev and Karshnukovi [2004). 
For strong magnetic fields the entanglement vanishes (the 
order is ferromagnetic) ; for large values of the anisotropy 
A the state is a classical Neel state with Ising order. Ex- 
cept for these cases, quantum fluctuations in the ground 
state lead to entangled ground states. 

As we discussed in Sec. IIII1 in low dimensional spin 
system there exists a particular choice of the cou- 
pling constairts_Jor wliich the ground state is factor- 
ized (jKurmann et a/.I .IT982). This is a special point also 
from the perspective of investigating the entanglement in 
the ground state. Several works were devoted to the char- 
acterization of the entanglement close to the factorizing 
point. It turns out that the point at which the state of 
the system becomes separable marks an exchange of par- 
allel and anti-parallel sector in the ground state concur- 
rence, see Eq.([7]). As this phenomenon involves a global 
(long-range) reorganization of the state of the system, 
the range of the concurrence diverges. (We notice that 
several definitions of characteristic lengths associated to 
entanglement decay exist). The concurrence is often ob- 
served to vanish when the two sites are farther than R 
sites apart: the distance R is then taken as the range of 
the concurrence. For the XY model it was found that 
this range is 



R cx In 



1-7 



In I A" 



A7Y 



(37) 



The divergence of R suggests, as a consequence of the 
monogamy of the entangleme nt ( Coffman et all l2000t 
lOsborne and Verstraetd . 120061 ). that the role of pairwise 



entanglemen t is enhanced while approa ching the sepa- 
rable point (|Roscilde et all 12004 l2005allbT ). In fact, for 
the Ising model (i.e. 7 = 1), one finds that the ratio 
t^/ti — > 1 when the magnetic field a pproaches the fac- 
torizing field hf (jAmico et M$, 120061) . For 7 / 1 and 
h[ < h z < h c it was found that t^/ti monotonically 
increases for h z — > hf and that the value {T%/T\)\ h + in- 
creases with 7 — ► 1. The existence of the factorizing has 
been also point ed out in other one dimensional systems 
both f or short (jAmico et 

Ml 120061: iRoscilde et ali, [2 004. 
2005a]) and long range interactions (|Dusuel and Vidal . 
l2005h . In all these cases the range of the two-site entan- 
glement diverges. The range of the concurrence was also 
studied for the XXZ (|Jin and Korepinl l2004ah where it 
was shown to vary as 



R 



2A(h z ) 
1 - 4M? 1 



(38) 



The exponent 9 = 2 for finite fields, while it is 6 = 1 for 
h = 0; the coeffi cient A(h) i s known exactly in the para- 
magn eti c phase (iLukvanovl Il999t Lukyanov and Terrasl 



2003]; iLukv anov and Zamolodchikovl Il997h (van- 
ishing magnetization) and in the saturation 
limit ( Vaidva and Tracvl . Il979alibh For generic h it 
was calculated numerically in ( Hikihar a and Furusakil . 
|2004| ). For the isotro pic Heisenberg antiferromagnet, 
R = l (|Gu et aLl . [2003h . 



In all the previous cases the increasing in the range 
of the pairwise entanglement means that all the pairs 
at distance smaller than R share a finite amount 
of entanglement (as quantified by the concurrence). 
There are one-dimensional spin systems where the pair- 
wise entanglement has qualitative different features as 
a function of the distance between the sites. An 
exa mple is the long-distance entanglement observed 
in ( Campos Venuti et all l2006bl ). Given a measure of 
entanglement E(pij), Campos Venuti et al showed that 



it is possible that E(pij) ^= when 



■3 



the 



ground state. Long-distance entanglement can be real- 
ized in various one-dimensional models as in the dimer- 
ized frustrated Heisenberg models or in the AKLT model. 
For these two models the entanglement is highly non uni- 
form and it is mainly concentrated in the end-to-end pair 
of the chain. 

Spontaneous symmetry breaking can influence the en- 
tanglement in the ground state. To see this, it is con- 
venient to introduce the "thermal ground state" po = 
3 (\9s°)(gs°\ + \gs e )(gs e \) = \ (\gs-)(gs-\ + \gs+)(gs+\) 
which is the T — > limit of the thermal state. In the 
previous expression gs + and gs~ are the symmetry bro- 
ken states which give the correct order parameter of the 
model. They are superpositions of the degenerate par- 
ity eigenstates gs° and gs e . Being con vex, the concur- 



rence m gs^ will be larger than in gs°^ e ( Osterloh et 



l2006h . The opposite is true for the concave entrop y of 



entanglement (see Ref. ( Osborne and Nielsen. l2002t) for 
the single spin von Neumann entropy) . The spontaneous 
parity symmetry breaking does not affect the concur- 
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rence in the ground state as long as it coincides with 
C 1 , Eq.®: that is, i f the spins are ent angled in an an- 
tiferromagnetic way ijSvliuasenl . l2003bh . For the quan- 
tum Ising model, the concurrence coincides with C 1 for 
all values of the magnetic field, and therefore, the con- 
currence is unaffected by the symmetry breaking, the 
hallmark of the present QPT. For generic anisotropics 
7 instead, also the parallel entanglement C 11 is observed 
precisely for magnetic fields smaller than the factoriz- 
ing field ( Osterloh et aUl2004| ): this interval excludes the 
critical point. This changes at 7 = 0, where the concur- 
rence indeed shows an infinite range. Below the critical 
field, the con currence is enhanced by the parity symme- 
try breaking (|Osterloh et all . l200d ) 



2. Pairwise entanglement and quantum phase transitions 

A great number of papers have been devoted to the 
study entanglement close to quantum phase transition 
(QPTs). QPT occur at zero temperature. They are in- 
duced by the change of an e xternal parameter or cou- 
pling constant ( Sachdevlll999f) . Examples are transitions 
occurring in quantum Hall systems, localization, the 
superconductor-insulator transition in two-dimensional 
systems. Close to the quantum critical point the sys- 
tem is characterized by a diverging correlation lenght 
£ which is responsible for the singular behavior of dif- 
ferent physical observables. The behavior of correlation 
functions however is not necessarily related to the behav- 
ior of quantum correlations present in the system. This 
question seem particularly interesting as quantum phase 
transitions are associated with drastic modifications of 
the ground state. 

The critical properties in the entanglement we are go- 
ing to summarize below allow for a screening of the qual- 
itative change of the state of the system experiencing 
a quantum phase transition. In order to avoid possible 
confusion, it is worth to stress that the study of entan- 
glement close to quantum critical points does not pro- 
vide new understanding to the scaling theory of quantum 
phase transitions. Rather it may be useful in a deeper 
characterization of the ground state wave function of the 
many-body system undergoing a phase transition. In this 
respect it is important to understand, for instance, how 
the entanglement depends on the order of the transition, 
or what is the role of the range of the interaction to es- 
tablish the entanglement in the ground state. In this 
section we discuss exclusively the pairwise entanglement 
while in the next section we approach the same problem 
by looking at the block entropy 5 . 

Pairwise entanglement close to quantum phase tran- 
sitions was originally analyzed in ( Osborne and Nielsenl . 



l2002t lOsterloh et all . 12002ft for the Ising model in one 
dimension. Below we summarize their results in this spe- 
cific case. The concurrence tends to zero for A ^ 1 and 
A -C 1, the ground state of the system is fully polar- 
ized along the cc-axes (z-axes). Moreover the concur- 
rence is zero unless the two sites are at most next-nearest 
neighbors, we therefore discuss only the nearest neigh- 
bor concurrence C(l) (see however Section IIV.A.1I for 
cases where there is a longer-range pairwise entangle- 
ment). The concurrence itself is a smooth function of the 
coupling with a maximum close to the critical point (see 
the right inset of Fig(2J) ; it was argued that the maximum 
in the pairwise entanglement does not occur at the quan- 
tum critical point because of the monogamy property (it 
is the global entanglement that should be maximal at the 
critical point). The critical properties of the ground state 
are captured by the derivatives of the concurrence as a 
function of A. The results for systems of different size 
(including the thermodynamic limit) are shown in Figj2] 
For the infinite chain d\C(l) diverges on approaching the 
critical value as 



dxC(l) 



3tt 2 



lnlA- A c 



(39) 



For finite system the precursors of the critical behavior 
can be analyzed by means of finite size scaling. In the 
critical region the concurrence depends only on the com- 
bination N X / U (X — X rn ) where v is the critical exponent 
governing the divergence of the correlation lenght and 
A m is the position of the minimum (see the left inset of 
FigEJ. In the case of log divergence the scaling ansatz 
has to be adapted and takes the form d\C(l)(N, X) — 
d x C(l){N,X ) ~ QlNV'SnX] - QlN^SmXo] where A 
is some non critical value, S m (X) = X — X m and Q{x) ~ 
Q(oo)lna; (for large x). Similar results have been ob- 



5 QPTs were also s t udied by looking at qua ntum fi- 
delity llCozzini et al I 120061: IZanardi et all )2006t) or the 
effect of single bit operations (Giampaolo et al., 2006a. b) 



taine d for the XY universality class (jOsterloh et all . 
l2002f ). Although the concurrence describes short-range 
properties, nevertheless scaling behavior typical of con- 
tinuous phase transition emerges. 

For this class of models the concurrence coincides with 
C 1 in Eq.® indicating that the spins can only be en- 
tangled in a antiparallel way (this is a peculiar case of 
7=1; for generic anisotropies also the parallel entangle- 
ment is observed). The analysis of the finite size scaling 
in the, so called, period-2 and period-3 chains where the 
exchange coupling varies every second and third lattice 
sites respecti vely, leads to the same scaling laws in the 
concurrence (Zha ng and Tonsi 120051 ). 

The concurrence was found to be discontinuous at the 
first o rder ferromagne tic t ransition A = — 1 in the XX Z 
chain ( Gu et a/ll2003l )(see lGlaser et all . 120031 for explicit 
formulas relating the concurrence and correlators for the 
XXZ model in various regimes). This result can be un- 
derstood in terms of the sudden change of the wave func- 
tion occurring because of the level crossing characteriz- 
ing these type of quantum critical points. The behavior 
of the two-site entanglement at the continuous quantum 
critical point of the Kosterlitz-Thouless type A = 1 sepa- 
rating the XY and the antiferromagnetic phases is more 
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formed ones). This behavior can be t raced back to th e 
properties of the finite size spectrum ( Gu et g?l l2005al) . 
In fact, at A = 1 the concurrence can be related to the 
eigenenergies. The maximum arises since both the trans- 
verse and longitudinal orders are power law decaying at 
this critical point, and therefore the excited states con- 
tribute to C(l) maximally. 

Studies of finite size energy spectru m of other mod- 
cls like the dimerized Heisenberg chain (jSun et "oE l2005h 
and Majumdar-Ghosh model (Eq.(|26|) with a = 1/2) 
show how level crossings in the energy spectrum affect 
the behavior of the bipartite e ntanglement occu rring at 
the quantum phase transition ( Gu et all , l2005at) . 



FIG. 2 The derivative of the nearest neighbor concurrence as 
a function of the reduced coupling strength. The curves corre- 
spond to different lattice sizes. On increasing the system size, 
the minimum gets more pronounced and the position of the 
minimum tends as (see the left side inset) towards the critical 
point where for an infinite system a logarithmic divergence 
is present. The right hand side inset shows the behavior of 
the concurrence for the infinite system. The maximum is not 
related to the critical properties of the Ising model. [From 
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FIG. 3 Nearest neigh bor concurrence of the XX Z model. 
[From (|Gu et all Wm ] 



complex. In this case the nearest neighbor concurrence 
(that is the only non vanishing one) reaches a maximum 
as shown in Fig(3J Further understanding on such be- 
havior can be achieved by analyzing the symmetries of 
the model. At the antiferromagnetic point the ground 
state is an su(2) singlet where nearest neighbor spins 
tend to form singlets; away from A = 1, this behav- 
ior is 'deformed' and the system has the tendency to to 
reach a state of the type <8>|<^) made of 'g-deformed sin- 
glets' correspon ding to the quantum algebr a su q (2) with 



2 A = q + q- 1 (jPasauier and Saleurl . Il99dh . This allows 



to rephrase the existence of the maximum in the con- 
currence as the loss of entanglement associated to the 
q-deformed symmetry of the system away from A = 1 
(note that q-singlets are less entangled than the unde- 



a. LMG model Because of the symmetry of the LMG 
models (see Eq |25]) any two spins are entangled in the 
same way. The concurrence C is independent on the two 
site indices, it can be obtained exploiting the explicit 
expression of the eigenstates. Due to the monogamy 
of entanglement the result must be rescaled by the co- 
ordination number, Cr = (N — 1)C, to have a fi- 
nite value in the thermodynami c. For the ferromag- 
netic model (|Vidal et all |2004al) . it was proved that 
close to the continuous QPT, A = 1 characterizing 
the ferromagnetic LMG model, the derivative of the 
concurrence diverges, but, differently from Ising case, 
with a power law. It is interesting that Cr can 
be related to, the_so called, spin sque e zing p arame- 
tcr S = 2 v / AS n± (IWang and Sanders! . 12003 ). mea- 
suring the spin fluctuations in a quantum correlated 
state (he subscript n± indicates a perpendicular axes 
to (S)). The relation reads £ = \/\ — Cr According 




FIG. 4 The rescaled concurrence of the antiferromagnetic 
LMG model. The fi rst order transition occurs at h = 0. [From 
(jVidal et all l2004al )] 



to ( Lewenstein and Sanperal . [l998) the two-spins re- 
duced density operator can be decomposed in a sepa- 
rable part and a pure entangled state p e with a cer- 
tain weight A. Such a decomposition leads to the re- 
lation C(p) = (1 — A)C(p e ). Critical spin fluctuations 
are related to the concurrence of the pure state C(p e ) 
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FIG. 5 The rescaled concurrence between two atoms in the 
Dicke mode. The concurrence is rescaled in both for finite N 
and in the thermod ynamic limit. The in set shows the finite 
size scaling. [From ([Lambert et qZ.l . l2004 )] 



while the diverging correlation lenght is related to the 
weight A (Shimizu and Kawaguchii I2006T ) . The analy- 
sis of critical entanglement at the first order quantum 
critic al point of the antif erromagnetic LMG model shows 
that ( Vidal et all l2004ah the discontinuity is observed 
directly in the concurrence also for spin interacting with 
a long-range, see FigdJ 



b. Pairwise entanglement in spin-boson models We first 
discuss the Tavis-Cummings model defined in Eq. (|33[) . 
In this model the spin S is proportional to the number 
of atoms, all interacting with a single mode radiation 
field. The pairwise entanglement between two different 
atoms undergoing the super-radiant quantum phase tran - 
sition ([Lambert et all |2004 120051 : iReslen et all 120051 ) 
can be investigated through the rescaled concurrence 
CV = NC, see FigEl similarly to what has been dis- 
cussed above for the LMG models. In the thermo- 
dynamic limit the spin-boson model can be mapped 
onto a quadratic bosonic system trough an Hqlstein - 
Primakoff transformation ( Emarv and Brandesl [2003). 
Many of the properties of the Tavis-Cummings model 
bear similarities with the ferromagnetic LMG model. 
In the thermodynamic limit the concurrence reaches a 
maximum value 1 — V2/2 at the super radiant quan- 
tum phase_transn]ir^ root singularity (see 
also ISchneider and Milburnl . 120021 ). The relationship be- 
tween the squee zing of the state and e ntang lement was 
highlighted in ( Sprensen and Molmerl. 1200 if) and ana- 
lyzed in more details in (jStockton et all 120031 ) where it 
was also suggested how to deal with entanglement be- 
tween arbitrary splits of symmetric Hilbert spaces (like 
the Dicke states span). 

Entanglement between qubits and a single mode 
and between two spins with an Heisenberg interaction 
of the XXZ type, additionally coupled to a single 
bosonic field wa s cons idered in ( Liberti et all 12006 
and in ( He et a/ll2006f ) respectively. 



3. Entanglement versus correlations in spin systems. 

From all the results summarized above it is clear 
that the anomalies characterizing the quantum critical 
points are reflected in the two-site entanglement. At 
a qualitative level this arises because of the formal re- 
lation between the correlation functions and the entan- 
glement. A way to put this observation on a quantita- 
tive groun d is provid e d by a generalized Hohenberg-Kohn 
theorem ( Wu et all , 120061 ). Accordingly, the ground 
state energy can be considered as a unique function of 
the expectation values of certain observables. These, 
in turn, can be related to (the various derivatives of) 



a given entangl e ment measure ([Campos Venuti et all 
l2006aHWu eraLl . 120041) . 

Specifically, for an Hamiltonian of the form H = 
Hq + XiAi with control parameters A; associated with 
operators Ai, it can be shown that the ground state re- 
duced operators of the system are well behaved functions 
of (Ai). Then, any entanglement measure related to re- 
duced density operators, M = M (p) is a function of (Ai) 
(in absense of ground state degeneracy) by the Hellmann- 
Feynman theorem: OE/dXi — (dE/d\i) = (Ai). There- 
fore it can be proved that 



M((A l )) = M(— ) 



(40) 



where E is the ground state energy. From this relation 
it emerges how the critical behavior of the system is re- 
flected in the anomalies of the entanglement. In partic- 
ular, first order phase transition are associated to the 
anomalies of M while second order phase transitions cor- 
respond to a singular behavior of the derivatives of M. 
Other singularities like those noticed in the concur rence 
for models with three-spin interactions ( Yand . 120051 ). are 
due to the non-analyticity intrinsic in the definition of 
the concurrence as a maximum of two analytic functions 
and the constant zero. 

The relation given in Eq. ([4U)l was constructed ex- 
plicitl y for th e quan tum Ising, XXZ, and LMG mod- 
els (|Wu et all l2006th For the Ising model: J2i XA i = 
h J^i 3f I the divergence of the first derivative of the con- 
currence is then determined by the non analytical be- 
havior of (S x S x ) $Wu et all 12004 ). For the XXZ model: 
J2i hAi — AJ2i SfSf +1 . At the transition point A = 1 
both the purity and the concurrence display a maximum. 
It was proved that such a maximum is reflected also 
in a stationary point of the ground state energy as a 
function of {S*Sf +1 )] the concurrence is continuous since 
the Berezinskii-Kosterlitz-Thouless transition is of infi- 
nite order. A relevant caveat to Eq. (|40| is constituted 
by the uniaxial-LMG model in a transverse field (with 
h y = and 7 = 0) that displays a first order QPT for 
h x = 0. The concurrence is continuous at the transition 
since it does not depend on t he discontinuous ele ments 
of the reduced density matrix ( Vidal et all l2004cj) . 

The relation between entanglement and criticality 
was also studied in the spin-1 XXZ with single ion 
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anisotropy. It was established that the critical anomalies 
in the entropy experienced at the Haldane-large- D (if an 
axial anisotropy Dj2i(Si ) 2 is added to the Hamiltonian 
in Eq. (|2"T|) ) transition fans out f rom the singularity of 
the loc al order parameter ((S z ) 2 ) (jCampos Venuti et all 
l2006ah . 

A way to study the general relation between entan- 
glement and crit i cal p henomena was also pursued in 
(jHaselgrove et alL\2004i . It was argued how for systems 
with finite range interaction a vanishing energy gap in 
the thermodynamic limit is an essential condition for the 
ground state to have non-local quantum correlations be- 
tween distant subsystems. 



4. Spin models with defects 

The problem to characterize entanglement in chains 
with defects was addressed first for the quantum XY 
models with a single defe ct in the exchange in teraction 
term of the Hamiltonian ( Osenda et all 120031 ). It was 
found that the effect of the impurity is to pin the entan- 
glement. Moreover the defect can induce a pairwise en- 
tanglement on the homogeneous part of the system that 
were disentangled in the pure system. Even at the quan- 
tum critical point the finite size scaling of the critical 
anomaly of the concurrence is affected by the distance 
from the impurity. This basic phenomenology was ob- 
served in a variety of different situation that we review 
below. 

The presence of two defec ts has been an alyzed in the 
XXZ chain. It turns out (|SantosL 120031) that various 
type of entangled states can be created in the chain by 
spin flip excitations located at the defects positions. The 
entanglement oscillates between the defects with a period 
that depend on their distance. Also the anisotropy A of 
the chain is a relevant parameter controlling the entan- 



glement between the def ects. Small anisotropies ca n sup- 
press the entanglement (|Santos and Rigolinl . [2005) . The 
way in which this kind of localization ca n be exploited 
for q uantum algorithms was studied in ( Santos et all 
2005). The entanglement was also studied in systems 
with defects in the presence of an external magnetic 
field (|Apollaro and Plastinal . l2006f) . It was demonstrated 
that such defect can lead to a entanglement localization 
within a typical lenght which coincides with the localiza- 
tion lenght. 

A possible way to mimick a defect is to change the 
boundary conditions. The concurrence was studied for 
the ferromagnetic spin 1/2 XXZ chain with antiparal- 
lel boundary magnetic field which give rise to a term 
in th e Hamiltonian of th e form Hboundary = h(Si — 
5jv) (|Alcaraz et all [20041 ). The boundary field triggers 
the presence a domain walls in the system that induces 
a first order phase transition between ferromagnetic and 
kink-type phases at h c = v A 2 — 1. In the ferromag- 
netic phase the pairwise entanglement vanishes. In the 
kink-type phase the concurrence acquires a finite value 



(for h = the ground state is factorized). For a finite 
chain of lenght L and it is enhanced at the center of 
the chain. In the gapless and in the antiferromagnetic 
regimes oscillation in nearest neighbor entanglement are 
established in the system resulting from the tendency to 
reach the antiferromagnetic order. The oscillations are 
more pronounced in the gapped phase. Finally a critical 
inflection point was noticed in the m easure of Meyer and 
Wallach (|Mever and Wallachl . 12002) for the global entan- 
glement at the transition point. The spin 1/2 XXX 
antiferromagnetic chain with open bou ndary cond itions 
with single defect was also studied in (| Wand . \2004 . It 
was proved that it exists a threshold value of the cou- 
pling between the impurity and the rest of the system at 
which the concurrence between them is switched on; for 
smaller values the entanglement is dimerized in such a 
way that the monogamy property prevent the impurity 
to be entangled with the rest of the chain. 

The case of many defects was al so studied. For the 
quantum Ising model it was found ( Huang et all l2004f ) 
that the disorder can shift the point at which the con- 
currence is maximum, eventually washing out the critical 
behavior (strong disorder). The concurrence tends to be 
suppressed at the lattice site corresponding to the center 
of the gaussian; such effect is more robust for near critical 
chains. Quantum XY and XXZ chains with a gaussian 
disorder in the exchange interaction, have been also stud- 
ied to investigate how the quantum criticality of the con- 
currence is ro bust by the insertion of the inhomogeneitie s 
in the chain ( Cai et al I l2006t iHovos and Rigolinf . l2006t ) . 



B. Two and three dimensional systems 

In higher dimensions nearly all the results were ob- 
tained by means of numerical simulations. The concur- 
rence and localizable entanglem ent in two dimen sional 
quantum XY and XXZ models ( Svliuasenl . l2003al) were 
considered. The calculations are based on Quantum 
Monte Carlo (QMC) simulations an d the use of stochastic 
series expansion for spin systems ( Sandvik and Kurkiil . 
119911 : ISvliuasen and Sandvikl . 120021 ). Although the con- 
currence for the 2d models results to be qualitatively very 
similar to the one-dimensional case, it is much smaller in 
magnitude. It is the monogamy that limits the entangle- 
ment shared among the number of neighbor sites (which 
is larger in two dimensions as compared with chains). Fi- 
nally, it is observed that the maximum in the concurrence 
occurs at a position very closer to the critical point than 
in the ID case. 

By studying appropriate bounds (concurrence of as- 
sistance and the largest singular value of the connected 
correlation functions), it was proved for the XXZ model 
that the localizable entanglement is long ranged in the 
XY region up to the isotropic antiferromagnetic point. 
Similarly to the case of the quantum XY chain, the 
bounds for the localizable entanglement are very tight 
in this case. 
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The pairwise entan glement in the rf -dimensional XXZ 
model was studied in ( Gu et "oil 12005ft . The concurrence 
reaches its highest value at the antiferromagnetic quan- 
tum critical point A = — 1. A spin- wave analysis, corrob- 
orated by numerical exact diagonalization indicates that 
the concurrence develops a cusp in the thermodynamic 
limit, only for d > 2. Such behavior can be explained 
by noting that the level crossing between the ground and 
the first excited states occurring at the antiferromagnetic 
point, causes a non-analyticity in the ground state en- 
ergy. The enhanced pairwise entanglement at the an- 
tiferromagnetic point together with its non analyticity 
support the conjectured existence of long-range order for 
two dimensional antifcrromagnets. Further support to 
this conjecture is the strong size dependence of the Von 
Neumann e ntropy that beco mes singular in thermody- 
namic limit |Gu et a/.l . 1200(1 . 

The ground state entangle ment in two dimensi onal 
XYZ model were analyzed in (jRoscilde et all l2005bft by 
means of quantum Monte Carlo simulations. The diver- 
gence of the derivative of the concurrence at the contin- 
uous phase transition, observed in d = 1, was confirmed; 
also in this case the range of the pairwise entanglement 
extends only to few lattice sites. By studying the one- 
and the two tangle of the system, it was proved that 
the QPT is characterized by a cusp-minimum in the en- 
tanglement ratio t\/t2- The cusp is ultimately due to 
the discontinuity of the derivative of t\ . The minimum 
in the ratio T1/V2 signals that the enhanced role of the 
multipartite entanglement in the mechanism driving the 
phase transition. Moreover by looking at the entangle- 
ment it was found that the ground state can be factor- 
ized at certain value of the magnetic field. The existence 
of the factorizing field in d = 2 was proved rigorously 
for any 2d — XYZ model in a bipartite lattice. Unex- 
pectedly enough the relation implying the factorization 
is very similar to that one found in d = 1 . 

Multiple spin exchange is believed to play an im- 
portant r ole in the physi c s of several bidimensional 
mag nets (|Schollwock et all l2004ft . Entanglement in 
the ground state of a two leg ladder with four spin 
ring e xchange wa s eval uated by means of the concur- 
rence (|Song et aLl . l2006ft . 



C. Pairwise entanglement in fermionic models 

1. Non interacting Fermions 

The site-based entanglement of spin degrees of freedom 
through the Jordan- Wigner transformation has been ex- 
ploited for calculating the concurrence of nearest neigh- 
bor sites and the single site von Neumann entropy (see 
Section ing for the one-dimensional tight-binding model 
in presence of a ch emical potential for spinless fermions 
111 (|Zanardil . [2002T) . This model is related to the isotropic 
XX model in a transverse magnetic field. In this specific 
case, no double occupancy can occur and the concept 



of entanglement coincides with that for spins 1/2. It 
was found that the nearest neighbor concurrence of the 
ground state at T = assumes its maximum at half filled 
chain. Due to particle-hole symmetry, the concurrence 
results symmetric respect to half-filling. At finite tem- 
peratures it was found that the threshold temperature 
for vanishing concurrence is independent of the chemical 
potential. Raising the chemical potential leads from a 
monotonically decreasing concurrence with raising tem- 
perature at low filling fraction to the formation of a max- 
imum at a certain temperature for high filling fractions. 

The continuous limit of the tight-binding fermion 
model is the ideal Fermi gas. In this system, the spin en- 
tangle ment be t ween two distant particles has been stud- 
ied in (|Vedrall . [2003ft . There, depending on the dimen- 
sionality, the pairwise spin-entanglement of two fermions 
has been found to decrease with their distance with a fi- 
nite range R of the concurrence. The two spin reduced 
density matrix is 



/I-/ 2 



P12 



4-2/2 



V 




1 

-P 






-f 

1 









I-/ 2 



(41) 



where f(x) = d^- with d e {2,3} being the space 
dimension and J\ th e (spherical for d=3) Bessel func - 
tion of the first kind (jOh and Kid , |2004 IVedrall 12003ft . 
This density matrix is entangled for f 2 > 1/2. As a 
consequence, there is spin entanglement for two fermions 
closer than dn « 0.65-jr- for d = 3 and dn w 0.55-7?- for 

d = 2 (kf is the Fermi momentum). A finite tempera- 
ture tends to diminish slightly the range of pairwise spin 
entanglement (|Oh and Kid . l2004ft . 

It should not be surprising that non-interacting parti- 
cles are spin-entangled up to some finite distance. It is 
true that the ground state and even an arbitrary ther- 
mal state of non-interacting fermions has vanishing en- 
tanglement among the particles (which should not be 
confus ed with the non-vanis hing entanglement of par- 
ticles (|Dowling et all l2006bft ). since the corresponding 
states are (convex combinations of) antisymmetrized 
product states. However, disentanglement in momentum 
space typically leads to entanglement in coordinate space. 
A monochromatic plane wave of a single particle for ex- 
ample corresponds to a W state, which contains exclu- 
sively pairwise entanglement in coordinate space for an 
arbitrary distance of the sites. Furthermore does a mo- 
mentum cut-off at kf correspond to a length scale of the 
order kj . 

It is interesting that a fuzzy detection of the 
particles in coordinate space increases the entangle- 
ment det ected by the measurem ent apparatus. To 

this end ( Cavalcanti et all , 12005ft calculated the two- 

(2) 

position reduced density matrix defined by p ss , tt , = 
<* t /(r') t *t(r) t *y(r')* s (r)) with blurred field opera- 
tors 4 r s (r) := J dr' dpip s (p)D(r — r')e lpr where D(r — 
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■ exp — L V-!r is a Gaussian distribution describ- 
ing the inaccuracy of the position measurement. This 
could be understood from the blurred field operators be- 
ing coherent sums of local field operators; the entangle- 
ment measured by the apparatus as described above then 
is the bipartite entanglement between the two regions of 
width a around r and r' . This entanglement is larger 
than the average of all pairwise contributions out of it 
due to the s uper- addit ivity of the entropy/negativity. An 
analysis in (jVedra j |2004rj) for the three fermion spin den- 
sity matrix reveale d that the state c arries entanglement 
within the W-class (Di ir et a/.l . 12000). provided the three 
particles are in a region with radius of the order of the 
inverse Fermi momentum; a similar reasoning applies to 
n ferm ions in such a region (jLunkes et al 1 120051: IVedral 
l2004bD . 



2. Pairing models 

Itinerant systems, where the focus of interest is the 
entanglement of degrees of freedom forming a repre- 
sentation of su(2) in terms of the fermionic operators 
have been also subject o f intense investigat i on. T his 
line has been followed in ((Shi |2004 IZanardil . 12002ft for 
analyzing a connection to BCS superconductivity and 
also to the phenomenon of 77-pairing, a possible scenario 
for high T c su perconductivity (see also iFan et all , |2004 
IVedrail2004aU bh. Such states appear as eigenstates of the 
Hubbard model with off diagonal long range order (see 
IIH. Bp . A simplified model of BCS -like pairing fo r spin- 
less fermions has been studied in (|ZanardAl2002h . The 
concurrence of the two qubits represented by the modes k 
and — k has been found to be a monotonically increasing 
function of the order parameter; it drops to zero signifi- 
cantly before the critical temperature is reached, though. 
For electrons with spin, a connection between the BCS 
order parameter and the local von Neumann entropy in 
the particle number pro jected B CS ground sta te has been 
proposed bv IShl 120041 (see also (iGedikl . 120021) 1. 

States with off diagonal long range order by virtue of 
77-pairing are defined in Sectior HlLBl These are symmet- 
ric states and consequently, their concurrence vanishes 
in the thermodynamic limit due to the sharing prop- 
erty of pairwise entanglement of su(2) degrees of free- 
dom. Consequently, a connection to the order parameter 
of off diagonal long range order O v = { 1 77^ ?7 fc | \I/ ) = 
N(L - N)/L{L - 1) — m(l - n) (with N,L — > 00 
and fixed filling fraction n) can not be established, not 
even for the rescaled concurrence, since C — > 1/L (see 
also the analysis for the LMG model in Section [iV.A.2[) . 
Nevertheless, the state is entangled, as can be seen from 
the entropy of en tanglement and the geom etric measure 
of entanglement (jWei and Goldbartl 120031) . The latter 
is tig htly connected to the relative entr opy dWei et all , 
12004 ) . Both have been calculated in (jVedrall . l2004af ) 
and clearly indicate the presence of multipartite entan- 
glement. 



3. Kondo models 

The Kondo models are paradigms to explore the quan- 
tum impurity problems. They identify a special class 
of physical systems whose macroscopic properties are 
dramatically influenced by the presence of fe w impuri- 
ties w ith quantum internal degrees of freedom ( Hewsor] . 
1997). In its simplest formulation, the effective Hamil- 
tonian describes a single impurity spin interacting with 
a band of free electrons. The many-body screening of 
the impurity spin provided by the electrons as a collec- 
tive effect, leads t he system from weak coupling to strong 
coupling regimes <|Andreil . fl995h : the lenght scale of the 
screening cloud is £k — v/Tk, where v is the speed of the 
low lying excitations and Tk is the Kondo temperature. 

In the first studies of the entanglement the charge 
degrees of freedom of the electrons were considered 
frozen. The pairwise entanglement of spin degrees 
of fre edom in the isotr opic Kondo model was ana- 
lyzed (lOh and Kiml . l2006l) within the variational formal- 
ism of Yosida where the Kondo singlet is described as 
= 75 (1^)1X1 ) - \4>\)\X\))- In the previous equa- 
tion \xa) denote the impurity spin states; \4> a ) represent 
the e lectronic states with an unbalanced spin a ( Hewsonl . 
119971) . In agreement with the common wisdom, the re- 
duced density operator of the impurity is found to be 
maximally mixed, meaning that the Fermi sea and the 
impurity spin are in a maximally entangled state (the 
Kondo singlet). The impurity spin and a single electron 
are in a Werner state made of a superposition of the 
back ground and the Kondo singlets. Due to the entan- 
glement monogamy (the electrons cooperatively form a 
singlet with the spin) two electrons cannot be entangled 
with each other within the Kondo cloud and the single- 
electron spin entanglement vanishes in the thermody- 
namical limit. Pairs of electronic spin can be nevertheless 
entangled in a finite system through the scattering with 
the spin impurity; this effect might be used to manip- 
ulate the electron-electron entanglement by performing 
a pro jective measure on the impurity spin ( Yang et all , 
2005). This suggests that some amount of electron- 
electron entanglement might be extracted even in the 
thermodynamical limit where it was demonstrated that 
the K ondo resonance i s wash ed out by the measure- 
ment ( Katsnelson et all , l2003l )-effectivelv removing the 
constraint of the entanglement sharing. 

The two impurity Kondo model was studied as well. 
The new feature here is the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) effective interaction between the impu- 
rity spins Si and S2, that competes with the Kondo 
mechanism (favouring non magnetic states); it is fer- 
romagnetic or antiferromagnetic depending on the dis- 
tance between the impurities. Because of such inter- 
play a quantum critical point emerges in the phase 
diagram separating the spin-spin interaction regime 
from the phase where the two spin are complete ly 
screened (jAffleck and Ludwid . Il992t bones et all , Il988l ) . 

As for the single impurity, the two impurity spins are 
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in a Werner state, for which the concurrence is character- 
ized by a single parameter p s , exhibiting a singlet type 
of entanglement between the two spin impurities. The 
concurrence is found to vanish at the critical point. 

For ferromagnetic RKKY interaction the concurrence 
between the impurity spin vanishes identically as the 
result of a S = 1 Kondo sc reening. It turns out 
that ( Cho and McKenzi e. 2006) the impurity spins can 
be entangled (with a finite concurrence) by the RKKY 
interaction only when certain amount of antiferromag- 
netic correlations f s = (S1S2) is established in the sys- 
tem; such value of the correlation function is that one 
reached at the quantum critical point. The entanglement 
between the conduction electrons and the Kondo impu- 
rities is quantified by a combined analysis of the Von 
Neumann entropies of the two impurities and of the sin- 
gle impurity (tracing out both the electronic spins and 
the remaining impurity). The latter quantity is maxi- 
mized independently on f s , meaning that the impurity 
spin is completely screened either by the Kondo cloud 
or by the other impurity spin. In the regime where the 
Kondo mechanism dominates, the concurrence cannot be 
finite because of the entanglement sharing. 

Entanglement in the Kondo physics of double quan- 
tum <fots_Jj2_an external magnetic field was studied 
in ( Ramsak et Ml, l2006h . The main phenomenology 
res ults to be consistent wi th the scenario depicted 
in ( Cho and McKenzid . 120061 ) especially if the dots are 
arranged in series (each dot is coupled to the leads exclu- 
sively, resembling the configuration of the Kondo spins 
embedded in the electrons). The concurrence switches 
to finite values for a certain threshold of the inter-dot 
coupling (for which the assumptions of negligible charge 
fluctuations results still valid) . The temperature weakens 
the entanglement between the qubits also at T > Tk ■ For 
the side- and parallel-coupled dots a more intense cou- 
pling amoung the qubits is required to entangle them. 
For the side-coupled dots this results because one has to 
win on the enhanced Kondo effect on the dot coupled to 
the leads ('two stage Kondo' effect); therefore the criti- 
cal inter-dots coupling is ~ T K . For the parallel-coupled 
arrangement the concurrence is zero because the effec- 
tive RKKY interaction turns out to be ferromagnetic up 
to a certain value that is the threshold to entangle the 
electrons. 

The RKKY interaction controls in an effective way 
the entanglement amoung the qubits also in the case 
of many im purity spins arranged as in the Kondo neck- 
lace model ( Saguia and Sarandvl . l2003h . The Hamilto- 
nian describes a Kondo lattice where the localized im- 
purity spins, displaced in every lattice site, interact with 
the (p seudo) spins of the electrons (see iTsunetsugu et all 
1997 for a review on Kondo lattice models). It results 
that the additional on-site spin-spin interaction impose 
a 'selective' monogamy of the entanglement, depending 
on whether the Heisenberg interactions is ferromagnetic 
or antiferromagnetic. The effects of finite temperature 
and magnetic field were considered also in anisotropic 



models. It emerges that a critical field exists separating 
different patterns in the thermal entanglement between 
the eigenstates of the model. 

The effect of the fluctuations of the charge de- 
grees of freedom of the electrons (frozen in the ref- 
erences cited above) in the RKKY mechanism was 
discussed quantifying the entanglement of particles in 
a small cluster described by the per iodic Anderson 
model ( Samuelsson and Verdozzil . 120061 ). It was evi- 
denced that the ground state of the system is charac- 
terized by a double occupancy of the electronic levels, 
whose entanglement can be only partially captured by 
assuming them as simple qubits. 



D. Entanglement in itinerant bosonic systems 

In contrast to the free Fermi gas, in bosonic systems 
the phenomenon of Bose-Einstein condensation (BEC) 
takes place at sufficiently low temperatures. Then, a 
macroscopic portion of the bosons is found in the single 
particle ground state of the system. This state clearly 
is a symmetrization of a product state in momentum 
space, and in fact all eigenstates of the ideal Bose gas 
are of that structure. Nevertheless, in principle entan- 
glement could be present when going in the coordinate 
basis. It results, however, this not the case. Neither 
do tw o distant bosons carry spin entanglement ( Vedrall . 
2003) nor are two distant groups o f n and m parti- 
cles entangled ( Dowling et all l2006bf ) when superselec- 
tion rul es for the particle number i s applied i n both 
regio ns ( Wiseman and Vaccarol l2003h . In Ref. ( Simonl . 
I2002D a very different notion of entanglement has been 
employed: the state of either part of a certain biparti- 
tion of the BEC has been viewed at as a qudit, or more 
precisely, a iV-level system. TV is the number of bosons 
in the condensate and the different states in both re- 
gions are labelled by its occupation number. The entropy 
of entanglement for a spatial bipartition of the BEC is 
then non-zero. A proposal for entangling internal atomic 
degrees of freedom in a weakly interacting BEC has 
been put forward in Ref. ( Helmerson and L.Youl . l200lt 
ISprensen et~al\ , \200l\ ) . 

The disentangled modes in a BEC naturally become 
entangled by means of interactions /scatt ering between 
these mo des. This h as been exemplified in (jVedrall . l2003f) 
(see also HE I2004D for the case where the scattering 
strength is independent of the momentum transfer q: 

H = J2p £ (p) a l% + v J2p.p'. q a l+ q a l'- q a P ' a p The Hamil- 
tonian can be diagonalized with for q — and p' = —p by 



means of a Bogoliubov transformation a p —: u p b p +Vpb_ p . 
We observe that such transformation entangle the two 
modes ±p. The cor responding entropy of entangle- 
ment is jvjkjl 111 Sp,^ = -(^) 2 ln(^) 2 - (1 - 

(^) 2 ) ln(l — (^-) 2 ) If Up = v p , the reduced state for mode 
p is maximally mixed, and hence the modes p and — p are 
maximally entangled. The entanglement entropy for a bi- 
partition in positive and negative modes is then given by 
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_ p . This constitutes a simple example on how 
mode mixing generates entanglement. Such a scenario is 
rather generic; a curious example being the entanglement 
of the accelerated vacuum and from the viewpoint of rel- 
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Further studies of entanglement in boson ic system in- 
clude the analysis in two- mode condensates (iHines et all 
|2002| ). in optical lattices ( Ng and Burnettl. 120061) and i n 



two-species spinor Bose condensates (|Shi and Niul .l 2006t) . 




B 



E. Entanglement of particles 

Studies which use measures for indistinguishable parti- 
cle entanglement (see Section III.F[) in the area of many- 
body systems are is still only few, particularly regard- 
ing the use of the fermionic concurrence, giving account 
for the possibility of double occupancy (with spin de- 
gree of freedom). The entanglement of particles and 
its differen ce with the usual spin entanglement, is dis- 
cussed in ( Dowling et all l2006bl ). starting with very 
small systems as two spinless fermions on four lattice 
sites and the Hubbard dimer, and then for the tight bind- 
ing model in one spatial dimension, in order to compare 
with exi sting results f or the spin entanglement elabo- 
rated in ( Vedrall . l2003h . For the Hubbard dimer (a two- 
site Hubbard model), the authors compare with their re- 
sults with those for the entanglement measured by the 
local von Neumann entropy w ithout supersel ection rule 
for the local particle numbers (jZanardil . 2002). Whereas 
the latter signals decreasing entanglement in the ground 
state with in creasing U /t, the ent anglement of parti- 
cles increases (jDowling et aZ.I . l2006bh ■ This demonstrates 
that imposing superselection rules may lead to qualita- 
tively different behaviour of the entanglement. Inter- 
estingly, an increase with U/t is observed also for the 
entan glement of modes without imposing superselection 
rules (|Deng et aLl . l2005l ). 

We would like to finish this section with the notice 
of a recent proposal of an experiment in order to decide 
whether even entanglement merely due to the statistics of 
the indistinguishable pa rticles can be use f ul for quantum 
information processing (jCavalcanti et "aE l200fjh . 



FIG. 6 The block entropy is evaluated after partitioning the 
system in regions A and B. For finite range correlations, it 
is intuitive that the wave function of the system is factorized 
\^a)\^b) by removing the region at the boundary (yellow). 
Accordingly the reduced entropy would vanish 



Wilczek ()Holzhev et all [l994) born in the context of 
black hole physics, the problem was first reanalyzed in 
the framework of quantum informa tion by Vidal and 
coworkers for quantum spin chains ( Vidal et all l2003f ) 
and by Audenaert and co workers for harmonic lat- 
tices ( Audenaert et aLl . l2002f h 

In studying the properties of block entropy it is im- 
portant to understand its dependence on the properties 
(topology, dimensions, • • • ) of the two regions A and B. 
A key property which is intensively explored to under- 
stan d its range of v alidity is what is known as the area 
law ( SrednickH 1 1 993f h When it holds the reduced entropy 



S would depend only on the surface of separation between 
the two regions A and B. In d-dimensional system this 
means that S ~ £ d_1 where I is of the order of the size 
of one of the block (see iRiera and Latorrel . 120061 for a re- 
cent discussion of various aspects of the area law). In 
the rest of the Section we discuss several different phys- 
ical systems in one- and higher dimensional lattices and 
see when the area law hold. Here we consider only many- 
body systems in their ground state, the thermal effect and 
the dynamical properties of the entropy will be discussed 
in Sections IVIII and |IX] respectively. 

The entanglement entropy is not a mere theoretical 
concept but it mig ht be measured. Following the pro- 
cedure proposed in ( Klich et~al\ , l200rj ) the measure of S 
can be related to the measurement of the distribution of 
suitably chosen observables. 



V. ENTANGLEMENT ENTROPY 



An important class of works analyzing the entangle- 
ment in many body system considered a bipartition of 
the system dividing it in two distinct regions A and B 
as shown in Fig|6] If the total system is in a pure state 
than a measure of the entanglement between A and B is 
given by the von Neumann entropy S associated to the 
reduced density matrix of one of the two blocks (pa/b) 

M otivated by the pioneering work of Fiola et 
al. ( Fiola et all Il994l ) and Holzhey, Larsen and 



A. One-dimensional spin systems 

We start our review on the properties of the block en- 
tropy by analyzing the case of one-dimensional spin sys- 
tems to which a large body of work has been devoted. By 
means of the Jordan- Wigner transformation it is possi- 
ble to map the models onto a lattice fermions, hence the 
results discussed here are applicable to fermionic models 
(after the appropriate mapping) as well. A particular im- 
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portant case which is amenable of an exact solution is the 
XY model (see Section IIII.A.1|) which can be mapped 
onto a free fermion model. For this case we discuss in 
more details the method to calculate the block entropy. 
In this section we consider only chains with short-range 
interaction. 



1. Spin chains 

In one dimension the surface separating the two re- 
gions is constituted by two points, therefore the area's 
law would imply that the reduced entropy is indepen- 
dent on the block size. This is indeed the case when 
the system is g apped and hence the correlation lenght 
£ is finite (see iHastingsl [2004 for a rigorous proof). 
In the gapless case, £ = oo, logarithmic corrections 
appear and the prefactor of the block entropy is uni- 
versal, related to the central charge of the underlyin 
confo rmal field theory. Holzey et al. ( Holzhev et al. 
Il994f ). b enefiting from an e arlier work of Cardy and 
Peschel ( Cardv and Peschel Il988l ). analyzed the block 
entropy of a 1 + 1 dimensional massless bosonic field. 
Vidal et al. studied n umerically one-dimensional Ising 
and Heisenberg chains ( Latorre et all [20041 IVidal et all 
I2003D and conjectured that the block entropy would sat- 
urate for non critical chain while would diverge logaryth- 
mically with a prefactor related t o the central charge o f 
the underlying conformal theory (jHolzhev et al .. 1994). 
Such violation of the area law in critical systems reflects 
how the mixedness of the state increases by the partial 
tracing operation, regardless of the spatial extension of 
the spin block. 

A calculation of the block entropy by means of con- 
formal field theory, generalizing in several respect the 
the results of Holzey et al. by including the case of free 
and periodic boundary conditions, different partitions, 
non critical systems and finite tem perature has been per- 
form ed by Calabrese and Cardy (|Calabrese and Cardvl . 
12004). Starting fr om the work on the XX model 
of Jin and Korcpin ( Jin and Korepinl 1200 4b). im- 
portant explicit analytic calculations for a number 
of one-dimensional XY spin ( free fermion) mod- 
els have been carried out in (lEisler and Zimboras 
20051: iFranchini et all l2006allbl lite et all 120051 12006 



Jin and KorepmC" 2004bt iKeating and Mezzadril 12005 



Korepinl. 12 004: Peschel, 2004, 2005; Popkov and Salerno . 



2005 : Westonl . 20061 ). Numerical calculations on 
the XX and XXZ mo d els were a l so p erformed 
in (|De Chiara et all 120051 : iLaflorenciel l2005fh The 

is of great inter- 



of 



est in a more general 


context [see Casini et al.. 


2005; Casini and Huerta. 


2005a; Fiola et all 1994: 


Holzhev et al.. 1994: Rvu and Takavanagi, 2006a.lb and 



references therein], Casini and coworkers, for example, 
evaluated the entangle ment entropy both in the mas- 
sive scalar field theory (|Casini and Huertal . l2005aj ) and 



are a law for gapped one dimensional systems was proved 
by (|Hastingsl . l2007al lbT). 

The main features of the reduced entropy in one- 
dimensional spin/fermi systems can be summarized as 
follows (for clarity we discuss only the long distance be- 
haviour as dictated by the underlying conformal field the- 
ory): 

• At criticality a one dimensional system has a block 
entropy which diverges logarithmically with the 
block size. If the block is of lenght I and the system 
is L long with periodic boundary condition then Se 
is given by 



St = g lOg 2 



L sir 
— sin I — I 



(42) 



where c is the central charge of the underlying con- 
formal field theory and a is an ultraviolet regular- 
ization cut-off ( for example the lattice spacing in 
spin systems). A is a non-universal constant. For 
the Ising model c = 1/2 while for the Heisenberg 
model c = 1 (see Fig. 7). 

• Slightly away from criticality, in the case in which 
the system has a large but finite correlation lenght 
£ ^> a and the Hamiltonian is short-ranged, the 
block entropy saturates to a finite value (see Fig. 7) 
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• An extension to finite temperature in the critical 
case ha s been obtained by means of conformal field 
theory ( Calabrese and Cardvl . [20041 ) and by confor- 
mal mapping together with t he second law of ther- 
modynamics (Korepin, 2004J) with the result 



Si = g lOg 2 



sinli ( —I 



A 



(44) 



for Dirac fields ( Casini et al. . 20051 )). Very recently the 



where A is a constant and j3 the inverse tempera- 
ture. In the finite temperature case however, the 
block entropy is not a measure of the entanglement 
between the two partitions as the state to start with 
is mixed. A comparison of numerical data with the 
CFT predictions is shown in Fig. 7 for the XXZ 
model. 

By now it is clear that the various measures of en- 
tanglement are sensitive to the presence of quantum 
phase transitions, the scaling of the entropy gives ex- 
cellent signatures as well. Recent works tried to con- 
struct efficient ways to detect quantum phase transi- 
tions by analyzing the reduced entropy for s mall clusters . 
One- s ite entropy has been considered in ( Chen et all 
120061 : iGu et all 120041) . Chen et al analysed the en- 
tanglement of the ground states in XXZ and dimer- 
ized Heisenberg spin chains as well as in a two-leg spin 
ladder suggesting that the phase boundaries might be 
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FIG. 7 The Block entropy Se for L = 200 in the case of an 
XXZ Heisenberg chain for a critical value A = 0.0 (circles) 
and non-critical value A = 1.8 (squares). The critical data 
compared with the conformal field theory prediction (dashed 
line). Lower inset: central charge extrapolated by fitting the 
numerical data Se with Eq. (142[) (with a factor 1/2 as in the 
numerical calculation the block is taken at the boundary) for 
different values of A. The data are for L = 1000. Upper 
inset: scaling of c extrapolated as a function of 1/L for the 
worst case A = 0.5 and compared to a quadratic fit (dashed 
line) . [From (|De Chiara et all 120051 )] 
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contiguous but A consists of I equally spaced spins, such 
that the spacing between the spins in this subsystem is p 
sites on the chain while B contains the remaining spins. 
For this choice of the partition, the "surface" separat- 
ing the two blocks grows with the system size (differ- 
ently from the case p — 1 where it is composed by two 
links). As a consequence non-local properties of entan- 
glement between the two blocks can be investigated. For 
p > 1 the leading contribution to the entropy scales lin- 
early with the block size Se(p) = £i(p)l + £"i(jp) ln£ In 
the limit p — » oo the coefficient £\{p) is, to leading or- 
der, sum from single spin contributions. The unexpected 
result is that the corrections vanish slowly, as l/p, differ- 
ently from other measures like concurrence where these 
long-range corrections are not present. 

We conclude this section by briefly discussing the 
single c opy entanglement E\ intr oduced by Eisert and 
Cramer ( Eisert and Cramerl 120051 ) a nd studied in details 
for one-dimensio n al sp i n systems ( Eisert and Cramerl . 
120051: lOrus et all 120061 : iPeschel and Zhaol 120051) . Sin- 
gle copy entanglement Ssc,e quantify the amount of en- 
tanglement that can be distilled from a single speci- 
men of a quantum system s. For spin chains it can be 
shown ( Eisert and Cramerl . 120051 ) that single copy en- 
tanglement asymptotically is half of the entanglement 
entropy lini£_ 1 . 00 Ssc,t/St = 1/2. This result was later 



generalized to conforma l ly inv ariant models (jOrus et all 
120061: IPeschel and Zhaol I2005D . 



2. XY chains and free fermion models 



FIG. 8 The idea of the comb-partition suggested by Keat- 
ing Mezzadri and Novaes is illustrated for three different val- 
ues of the spacing p. In all the case the subsystem A, de- 
noted by black circles) contains I spins while the other (de- 
noted by empty circles) contains the rest of the chain. The 
case, p — 1, corresponds to t he well known 'block' divi- 
sion. [From l|Keating et all l2006ft ] 



identified based on the ana lysis on the l ocal e xtreme of 
the entangleme nt entropy (fchen et all 20061). Legeza 
and c oworkers ( Legeza and Solvoml 20061 :" Legeza et all 
120061 ) pointed out that in the biquadratic spin-1 Heisen- 
berg chain, see Eq.(|27|). the two-site entropy is ideal to 
highlight the presence of a dimerized phase. They also 
considered t he two-site entropy also for t he ionic Hub- 
bard model ( Hubbard and TorranceL [T"98lh - 

The idea of partitioning the system in a more elab- 
orated way in order to analyze additional properties of 
entanglement lead t o the introduct i on of the concept of 
comb entanglement ( Keating et aLl .l2006). This is illus- 
trated in Fig(5) The two blocks A and B are not chosen 



As we saw in Section IIII.A.l.al the XY chain (A = 

in Eq. (|2li|) ) can be mapped onto a model of free 
fermions. As a result an analytical (albeit non-trivial) 
approach is possible for the calculation of the block en- 
tropy. An analytic proof of the logarithmic dependence 
of the block entropy in the isotrop ic XY mod e l was 
obtained in (iFranchini et all l2006al lbl: llts et all 120051: 

1 Jin and Korepinf 2004b ) . The relation between the en- 
tanglement entropy of this model and the corner transfer 
ma trices of the tr iangular Ising lattice has been derived 
in ( Peschell |2004 ). Keating and Mezzadri considered a 
more general free fermion Hamiltonian in which the ma- 
trices A and B (defined in Eq. (|2"4")) ) do not have the tridi- 
agonal structure which appear in t he case of the XY 
model ( Keating and Mezzadril . [2005) . They showed that 
under certain conditions the entropy can be expressed 
in terms of averages over ensembles of random matrices. 
In this section we recall the main steps of the derivation 
leading to the evaluation of the entanglement entropy , 
more detail s can be found in th e review of Latorre, Rico 
and Vidal (lLatorre et all I2004D and in the above men- 
tioned papers (IFranchini et all l2006allbl: llts et al. 1. 120051 



I Jin and Korepinll2004rilKeating and Mezzadril |2005| ). 

The reduced density matrix of a block of I spins can 
be expressed in terms of averages of strings of I spin 
operators with weights given by the averages of these 
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strings on the ground states. By means of a map- 
ping of the spin operators in terms of the Majorana 

fermions, a 2 ^i = (Il m <j CT m) °f ; a 2/ = (II m <z CT m W 
and given the fact that the resulting fermionic Hamilto- 
nian is quadratic (Wick theorems holds), it is possible to 
express the block entropy in terms of the elements of the 
correlation matrix Bg 
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for L — > oo, by 



with real coefficients gi given as, 
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(46) 



It is crucial to notice that the matrix Bg are block Toeplitz 
matrices, that can be thought as usual Toeplitz matrices 
but with non commuting entries. After the transforma- 
tion of Bi into a canonical form = v m ia y , the 
system is described by a set of i independent two-level 
systems. Therefore the entanglement entropy is given by 



rn — 1 
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■log 
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1 - v„ 



(47) 

Numerical and analytical analysis of (|47|) has been per- 
formed leading to the behaviour described in the preced- 
ing subsection. In order to obtain the analytical formula 
for the asymptotics of the entropy ([47)) , the first step is to 
recast it, by the Cauchy formula, into a contour integral 
(|lts et all 120051 : Ijin and Korephl . l2004bh 



s = hm lim -L/ e (l + 6,A) dln ^ A(A) dA, (48) 



where e(2x, 2y) — — ^2±(xiy) log 2 (x ± y). The contour 
of integration c(e,S) approaches the interval [—1,1] as e 
and S tend to zero without enclosing the branch points 
of e(l + e, A). The matrix in D A (\) = det(A7 - Bt) is 
again of the block Toeplitz type. The asymptotics of 
the entropy can then be obtained from the asymptotics 
of Da (A). This, in turn can be done resorting to the 
Riemann-Hilbert approach to the theory of Fredholm in- 
tegral equations ( Bogoliubov et aLl .[T993). This allows to 
generalize the standard Szego theorem for scalar Toeplitz 
matrices to obtain the leading formula for the determi- 
nant of the block Toeplitz matrix Da(X). This leads to 
the asymptotics, I — > oo, for the entropy 



M«.t) + f ]9«[«i) 
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dx (49) 



where 9z{s,t) is one of the Jacoby elliptic function, 
(3{x) = (l/27ri)[ln(x + 1) - ln(x - 1), a = 1,0 for h z < 



, > h c and r ~ [log(| 1 — A _1 |)]. The critical behaviour 
of Eq. (|4*5|) can be obtained by the asymptotic properties 
of 03(s,t) for small r, and the leading term of the crit- 
ical entropy results S = -(1/6) log (|1 - A" 1 |/(47) 2 ) + 
0[\1- A -1 1 (log |1 - A" 1 1) 2 ] . The previous expres- 
sion can be obtained in a more direct way resort- 
ing a duality relation connecting the quantum Ising 
chain with the Ising model on a square or tr i angu- 
lar lattice (ICalabrese and Cardvl . 12004 iPeschell 12004 
iPeschel et all 1 19991) . In particular the reduced density 
matrix can be written as the trace of the corresponding 
(corner) transfer matrix. In this way, however the ex- 
pression for the critical entropy can be obtained only for 
j 2 + h 2 z >l. 

In the isotropic case 7 = 0, c = 1 and then the prefac- 
tor of the log divergence is 1/3. In this case the critical 
entropy can be obtained as an average, by realizing that 
the block Toeplitz matrices (for 7 = 0) are unitary. Then 
the contour integral can be recast into an integral in the 
ensemb le of unitary matrices . Thi s remarkable obser- 
vation ([Keating and Mezzadril |2005|) allows to relate the 
spectral statistics of the model with the entanglement en- 
coded in the ground state, following a reasoning that has 
many analogies in spirit with random matrix theory. The 
entanglement entropy was obtained explicitly for matri- 
ces B, in Eq. ([24|) . being elements of classical groups. It 
is interestingly that only B affects the the prefactor of 
the logarithm in entropy. It is proportional to 2 WG where 
wq is a universal quantity related solely to the classical 
group establishing the symmetries, the constant of pro- 
portionality is Hamiltonian dependent. 



3. Disordered chains 

Conformal invariance implies universal properties for 
the entanglement entropy. What happens when con- 
formal invariance is lost as in the case of certain 
one-di mensional disordered spin systems? Refael and 
Moore (Refael and MooreL 12004) were the first to look 
at this question by computing the block entropy for the 
Heisenberg, XX, and quantum Ising chains with random 
nearest-neighbor coupling. Their approach was based 
on a real space renormalization group developed ear- 
lier (|Fisherl . ll99l lM a et al. , 1979) for random spin chains 
where disorder is relevant and drives the system at low 
energies in the so-called random singlet phase which can 
be thought as a collection of singlet bonds of arbitrary 
length. Consequently the entropy of a given segment of 
the chain is just In 2 times the number of singlets cross- 
ing the boundary between the two regions in which the 
systems is partitioned. Refael and Moore showed that 
the entropy, as in the case of clean critical chain, grows 
as the logarithm of the block size Si ~ c log€ with a 
"renormalized central charge" c = cln2. 



A numeric al test of this 
for the XX (lLaflorencie 
models (De Chiara et al 



predi ction was performed both 



2005) and for the Heisenberg 
2005). In Figdwe report the 
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data of Lafiorencie, the two curves represent the results 
for the clean and disordered case and fully confirm the 
prediction of Refael and Moore. The result that the ra- 




FIG. 9 (color online) Entanglement entropy of a subsystem 
of size I embedded in a closed ring of size L, shown vs I 
in a log-linear plot. Numerical results obtained by exact 
diagonalization performed at the XX point. For clean non 
random systems with L — 500 and L = 2000 (open cir- 
cles), Se is in agreement with Eq. (1421 (red and blue curves). 
A = 0.8595 + i|2lni, Y = 0.72602 + | ln(^ sin ^ ), 
Z = 0.72567 + 1 ln(52£ sin [from (jLaflorencid . 120051 )] 



tio between the random and pure values of the prefac- 
tor of the block entropy is the same for a ll the different 
chains studied in (jRefael and Moord . [20041 might suggest 
that this value might be determined by the flow from the 
pure to the random fixed point. This conjecture was re- 
cently questioned by analyzing the entanglement entropy 
for a family of models which includes the N-states ran- 
dom Potts chain and the Zn clock model. In this case it 
was shown that the ratio between the entanglement en- 
tropy in th e pure and in the d isordered system is model 
dependent ( Santachiaral . 1200(1 ). 



4. Boundary effects 

Boundaries or impurities may alter in a significant way 
the the entanglement entropy. 

The result given in Eq.([42| was obtained for pe- 
riodic boundary conditions. If the block is at the 
boundary of the chain then the prefactor is modified 
and the block entropy is one half of the one given in 
Eq.(|42} (ICalabrese and Cardvi . l2fK)l . 
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where A is the non-universal const ant given in Eq. 
and q is the boundary entropy ([Affleck and Ludwiel . 
Il99lh . The case of open boundary conditions in 
cri tical XXZ chains wa s also recently reconsidered 
in (Lafior encie et al ., 200G). In addition to the log diver- 
gence there is a parity effect depending on the number 
of spins of the block being even or odd. The amplitude 



of the resulting oscillating term decays as a power law 
with the distance from the boundary. The origin of this 
oscillating term is easy to understand qualitatively as an 
alternation of strong and weak bonds along the chain. 
The boundary spin has a strong tendency to form a sin- 
glet pair with its nearest neighbour on the right-hand 
side; due to the monogamy of the entanglement this last 
spin will be consequently less entangled with its part- 
ner on the third site of the chain. Furthermore it was 
also shown that the alternating contribution to the en- 
tanglement entropy is proportional to a similar term in 
the energy density (the constant of proportionality be- 
ing related to the lattice constant and to the velocity of 
the excitations) . The effect of open boundary conditions 
on the entanglement entr opy of a resonant valence bond 
solid was studied as well ( Fan et al\ , 2006). In this case 
however the corrections due to the open ends decay ex- 
ponentially. 

Different type of boundaries can appear in the AKLT 
quantum spin chain, with bulk spin-1 and two spin-1/2 
at the ends. The en tanglement entropy has been studied 
in ( Fan et a/] . [2004). They showed that the block entropy 
approaches to a constant value exponentially fast with i 

The entanglement entropy of one-dimensional sys- 
tems is affected by the pr e sence of impurities in the 
bulk (jLevind . 12004 [Peschell 120051: IZhao et all l2006f ) or 
aperiodic couplings ( Igloi et all 2007 ). In these cases the 



entanglement entropy has the same form as in Eq. ([42|) 
but with an effective value which depends on the strength 
of the defect. The entanglement properties of anisotropic 
open spin one-half Heisenberg ch ains with a modifi ed cen- 
tral bond were considered in ( Zhao et~al\ , 120061 ) where 
the entanglement entropy between the two half-chains 
was calculated using the DMRG approach. They find 
a logarithmic behavior with an effective central charge 
varying with th e length of t he sy stem. The numerical 
simulations of ( Zhao et al .], 120061) show that by going 
from the antifcrromagnctic to the ferromagnetic case the 
effective cen tral charge gr ows from zero to one in agree- 
ment with ( Levind . l2004h . The combined presence of 
interaction between the excitation and a local impurity 
modifies in an important way the properties of a one- 
dimensional sy^tejn^_Stardng_from the work of Kane 
and Fisher ( Kane and Fisheiill992l ) is by now understood 
that at low energies the scattering with the impurity is 
enhanced or suppressed depending on the interaction be- 
ing repulsive or attractive. It is therefore expected that 
the entanglement entropy is affected as well. Levine, 
by means of bosonization, studied the entanglement en- 
tropy in a Luttinger liquid interrupted by an impurity 
and found that there is a correction to S, due to the im- 
purity which scales as SSi mp ~ — Vb log (I / a) where Vb i s 
the renormalized backscattering constant (jLevind . 120041 ). 
In the repulsive case the backscattering flows to large val- 
ues suggesting that the total entropy would vanish (the 
correction is negative). In the opposite case of attractive 
interactions, the impurity potential is shielded at large 
distances and the entropy would approach the value of 
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the homogeneous liquid. 

The single copy entanglement i n the presence of b ound- 
aries has been considered as well ( Zhou et flZl l200rfh Dif- 
ferently from the bulk contribution here the boundary 
contribution to the von Neumann entropy equals that of 
the single copy entanglement. 

Some of these results provided the fertile ground to 
study the entanglement encoded in the Kondo cloud. 
Specifically, the block entropy Si mp of a spin cloud of 
radius r aro und the impurity with the rest of the system 
is analyzed (jSprensen et all [20061 ). By using a combina- 
tion of Bcthc ansatz results, conformal field theory and 
DMRG methods, the authors demonstrated that Si„ lp is 
a universal scaling function of r/^K- 



B. Harmonic chains 

Static systems of harmonic ch ains have been first ana- 
lyzed in ( Audenaert et a/.l . l2002f ). where periodic arrange- 
ments of harmonic oscillator modes have been considered. 
The oscillators have been coupled in the standard way via 
their coordinate variables and the Hamiltonian has been 
chosen to be translational invariant. The entanglement 
in these systems has been analyzed for both the ground 
state and of thermal states; both belong to the class of 
G aussian states. Here w e review the results obtained 
in lAudenaert et all 120021 . For work on higher dimen- 
sional lattices and emphasize on the entropy area law we 
refer to Section IVXI21 

Using the covariance matrix defined in Eq. (|35[) . the log- 



to a symplectic eigenvalue of V that diverges with the 
system size with a consequent divergence of the negativ- 
ity. 



arithmic negativity can b e expressed direct l y in t erms of 
the interaction matrix U (lAudenaert et all 120 02) En = 



tr log 2 U _1//2 PU 1 / 2 P where P is a diagonal matrix, with 
non-zero entries Pjj = — 1 where the partial transposition 
is performed and Pjj = 1 elsewhere. This entanglement 
monotone has been analyzed for bipartitions of a ring 
containing an even number of oscillators. It is convenient 
f U' U"\ 

to define U = I jj„ jj, I For the symmetric bisection 

into equally large connected parts, a lower bound for the 
logarithmic negativity has been obtained as 

En > ||tr F n log 2 U| = \ log 2 (l + 4 £ a 2m+1 ) (51) 

m=0 

where the coefficients a n have been defined in Sec- 
tions IIII.DI F n is the n x n flip matrix with "1" in 
the cross-diagonal and "0" elsewhere. Equality holds 
if FikU" is semi-definite, which is the case for nearest 
neighbor interaction. For this case one obtains E N ' n ' = 
|log 2 (l + 4a i) . Remarkably, this result is independent 
of the size of the ring. This also tells us that the neg- 
ativity of the symmetric bisection for a model including 
couplings ad of arbitrary range is higher than that of the 
corresponding chain with only nearest neighbor coupling 
and coupling strength Yld a d- ^ is interesting to antici- 
pate here that for critical systems, the lowest eigenvalue 
of U tends to zero with growing system size. This leads 




FIG. 10 Logarithmic negativity for harmonic chains. Top 
panel: The logarithmic negativity for coupling strength a = 
20 and a bipartition in m and rii oscillators (total number 
of oscillators N = n\ + n^). For sufficiently large parts, a 
plateau is reached. For one part consisted of only few oscil- 
lators, the negativity decreases as a function of the system 
size. Bottom panel: The logarithmic negativity relative to its 
plateau value En (oo) for ri2 = 20 as a function of n\ and 
varying coupling strength a. For small coupling the plateau 
limit is reached faster; since the plateau limit is connected to 
a local quantity (the average energy per oscillator), this can 
be this can be explained by a correlat ion length that grows 
with a. [From ijAudenaert et all l2002h ] 



The analysis for general bisections revealed that for 
nearest neighbor couplings the negativity of a single os- 
cillator with the rest of the chain monotonically decreases 
with the size of the chain. This single-oscillator negativ- 
ity turned out to establish also a lower bound for the 
negativity of any connected set of two or more oscilla- 
tors with the rest of the same chain (see top panel of 
Fig. UOp . In all cases the maximum negativity has been 
observed for the symmetric bisection. Both features are 
expected t o be generic to coupled e nsembles of harmonic 
oscillators (] Audenaert et aLl . I2002T ). In particular should 
the infinite size limit of the symmetric bipartition nega- 
tivity establish an upper bound: limm—Kx, En (m, m) > 
EN{ni,n-2). This upper bound appears as a plateau in 
the top panel of Fig. [TUl which as a function of m and 
n 2 is reached already for not too small n\ and n%. The 
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plateau value is essentially proportional to the average 
energy per oscillator. With increasing nearest neighbor 
coupling strength, a more shallow approaching of the 
plateau value is observed (see bottom panel of Fig. [TU)) . 

The situation changes when the negativity of two dis- 
connected parts of the chain is considered. The particu- 
lar limiting case of an alternating bipartition, consisting 
in all the oscillators located at odd sites and the com- 
plement being all the oscillators at even sites, has been 
analy zed in the presence of nearest neighbor couplings 
only ( Audenaert et all 120021 ). In this setting the loga- 
rithmic negativity is no longer limited by the average en- 
ergy per oscillator but instead it grows linearly with the 
system size n. The conclusion is that one part of the neg- 
ativity can be related to the energy; the second part can 
be seen as a surface term, proportional to the area of the 
boundary forming the contact between both bipartitions. 
In one spatial dimension this "area" is the number of con- 
tacts between both parts, which in the periodic setting 
described above equals n. This interpretation finds fur- 
ther support in the result for the logarithmic negativity 
of a symmetric bisection in an open chain of oscillators, 
which is then roughly half the logarithmic negativity of 
the corres ponding chain wi t h per iodic closure (see dis- 
cussion in lAudenaert et al l l2002h . 

An interesting puzzle is given by the analysis of the 
negativity of two connected parts of the chain that are 
separated from each other by a finite number of sites. Be- 
sides an expected all-over exponential decay of the log- 
arithmic negativity with the distance between the two 
equally large groups, the negativity is also of limited 
range (see top panel of Fig. [TT|l . This range increases 
with the size of the two parts. In particular is there no 
pairwise negativity between two single oscillators unless 
they are neighbors. This implies that the entanglement 
of distant groups of oscillators cannot be due to "free" 
pairwise entanglement of single oscillators (as opposed to 
"bound" entanglement not detected by the negativity). 
So either "bound" pairwise entanglement is responsible 
for the entanglement present in distant groups or multi- 
partite entanglement might play an important part. 

It is instructive to mention that pairwise correlations 
between single oscillators do exist notwithstanding a van- 
ishing pairwise negativity ( Audenaert et al\ . l2002h . In 
any case the presence of correlations is necessary for 
quantum entanglement of the parts. It is worth notic- 
ing that both the plateau exhibited by the negativity 
for not too small size of the parts, and the decrease of 
the single oscillators' negativity with the system size find 
a plausible explanation merely in terms of the correla- 
tion length. The same plausibility arguments predict 
the single oscillators' negativity to increase with grow- 
ing system size when open chains are considered; in fact 
this is what the autho rs observed (see the discussion in 
lAudenaert et aLl . l2002t ). The observed short range of the 
negativity in particular for small connected sets of os- 
cillators overstretches this simple reasoning and demon- 
strates that the connection between entanglement and 
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FIG. 11 Logarithmic negativity for harmonic chains. Top 
panel: The logarithmic negativity of two connected groups of 
s oscillators as a function of the distance between both groups. 
It is seen that the range of the negativity is finite but grows 
with the group size s. For s — 1, its range is zero, mean- 
ing that no two single oscillators share pairwise entanglement 
unless they are neighbors. Bottom panel: Logarithmic neg- 
ativity as in the top panel of fig. [TUJ but for an open chain. 
The value of the negativity is roughly half as large as for the 
periodic chain since the number of connection points of the 
parts is halved. In contrast to the periodic chain, the nega- 
tivity of one small part with the rest of the ope n chain grows 
with system size [From ([Audenaert et all l2002h ] 



correlations is indeed more subtle. 



C. Systems in d > 1 and the validity of the area law 

The scaling of the entanglement entropy in systems 
of higher dimensions has been subject of intense inves- 
tigation in various fields of research. In the context of 
quantum information the understanding of the scaling 
of the entropy as a function of the block size has im- 
portant consequences on the simulability of a quantum 
system by a classical algorithm and therefore is attract- 
ing a lot of interest. The picture that emerged from the 
analysis of the one-dimensional case, i.e. the violation 
of the area law when the system is critical, does not 
seem to hold in higher dimensions. The situations ap- 
pears more complex. The higher dimensional generaliza- 
tion of matrix product states, the projected entangled 
pair states, satisfy the area law ( Verstaete and Cirad . 
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12004 IVerstaete etaH l2006f) despite a divergent corre- 
lation length. Another example is ground-state of an- 
tiferromagneti c Ising-spin networks emb edded on planar 
cubic lattices (jWellard and Orusl |2004) where the area 
law is fullfilled also at the critical point. In the following 
of this section we concentrate on the ground state of some 
many-body Hamiltonian for which the block entropy has 
been recently computed. We first consider the case of 
hopping Hamiltonian of (free) fermions and bosons and 
then we discuss the case of harmonic lattices (note that, 
hysterically the bosonic case was considered first). 



1. Fermi systems 

In one dimension, by virtue of the Jordan- Wigner 
transformation, the block entropy of a system of interact- 
ing spins is tightly connected to that of a (free) Fermi gas. 
It is of great interest to understand what are the proper- 
ties of the block entropy for free fermio ns in d dimensions . 
This quest i on ha s been studied in ( Gioev and Klichl . 
120061: IWoli I2006T) where it was shown that logarithmic 
corrections persist also in higher dimensions 



(52) 



The expression of the constant of pro portionality in the 
equat ion above has been obtained by (IGioev and Klich 



I2006D resorting the Widom conjecture. IWo lf. 2006 



explo ited the quadratic lower bound of (jFannes et all 
2003). The corrections to the area law are a Fermi sur- 
face effect. In the case of fractal dimension of either 
the Fermi or the block surface, the scaling is modified 
into Sl T ~ l d ~P log 2 I where 1 — j3 is the maximum frac- 
tal enhancement of dimension of either the Fermi or the 
block boundary. An interesting case where the system 
u ndergoes a Lifh i tz ph ase transition has been considered 
in ICramer~e£ ~al\. 120071 Indeed as it was pointed out by 
Cramer et al these transitions related to a change in the 
Fermi surface manifest in a non-analytic behaviour of the 
prefactor of the leading order term entanglement entropy. 

For regular block and Fermi surface, numerical anal- 
ysis has confirmed t he modified a rea-law for crit- 
ical two-dimensional dBarthel et all l2006al: iLi et all 
|2006[ ) and three-dimension al ( Li et all l2006h mod- 
els. iBarthel et all l2006al study the tight binding 
model as an example for a two-dimensional model 
with a connected Fermi surface as well as the model 

H = -Ez.yK 1 + {- l ) V )<i,yC x ,y + l + 4,y C x+l,3/+l + 

f -1 ™fith a disconnected Fermi surface and 



+ (1 + (-l)*+»)ct c 



^x,y^x 

H — ~ Z->x,y l' b ^x,y^x+l,y T l iT V~ V " J^x,y^x,y+1} 

with a zero dimensional Fermi surface, as for the one- 
dimensional case. Whereas in the first two cases, the 
entropy is found to obey the modified area law, this is 
no longer true for the third model with zero-dimensional 
Fermi surface. There, the corrections to the standard 
area law S ~ l d ~ x are sub-logarithmic. 

The same feature has been observed by ( Li et aZ.l . l2Q06h 
studying the spin-less fermionic in two and three spatial 



dimensions. The authors conjecture an interesting con- 
nection between the modified area law to be observed 
and the density of states at the Fermi energy. They for- 
mulate this in terms of the co-dimension at the Fermi 
energy: i.e. the dimension of momentum space minus 
the dimension of the degeneracy at the Fermi energy in 
momentum space, providing a measure of the relative 
portion of the gapped excitations in the low lying spec- 
trum of the model. In agreement with the findings for 
the model zero-dimensional Fermi surface, the authors 
observe only sub-logarithmic corrections to the area law 
if the co-dimension at the Fermi energy is 2. The authors 
conjectured from this that in two spatial dimensions a co- 
dimension less or equal to 1 is necessary for the modified 
area law S ~ £ log 2 ( to apply. They do not mention im- 
plications of fractal co-dimension due to a fractal Fermi 
or block surface. This finding would be worth further 
investigation in direction to higher dimensions in order 
to fix the connection between area law and co-dimension 
at the degeneracy point. 

We finally mention the interesting connection between 
the block entropy and the Berr y phase in lattice models o f 
fermions recently discussed in (|Rvu and HatsugaiL l2006h . 



2. Harmonic systems 

Harmonic systems have been also investigated to un- 
derstand the validity of the area law. They provide one 
of the few physical systems fo r which exact anal ytical 
treatements are available (see ICramer et al. I, 120061 and 
references therein). 

We first consider a system of non-critical harmonic os- 
cillators with nearest neighbor interaction and periodic 
boundary conditions. Non-criticality implies that the 
lowest eigenvalue of the interaction matrix U, A m i„(U), 
is well separated from zero. Further peculiar characteris- 
tics of the covariance matrix, in particular its symmetric 
and circulant form, allowed to give estimates for upper 
and lower bounds of the block entropy of some compact 
rf-dimensional hyper-cubic re gion with edge len gth £ and 
surface proportional to i d ~ x ( Plenio et all [20051) . 

The upper bound has been established direct l y from 
the logarithmic negativity (see lAudenaert et al. ,'2002). 
whereas for the lower bound several estimates for the 
dominant eigenvalue of the reduced density matrix have 
been employed. 

Key ingredients to the problem are the largest eigen- 
value of the covariance matrix and the uncertainty re- 
lation which constrains all eigenvalues of the covariance 
matrix to lie above H/2. The result is that both bounds 
go proportional to hence the entanglement entropy 
is indeed proportional to the surface of the block. An 
extension to general block s hapes has bee n form ulated 
as well for Gaussian states (ICramer et all [2006) . For 



pure states, lower and upper bounds are given as be- 
fore, which both scale linearly with the surface area of 
the block. For a given finite range of interaction beyond 
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nearest neighbors, the area law could be stated on ly in 
terms of an upper bound (as in ( Plenio et aZ.I . l2005t ) . this 
upper bound has b een obtained from th e logarithmic neg- 
ativity). Therefore lCramer et ]iE|2006 could not exclude 
the block entropy to scale with lower dimensionality than 
the block-area. The area-law can even be present in dis- 
ordered systems; the crucial requirement for the area law 



to hold in this case is that K := 



is bounded for 



A„„„(U) 

all admissible disorder realizations. The presence of a 
finite temperature T enters only in the proportionality 
factor, when thermal equilibrium states (which are still 
Gaussian) are considered. 

It is interesting to notice that an algebraically diverg- 
ing correlation length d oes not automaticall y imply a vi- 
olation of the area law lCramer et al. 1. 120061 find that for 
non-critical systems of harmonic oscillators and algebraic 
decay r~ a of the correlations with the distance r still 
leads to the area-law as long as a > 2d. Essential for 
this conclusion is that k is bounded. The authors conjec- 
tured d as a tighter bound at least for hyper-cubic blocks 
(When compared to the one-dimensional Ising model, al- 
ready at a = d = 1 logarithmic corrections appear) . 

For critical systems the situation is different: critical- 
ity leads to eigenvalu es of U arbitrarily close to zero an d 
hence unbound k. In (jUnanvan and Fleischhauerl I2005T) . 
a one-to-one connection between criticality and logarith- 
mic corrections to the area law has been reported for 
one-dimensional systems with finite range interaction. 
The authors also report on evidence for this connec- 
tion to hold also in higher dimensional setups. How- 
ever the absence of logarithmic corrections for critical 
two-dimensional arrays of harmonic oscillators and near- 
est neighbor interaction h as been evidenced numerically 

. The analytical calculations 
conclude the area law to ap- 



in (Barthcl et al. 



of (jCramer et 




ply for arbitrar y number of dimensions, opposed to the 
conslusions of ( Unanvan and Fleischhauerl 20051) . As 
support for their claim. lUnanvan and Fleischhauerl . 120051 
quote a factoring interaction matrix, which however cor- 
responds to a noninteracting array of one-dimensional 
harmonic chains; this can be seen as the limiting case 
of an anisotropic interaction with finite range and does 
not give support to the original claim. It is indeed fun- 
damentally different from harmonic <i-dimensional lat- 
tice s with isotrop i c finite range interaction, as discussed 
in (jBarthel et all l2006al ICramer et all . l2007f ). Never- 
theless it raises the question for a critical anisotropy 
for the coupling of the harmonic oscillators, which on 
the background of t he findings in (jBarthel et all l2006al: 
ICramer et all 120071 ) could be phrased as: "Does a finite 
critical anisotropy exist beyond which the harmonic lat- 
tice is quasi one-dimensional?" 

Logarithmic corrections are also being observed, when 
infinite range interactions ar e considered, which drive 
the system towards criticality ( Unanvan et all [20071 ). To 
this end, the authors consider a two-dimensional array of 
harmonic oscillators, with an interaction of finite range in 
x-direction and an infinite range in the y-direction. This 




FIG. 12 Entanglement entropy for TV = 500 and £ = 125 as a 
function of the parameters h z and 7 of the ferromagnetic LMG 
model. For 7 / 1, there is an anomaly at the critical point 
h z — 1, whereas the entropy goes to zero at large h z since the 
ground state is a fully polarized state in the field direction. 
In the zero field limit, the entropy saturates when the size of 
the system increases. For 7 = 1, the entropy increases with 
the size of the system i n the region < h z < 1 and jumps to 
zero at h z = 1. [From l|Latorre et all [2005bl )] 



is a very instructive example in that it leads to the loga- 
rithmic correction ~ l x lnl y . Interestingly, the logarithm 
contains the length l y of the block, where the interaction 
has infinite range. The prefactor of the logarithm is half 
the length of the block in x direction (with finite range 
interaction). 



D. LMG Model 

The logarithmic divergence with the block size of the 
entanglement entropy is not exclusive of one-dimensional 
systems. The block entropy of the LMG model was 
studied both i n the ferromagnetic fearthel et all 2006b; 



Latorre et ~al\, l2005bl : [Stockton et al. 



20031: [Vidal eTal. 
20071 ) and antiferromagnetic (jUnanvan et aZ.I . I2005T ) case. 
In the LMG model each spin is interacting with all the 
other spin in the network therefore the idea of a block 
as depicted in FigJS] does not fit very well. Nevertheless 
it is perfectly legitimate to define the reduced entropy 
of £ spins once the other N — £ (N is the total number 
of spins) have been traced out. Evidently the entropy 
is independent on which spins have been selected to be 
part of the block. In Fig[T2]the representative behaviour 
of the entanglement entropy as a function of the various 
regions of the phase diagram is shown for a fixed value 
of the block and system sizes. Below we summarize the 
main findings related to the LMG model. In the ferro- 
magnetic case and in the case of 7 = 1 (isotropic model) , 
the entropy diverges logarithmically with the block size 
S ~ (1/2) log£ while at fixed £ and N diverges when the 
external magnetic field app roaches 1 from below. Also in 
the antiferromagnetic case ( Unanvan et all 120051 ). when 
no transition as a function of the field is present, the en- 
tropy grows logarithmically with the size of the block in 
the isotropic limit. Differently from the one-dimensional 
case where the prefactor is universal and related to the 
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FIG. 13 Entropy S between two-level atom and cavity 
field as a function of the atom-radiation coupling for both 
N —* oo and finite N. In the inset the scaling of the 
value o f the entanglement m aximum as a function of logiV. 
[From (|Lambert et all\2004 )] 



central charge, here the origin of the prefactor of the log- 
arithmic divergence is related to the presence of the Gol- 
stone modes and to the number of vanishing gaps. The 
recent work of Vidal et al clarifies this issue by study- 
ing a number of collective spin model by means o f 1/JV- 
expansion and scaling analysis ( Vidal et~al\ , 120071) . 



E. Spin-boson systems 

The entropy in models of spins interacting with 
harmonic oscillators have been analyzed as well. 
Here the separation between spin and bosonic de- 
grees of freedom is natural and the partition leads 
to study the reduced entropy of one subsystem 
(say the spin). The entangle ment entropy wa s 
studied for the Jaynes-C ummings dBose et all l200l[) 



Tavis-Cummings model ([Lambert et al. 
and for the spin-boson model (ICosti and McKenzie 



12004 12005 1 



20031: iJordan and Buttikerl . 12004 iKopp et all 12006 : 



Stauber and Guineal . l2006af) . Lambert and coworkers an- 
alyzed how the super-radiant quantum phase transition 
manifests in the entanglement between the atomic en- 
semble and the field mode. They compute the von Neu- 
mann entropy numerically at finite N and analytically in 
the thermodynamical limit. They found that the entropy 
diverges at the phase transition as (see Fig. 13) 



S~-(l/4)log|A-A c | 



(53) 



where A is the coupling between the spins and the boson 
field and A c is the value at which the transition takes 
place (see Fig. 13). The entropy of the spin-boson model 
(see Eg. (1321) ) was studiedby numer ical renormalization 
group in (jCosti and McKenzid . 120031 ). An analytic calcu- 
lation, inc luding other dissipat i ye mo dels, has been pre- 
sen ted in (IStauber and Guineal. 12004!) and more re cently 
in (|Kopp et aU^OOetlStauber and Guineal . l2006al lbl). In 
the broken-symmetry state has an effective classical de- 
scription and the corresponding von Neumann entropy 



is zero. In the symmetric phase the entropy can be 
easily expressed as a function of the ground state en- 
ergy, and 8 (defined in Eq. (|52")) ). The previous ex- 
pression allows At the transition p oint the entropy is 
discontinuous with a jump given by (jKopp et all 120061 ) 
AS = In 2 + 8/4S c \n(S/S c ) (S c is a high energy cutoff). 
A systematic analisys of the entropy in the spin-boson 
mo del for different coupling regimes was pursued re cently 
in (|K.Le Hur et all l2007t IKopp and Le Hurl . 120071) . 

We finally mention the interesting connection be- 
twe en entanglement and energy fluctuations introduced 
by ( Jordan and Buttikerl 120041 ) and exploited in details 
both for a spin and for an harmonic oscillators cou- 
pled to a bath. This connection might be useful in the 
light of p ossible experiment al measure of entanglement 



(see also iKlich et all 120061). For example, as pointed 



by ( Jordan and Buttikerl . 120041 ) . in certain mesoscopic 
realization of qubits as metallic rings of superconduct- 
ing nano circuits, a measurement of persistent current 
can be directly related to a measurement of the entropy. 



F. Local entropy in Hubbard-type models 

A very important class of interacting fermion models is 
that of Hubbard type models (sec lIII.Bj) . First studies of 
ent anglement in th e one-di mensional case h ave appeared 
in (jKorepinl 120041) and in (|Gu et aU 120041) . Most of the 
studies in this type os systems analyzed the properties of 
the local entropy. 

Gu et al analyzed the local entropy for the one- 
dimensional extended Hubbard model for fermions with 
spin 1/2. Due to the conservation of particle number and 
^-projection of the spin, the local density matrix of the 
system takes the simple form 



= z|0) <o| + ti+|T><TI + «"U>UI + <»IU><ni 



independent of the site number j because of translational 
symmetry. The broken translational invariance in the 
charge density wave phase has not been taken into ac- 
count in this work. This does not affect the central re- 
sult but might affect the entropy within the charge den- 
sity wave phase. Except the superconducting phase, the 
phase diagram at half filling (for /i = 0) of this model 
has been nicely reproduced by the contour plot of the 
local entropy (see top panel of Fig. [14]), where the phase 
transition coincides with its crest. This turns out to be a 
general feature of local entropies - also for spin models - as 
opposed to entanglement class specific measures, as e.g. 
the concurrence for pairwise entanglement whose maxima 
in general appear at a certain distance to quantum crit- 
ical points and hence are not associated to the quantum 
phase transition. In view of the monogamy of entangle- 
ment this is interpreted as evidence of dominant multi- 
partite entanglement in the vicinity of quantum phase 
transitions. 

For the Hubbard model Eq. (|3"0")) and fixed U/t, the local 
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FIG. 14 Upper panel: The local entropy as a function of 
the on-site coupling u and nearest neighbor coupling v. The 
contour plot below includes the known phase diagram of the 
model (full black line). Except the superconducting phase, 
the phase diagram is nicely reproduced. Lower panel: The 
local entropy for the standard Hubbard model and different 
on-site couplings u as a function of the filling (1 corresponds to 
half filling). Except for U/t = 0, where a maximum appears, 
this quantity shows a dip at half fill i ng, i. e. at the metal- 
insulator transition. [From (|Gu et all [2004)] 



entropy as a function of n shows a dip at the critical fill- 
ing fraction n c = 1, where a metal-insulator transition 
occurs (for U > 0) (see bottom panel of Fig. [T4")l . For 
the two limiting cases U = 0,oo the maximum instead 
is located at filling fractions, where the ground state is a 
singlet of the largest symmetry group. Gu et al conjec- 
ture that this was true for general U > and then the 
presence of an unknown phase transition at these max- 
ima. 

This analysis clearly points out that the local entropy 
indicates different phase transitions in different ways, es- 
sentially depending on whether the quantity is sensitive 
to its order parameter or not. Due to the u(l) sym- 
metry of the model, the single site reduced density ma- 
trix is a functional of occupation numbers only. These 
operators cannot, however, describe order parameters of 
superconductivity or some order parameter of the metal- 
insulator transition. Indeed, the superconducting phase 
can be predicted if the entropy of entanglement is cal- 
cula ted for a block of spins, instead of for just a single 
site ( Deng et~al\ , l2006h . A reduced density matrix of at 
least two sites is necessary for being sensitive to supercon- 



ducting correlations [see also Ref. iLegeza and Solvoml . 
2006 for a similar result obtained for the ionic Hubbard 
model] to be seen. 

Another model studied with this respect is the so called 
bond-charge extended Hubbard model (see IIII.B|) . In 
phases II and III of fig. 15 there are superconducting 
correlations which are due to 77-pairing and hence indi- 
cate the prese nce of multipa rtite entanglement, as dis- 
cussed before ( Vedrall . [2004al) . At the bond-charge cou- 
pling corresponding to x = 1, th e entanglemen t of the 
model has been analyzed in Ref. (lAnfossi et all [2 005b) 
and for gen eral x and n = 1 in Ref. ( Anfossi et all 
l2005al 120061 ). Besides the local e ntropy of entangle- 
ment |_SVj_Jn_^^foss^i^L j_p005t]) , also the negativ- 
ity ( Vidal and WernerT" 20021) and th e quantum mu- 
tual information ( Grois man et ~al\. 120051) have been used. 
While Si measures all (pairwise and multipartite) quan- 
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FIG. 15 Upper panel: The ground state phase diagram of 
the Hirsch model at x = 1. Empty, slashed and full cir- 
cles indicate the presence of empty, singly and doubly occu- 
pied sites, respectively. Bottom panel: Except the insulat- 
ing lin e IV, the phas e dia gram is nicely reproduced by d u Si. 
[From (|Anfossi et all l2005al )] 

turn correlations involving this specific site, the negativ- 
ity offer a lower bound for the quantum correlation of 
two specific sites, and the mutual information accounts 
for pairwise quantum and classical correlations. There- 
fore, this combination of correlation measures opens the 
possibility to decide, what type of correlation is relevant 
at a quantum phase transition. The results are shown 
in the upper panel Fig. 1151 The different phases are 
shown of Fig ll5l they are discriminated by local occupa- 
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tion numbers as described in the top panel; consequently, 
the entropy Si bears the information on all the phase di- 
agram except the insulating line IV. This is seen from the 
plot of d u Si (with u = U/t) as a function of the on-site 
Coulomb coupling u and the filling fraction n. A com- 
parison of first derivatives respect to y — n, u (depending 
on the phase transition) of all three correlation measures 
reveals common singularities for d y Si and d y Xij only for 
the transitions II-III and II-IV; furthermore, for both of 
them it was proven that the range of the concurrence R 
diverges! Anfossi et al\ .\2006). These facts allow to char- 
acterize the transitions II-III and II-IV (at n = 1 and 
arbitrary x) as governed by pairwise entanglement, that 
is the more spread along the chain the closer the transi- 
tion is got. For the transitions II-I and II-I' instead the 
multipartite entanglement is relevant, with a finite range 
of the concurrence. A similar behaviour was encountered 
for non-critical spin models where the divergence of R 
is accompanied by the emergence of a fully factorizcd 
ground state (see Sec. IIIII and IIV.AT]) . Here R diverges 
close to QPT; it was also noticed by Anfossi et al. that, 
while the ground state in IV is indeed factorized, the 
phase III is made of a superposition of doubly occupaied 
and empty states. Observing that the pairwise entan- 
glement is vanishing in both the phases III and IV, the 
authors conjectured that the divergence of R marks an 
'entanglement transition' in solely in the pairwise entan- 
glement. 

In order to detect the transition II-IV at n = 1 
and any x, d x Sj has been calculated by means of 
DMRG (I Anfossi et all l2005al) . Its singularities allow to 
accurately determine the charge gap as a function of the 
bond-charge coupling x. 

We now proceed with the Hubbard model in a mag- 
netic field (see IIII.Bj) Also here, the local entropy Si has 
been looked at in order to analyze its entanglement. As in 
the examples before, Si indicates the second order phase 
transitions in terms of divergences of its derivatives dhSi 
and d^Si, respectively. Indeed, it has been demonstrated 
that dhSi and d^Si can be expressed in terms of spin an d 
charge susceptibilities (jLarsson and Johannesso nl l2005h . 
hence bridging explicitly between the standard method in 
condensed matter physics for studying phase transitions 
and the approach from quantum information theory. 

The local entropy for thebosonic version of the Hub- 
bar d model, the Bose-Hubbard model, was computed 
in dBuonsante and Vezzanll . l2007t iGiorda and ZanardU 
I2004E " for different graph topologies. The attention in this 
work was on the dependence of the entanglement on the 
hopping amplitude. The authors showed that for certain 
classes of graphs the local entropy is a non-monotonic 
function of the hopping. Also for the bosonic case the 
local entropy is a good indicator for the presence of a 
(superfluid-insulator) quantum phase transition. 

Summarizing, the body of work developed sug- 
gest to conclude that that local entropies can detect 
QPTs in systems of itinerant fermions, particularly 
if the transition itself is well predicted by a mean 



field approach for local observables o f the model, [see 
also ( Larsson and Johannessonl l2006t )] . Furthermore, 
translational invariance is necessary for the prediction 
to being independent of the site, the local entropy is cal- 
culated for. If this symmetry is absent, it might prove 
useful to average over the sites; the resulting measure is 
then equivalent to the Q- measure ( Meyer and Wallachl . 
I2002D . 



G. Topological entanglement entropy 

We close this section by summarizing the ongoing 
research activity studying the subleading corrections 
to the block entropy in the two-dimensional systems. 
Most of the results were demonstrated for quantum two- 
dimensional lattices (though generalizations two higher 
dimension is straightforward). 

Fradkin and Moore considered quantum critical points 
in two spatial dimensions with scale invariant ground 
state wave functions, characterizing for example the scal- 
ing limit of quantu m eight vertex models and non abelian 
gauge theories (see lArdonne et al. I. b M ill for a recent ref- 
erence). The main result is that a universal logarith- 
mically divergent correction, determined by the geome- 
try of the partition of the system, emerges in addition 
to the area's law term in the e ntanglement block en- 
tropy (jFradkin and Moorell2006l ) . 

Such work benefit of the earlier seminal con- 
tribution of Kitaev and Preskill, and Levin and 



Wen (|Kitaev and Preskilll l2006t iLevin and Wenl l2006f) 
demonstrating that the correction to the area law is of 
topological origin. Namely the entanglement entropy was 
demonstrated to scale as 



s = ae-^ + oir 1 ) 



(54) 



The coefficient a is non universal and ultraviolet diver- 
gent. In contrast, the quantity 7 can be expressed as 
log -D, where D is known as the total quantum dimen- 
sion^ is universal, and related to the topological winding 
number of the theory. The calculations were pursued 
with methods of topological field theory, giving an ex- 
plicit expression for 7. 

We remark that the result of Kitaev and Preskill, and 
Levin and Wen could provide an a l ternat ive avenue to 
lattice gauge theory methods! Wenl . 120041 ). detecting a 
genuine topological order in the system (when descrip- 
tions based on local order parameters fails) by direct in- 
spection of the wave function. 

When a topological order is present, the ground state 
of the system acquires a peculiar degeneracy when its lat- 
tice has a non trivial genus g. Hamma et al. studied the 
entropy of the Kitaev model on a two dime nsional lattice 
with generic genus g. The Kitaev model ( Kitaevl 120031 ) 
is a two dimensional exactly solvable lattice model with 
double periodic boundary conditions, whose Hamiltonian 
can be realized with a set of spins in a square lattice with 
ring exchange and vertex interactions. The ground state 
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of the Kitaev model is characterized by the presence of 
a topological order 6 . For such systems Hamma et al. 
related the degenerac y of the ground state (that is 4 s ), 
to the block entropy (jHamma et al. I. l2005alfbh . 

As a step toward models with more generic topologi- 
cal orders, the topological entropy was studied numeri- 
cally for the quantum dimer model in a triangular lat- 
tice (jFurukawa and Miseuichl. 1200a) and for fractional 



quantum Hall states ( Hague et all [2006) 



We finally mention that the topological term in the 
entanglement entr opy in th e context of quantum gravity 
was ev idenced in ( Fursaevl 120061 : I Ryu and Takavanagil 
l2006bri . 



H. Entanglement along renormalization group flow 

One of the original motivations put forward by 
Preski ll to investigat e entanglement in many-body sys- 
tems (Pr eskillL l2000h was the idea that quantum infor- 
mation could elucidate some features of the renormal- 
ization group which is a cornerstone method in mod- 
ern physics. It is natural to think that the proce- 
dure of tracing out high energy modes in a renormal- 
ization group step has some kind of irreversibility built 
in. Quantum information concepts could prove to be 
useful in elucidating issues related to this irreversibility 
and possible help could come from relating the celebrated 
c— theorem ( Zamolodchikovl Il986l ) to the loss of informa- 
tion. 

Entanglement loss along a reno rmalization group tra - 
jectory was studied in spin chains (jLatorre et a/.l . l2~005aF ). 
More recently a number of relations relating renormal- 
ization group, conformal fiel d inva r iance and entangle- 
ment loss were derived in (|Orusl . I2005T ). According 
to ( Latorre et al .], l2005al) entanglement loss it can be 
characterized at three different levels: 

- Global entanglement loss - By using the block entropy 
as a measure of entanglement, for which we know the re- 
sult of Eq. (|4"2"|) , and an inequality on the central charges 
which derives from the c— theorem, it follows that S uv > 
S IR . The block entropy at the ultraviolet fixed point can- 
not be smaller than that at the infrared fixed point. 

- Monotonous entanglement loss - It is also possible to 
follow the entanglement along the whole transformation. 
Also in this case the entropy is a non-decreasing function 
along the flow. As a simple example one can consider the 
block entropy of and Ising chain close to the critical point 
which goes as S ~ In | A — 1 1 , from which monotonicity 
follows. 

- Fine-grained entanglement loss - The monotonicity of 
the entanglement seems to be present at a deeper level 



1 The Kitaev model was suggested to provide a realization of the 
so caaled 'toric code', namely a topological quantum computers 
made by a physical system with anionic excitations (see also 
Castagnoli and Rasctti, 1993) 



in the structure of the dens i ty matrix. I t is p ossible to 
show (jLatorre et all l2005ar . lOrus et all l2006h through 
majorization relations that the spectrum of the reduced 
density matrix gets more ordered along the flow. By 
denoting with r m the eigenvalues of the reduced den- 
sity matrix p, majorization relation between the two 
set of spectra (corresponding to two different parame- 
ters) means that there is a set of relations for which 



> 



for n 



d (d is the di- 



mension of p). 

Motivated by ultraviolet divergencies of the entropy 
of entanglement in quantum field theory, Casini and 
Huerta introduced a quantity F (A, B) related to the en- 
tropy measuring the degree of entanglement between the 
two regions A and B. The function F is defined as 
F(A, B) = S(A) + S(B)-S(AnB)-S(AUB) which co- 
incides with the mutual information Eq. §§§ in the case of 
non-intersecting regions. In two dimensions, it is a finite 
positive function with the property F(A, B) < F(A, C) 
for B C C if AnC = 0. Then, for sets with a single (path- 
connected) component in two dimensional conformal field 
theories they showed that it allows to prove an alternative 
entrop ic version of the c-theorem ( Casini and Huertal . 
l2005bri . 



VI. LOCALIZABLE ENTANGLEMENT 

A. Localizable entanglement and quantum criticality 

The study of localizable entanglement (see Sec- 
tion Ill.Cj) in spin chains allows to find a tighter con- 
nection between the scales over which entanglement and 
correlations decay (we saw in the previous sections that 
the two spin entanglement, expressed by the concur- 
rence, does not decay on a the same range of correla- 
tions) (jPopp et aLl . l2005l . 120061 1 Verstraete et a/.l . [2004ah . 
One expects that the procedure of entangling distant sites 
by a set of local measurements will be less effective as the 
distance between the two particles increases thus leading 
to a definition of entanglement length £e. For a transla- 
tional invariant system £g can be defined in analogy of 
the standard correlation lenght 



= - lim log 

\i-j\->c 



E loc (\i-j\) 



J\ 



(55) 



By definition the entanglement length cannot be smaller 
than the correlation length, £e > £, therefore at a sec- 
ond order phase transition the localizable entanglement 
length diverges. In addition there may also appear "tran- 
sition points" associated solely to a divergence in £g. In 
order to avoid misinterpretations, it must be stressed that 
the localizable "classical" two-point correlations then di- 
verge as well. Thus, the essence of the phenomenon 
is that correlations can be localized between arbitrar- 
ily distant sites by means of suitable local operations 
and classical communication despite a finite correlation 
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length; n ecessary for this is the presen ce of global entan- 
glement ( Popescu and RohrlichL Il992l ) . 

Fo r the Ising model in a tran sverse field it can be shown 
that (|Verstraete et all l2004b l ) 



max \Q%\ 

y.—x,y,z 



< Ei oc {i 



(56) 



(1 ± (SfS?)) 2 - ((SI) ± (Sp) 2 and Qt 



where s± 

(S?S?)- (S?)(S? , . 

In this case, the lower bound in Eq. (|56p is determined 
by the two-point correlation function in the x-direction. 
In the disordered phase (A < 1) the ground state pos- 
sesses a small degree of entanglement and consequently 
its entanglement length is finite. The situation changes 
at the other side of the critical point. Here, although 
the correlation length is finite, the entanglement length 
is infinite as asymptotically the correlation tends to a 
finite values. The divergence of £g indicates that the 
ground state is a globally entangled state, supporting the 
general idea that multiparti te entanglement is most rel- 
evant at the critical point ( Osborne and Nielsen! . l2002t 
iRoscilde et al. 1. 120041) . 

The properties of localizable entanglement were 
furt her investigated for a spin- 1/2 XX Z - chain 
in ( Jin and Korepinl l2004al : IPopp et al. , 2005) as a 
function of the anisotropy parameter A and of an 
externally applied magnetic field h. The authors used 
exact results for correlation functions relying on the 
integrability of the models to find the bounds in Eq. (|5"6"|) . 
For the antiferromagnetic XXX-case they provided the 



following lower bound Ei oc (i — j) > 



In H—j\ 



The 



presence of the anisotropy increases the lower bound 
of the localizable entanglement. At the Berezinskii- 
Kosterlitz-Thouless critical point (A = 1) the lower 
bound of the n earest-neighbour lo calizable entanglement 
shows a kink I Popp etal\ . l2005h . As pointed out by 
the authors this might have implications in the general 
understanding of the Berezinskii-Kosterlitz-Thouless 
phase transitions where the ground state energy and its 
derivatives are continuous as well as the concurrence 
(seellVAJJand Fig©. 

The localizable entanglement in two-dim e nsiona l 
XX Z model was discussed as well ( Svliuasenl . [2003b) 
by means of quantum Monte Carlo simulations. A lower 
bound has been determined by studying the maximum 
correlation function which for A > — 1 is Q x , the long- 
range (power law) decay of the correlation implying a 
long ranged localizable entanglement. 

The definition of localizable entanglement has an inter- 
esting connecti on with the concept of quantum repeaters 
introduced in ( Briegel et all [l998h . Quantum repeaters 
have been designed to enhance the transmission of entan- 
glement through noisy channels. The idea is to distribute 
along the channels a number of intermediate sites where 
a certain number of local operations are allowed in order 
to maximize the entanglement between the transmitter 
and a receiver. This is the very definition of localizable 



entanglement. 

Localizable entanglement has been defined in Scc lII.CI 
as an average over all possible measuring processes it 
is of interest to understand also the statist ical fluctua- 
tions around this average value. To this end ( Popp et all 



2005) analyzed also the variance associated to the en 
tanglement fluctuations: 5Ef oc — J2 s Ps^ 2 (\' l l } s)) — Ef oc , 
where E is a measure of pairwise entanglement. The 
fluctuations of the entanglement increase in the vicinity 
of a critical point. This was checked explicitly for the 
one-dimensional Ising model. 

As detailed below, additional interesting results were 
obtained for spin-1 systems where a true transition in the 
entanglement (with a diverging £e but finite correlation 
length) has been found. 



B. Localizable entanglement in valence bond ground states 

For half-integer spins, gapped non degenerate ground 
states are characteristic for systems in a disordered phase 
(consider paramagnets for example). A finite gap in 
the excitation spectrum of the system in the thermo- 
dynamic limit makes the correlations decaying exponen- 
tially. This is the Lieb-Schultz-Mattis theorem, estab- 
lishing that, under general hypothesis, the ground state 
of a spin syste m is either unique and gapless or gapp ed 
and degenerate (jLieb et "aZl. ll96ll ) fsee lHastingsll2004l for 
recent results). It was a surprise, when Haldane discov- 
ered t hat systems of in teger spins can violate this the- 
orem (jHaldand . Il983al lbl). This suggests to investigate 
whether the entanglement in the ground state might play 
some role in establishing the hidden order characteristic 
for the Haldane phases. An aspect that might be relevant 
to this aim was recently addressed by s tudying the local- 
izable entanglement in AKLT models ( Verstraete et all 
2004a). The ground state of this class of models is of the 
valence bond type as discussed in IIII. A.4~1 



Iff*} 



AKLT 



= (®kA k i) |/)oi|/)l 2 ...|I)jV 



NN+1 



(57) 



\I) are singlets and A are 3x4 operators projecting the 
Hilbert space of two combined spins on its symmetric 
part, at the given site. This is nothing else than a matrix 
product state (see Section HlI.A.4|) . For this state it was 
demonstrated that a singlet state made of two spins-1/2 
located at the ends of the chain can be always realized 
(see FigUJ. This implies that the localizable entangle- 
ment is lon g ranged despite the exp onentially decaying 
correlation (jVerstraete et al 1. l2004aD . 

The localizable entanglement can be related to the 
string order parameter 0" tring defined in Eq. (f2"8")l . 
The key to make explicit this relation is to observe 
that the localizable en tanglement can be calculated 
as an expectation value (jCampos Venuti and Roncagliai 
120051) as L(\ip)) = (ip\a G s (ip)a N \tp) where G s = 
S s \s)(s\sign((vp\<j y ®<T y \ip}), and |s) is the optimal basis 
which maximize the entanglement of assistance. For the 
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AKLT model, the expression above reads 



N 



L{\9s) 



AKLT 



)=<n< 



(58) 



In this case, both the localizable entanglement and all the 
three components of the string order parameter saturate. 
Perturbing the AKLT ground state, namely making the 
resonating valence bonds with non-maximally entangled 
states (| 10) - e"^ | 01))/ a/2, the relation between the 
hidden order and the localizable entanglement is weakned 
(as compared to the AKLT model): it can be demon- 
strated that the string order parameters for a — x, z are 
finite, while the localizable entan glement vanishes expo- 
nentially with the deformation </> ( Popp e£~aZl . l2005l ). but 
a tight connection of the localizable entanglement and 
0^™„ is observed to persist also for the deformed 



string 



ground state ( Campos Venuti and Roncaglial [2005h . 

The valence-bond-solid phase order was further stud- 
ied by looking at the hidden order in chains with more 
complicated topology. The von Neumann entropy was 
studied in spin-1 XXZ model with b iquadratic inter- 
action and single ion a nisotropv in (iGu et all 2006t 
IWang et a/.l . l2005l ) and in (|Campos Venuti et al.l . l2006al) . 
Some of the features of the corresponding phase dia- 
gram are captured. The Haldane transitions exhibited 
in the phase diagrams are marked by anomalies in the 
Von Neumann entropy; its maximum at the isotropic 
point is not related to any critical phenomenon (the sys- 
tem is gapped around such a point), but it is due to 
the equi-proba bility of the three spin-1 stat es occurring 
at that point ( Campos Venuti et all l2006al ). Since the 
Berezinskii-Kosterlitz-Thouless transition separating the 
XY from the Haldane or large- D phases connects a gap- 
less with a gapped regime, it was speculated that an 
anoma ly in the entang lement should highlight such tran- 
sition (|Gu et aUl200l . 



VII. THERMAL ENTANGLEMENT 

Though the very nature of entanglement is purely 
quantum mechanical, we saw that it can persist for 
macroscopic systems and will survive even in the ther- 
modynamical limit. Entanglement survives also at fi- 
nite temperatures. This temperature could be as high as 
100 Kelvin in hi gh-temperature super conductors (jVedrall . 
l2004al) (see also lFan and Llovdl . l2005l ). In this section we 
review the properties of entanglemen t in many- body sys- 
tems at finite temperatures (see also lAnders and Vedral . 
2007). We will see that the analysis could shed new 
light on the interplay between the quantum nature of the 
system and its thermodynamics. Moreover, somewhat 
surprisingly, macroscopic state variables can be used to 
detect entanglement. Thermodynamics describes large 
scale systems with macroscopic properties, its state vari- 
ables T, N, V, p, its external fields h and its response func- 
tions the susceptibility and heat capacity, respectively 



X, C, .... Addressing entanglement as a thermodynamical 
property raised significant amount of interest in various 
communities. One wants to know, for example, under 
what conditions can we detect and extract entanglement. 
Can we see entanglement itself as a state variable, just 
like pressure is for a collection of atoms in a gas? What 
could be the corresponding thermodynamical potential? 
Is entanglement extensive? Since entanglement is closely 
related to entropy, we would expect the answer to the 
last question to be "yes" . 

The states describing a system in thermal equilibrium 
states, are determined by the Hamiltonian and the in- 
verse temperature /3 = (1/T). The density matrix is 

p = Z~ 1 e~ ,3H . where Z — tr[e~P H ] is the partition func- 
tion of the system. The thermal states expanded in the 
energy eigenbasis \ei),i = 0, 1, ... are then 



-0E O 



-pEi 



|eo)(e | + 



-0E 1 



l e i)( e i 



(59) 



Any separable state, or classical state, with respect 
to this split into subsystems A,B,C,D,.. (for exam- 
ple the sites of a spin system) can then be written as 
convex mixture of tensor products of states of the re- 
spective subsystems A, B,C, D, ... with probabilities pi 



'Pi 



If the state in Eq. ljSl?]) 



cannot be written in the form given above then it is en- 
tangled. In this section we will discuss the properties of 
this thermal entanglement. 



A. Thermal pairwise entanglement 

Extensive efforts have been made to understand 
how to quantify thermal entanglement in many-body 



systems starting from t he init i al paper s ( Arnesen et 



2001) 



Sev- 



120011 : iGunlvcke etHll . 120011 : iNielsenl . 
cral models of interacting spins in arrays were dis- 
cussed. Entanglement as measured by concurrence 
was shown to exist at nonzero temperatures in the 

oisen- 



was snown to e xist at nonzero te mperatures m 
trans yerse Ising (lOsborne and Nielsen!. 20021). He 



Wand. 12 002b 



120061 : 

XXZ (ICanosa a nd Rossignoli, 



berg ( Asoudeh and Karimipour . 2004 : Tribedi and Bosel . 



Wang and Zanardil l2002h . 
20051 12001 . XYZ 
models (Rigolinl. l200l" IZhang and Zhul l2006h . ferrimag- 



netic t ransition (Wang and Wand . 2006 ) and spin-one 
chains ( Zhang and Li . 2005h . Sevaral non-trivial aspects 
of the behaviour of the pairwise entanglement at finite 
temperatures can be illustrated by considering the 
simple case of two-sites systems. 



a. XXX Model We start our discussion on thermal en- 
tanglement by considering the XXX antiferromagnet 
(see Eg . (|2"51 ) . In this case the thermal state of this sys- 
tem can be written as a Boltzmann mixture of the singlet 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 

h z 

FIG. 16 Entanglement between two qubits interacting ac- 
cording to the antiferromagnetic Heisenberg model as a func- 
tion of the external field h z and temperature (multiplied by 
the Boltzmann's constant) T with coupling J = 1. The top- 
most plot shows the variation of nearest-neighbor entangle- 
ment with the magnetic field. The middle and the bottom- 
most plot show the same for next n earest and next to ne xt 
nearest neighbors, respectively. From lArnesen et aUl200ll 



and the triplet states: 



p T = -e 3f}J \^-)(ip-\+e-^ J+2h *>\00)(00\ 

+ e- /3 ( J - 2 ' i >|ll)(ll| + e- /3J |^+)(V; + | (60) 



where Z = e 3f3J + e -^ J+2h '> + e -^ J ~ 2 ^ + e~P J is the 
partition function of the system and ) = |01) ± |10). 
For the sake of simplicity we focus our attention only on 
two regimes 

The first regime is when the coupling J is large com- 
pared to the external field h. The ground state is then 
the singlet state, and at low temperature the system is 
therefore entangled. At higher temperatures the triplet 
becomes mixed into the singlet, and when (roughly) 
T > J/k, the entanglement completely disappears (when 
the external field is zero). Therefore, in order to have 
high-temperature entanglement in dimers we need a large 
value of the coupling constant J. 

When J is fixed, the second regime occurs when we can 
vary the value of the external field h. When h is large 
(greater than 2 J), the ground state is |11) and at zero 
temperature the dimers are therefore not entangled. The 
point where the singlet state becomes replaced by 1 11) as 
the ground state corresponds to the quantum phase tran- 
sition (occurring in the thermodynamical limit). How- 
ever, if we start to increase the temperature, the sin- 
glet state — which is the first excited state under these 
circumstances — starts to become populated. Entangle- 
ment can then be generated by increasing the tempera- 
ture. The behaviour of entanglement as a function of the 
magnetic field is shown in FigrTBI 



b. Ising model Another interesting case is that of a Ising 
coupling. At zero temperature only the lowest energy 
level is populated. In the case N = 2 the tangle of this 
pure state can easily be calculated from the density ma- 
trix, for h > 0, n = J 2 /{J 2 + hi). It is clear that 
the entanglement is highest for nearly vanishing magnetic 
fields and decreases with increasing field amplitude (this 
expression however is not valid for strictly h = 0, where 
no entanglement is present). Let us now turn our atten- 
tion to the case of non-zero temperatures. For a general 
pure state only one of the eigenvalues of the Hamilto- 
nian weight Eq. (|59p is non-zero and therefore equal to 
the tangle. For low temperature and magnetic field, i.e. 
h, T <C J, it is a good approximation to assume that only 
the two lowest energy levels are populated. In this case, 
the combination of the two lowest states also combines 



their concurrences in the following way (jGunlvcke et all , 
120011) : 

C = max{\w C -w 1 C 1 \,0}, (61) 

where the index refers to the ground state, while 1 
refers to the excited state and wq and u>i are the thermal 
weights of the ground and excited states respectively, see 
Eq. (|59[) . This phenomenon has been named as concur- 
rence mixing. In this case, the first excited state is the 
Bell state, |*") = (|01) - |10))/\/2, and Eq. (JBTJ reduces 



to C 



\wo- 



— Wi\. In general, the first term in 



the above equation is larger than the second, and in this 
case the concurrence decreases with temperature as wo 
decreases and wi increases. Moreover it is also possible 
to see that, for a given temperature, the entanglement 
can be increased by adjusting the magnetic field and is 
generally largest for some intermediate value of the mag- 
netic field. This effect can be understood by noting that 
wo increases with increasing h as the energy separation 
between the levels increase, but J / y/ J 2 + h 2 decreases. 
As a result the combined function reaches a peak as we 
vary h and decreases subsequently, inducing analogous 
behavior for the concurrence. 



B. Pairwise entanglement in the T critical region 

At finite temperatures but close to a quantum criti- 
cal points, quantum fluctuatio ns are essential to describe 
the properties of the systems ( Chakravartv et all , Il989t 
ISachdevl . [1999). In order to study the interplay between 
the thermal entanglement and the quantum fluctuations 
caused by the critical point at T = the analysis on small 
clusters is no longer sufficient. For the presentation we 
consider the one-dimensional quantum XY models. Such 
systems cannot exhibits any phase transitions at finite 
temperature, but the very existence of the quantum crit- 
ical point is reflected in the crossover behaviour at T ^ 0. 
The renormalized-classical crosses-over the quantum dis- 
ordered regimes throu gh the so called quantum critical 
region (|SachdevL fl999h . In the T — h plane a F-shaped 
phase diagram emerges, characterized by the cross-over 
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temperature customarily defined as T cross = |A _1 — 
For T <C T cross the thermal De Broglie length is much 
smaller than the average spacing of the excitations; there- 
fore the correlation functions factorize in two contribu- 
tions coming from quantum and thermal fluctuations sep- 
arately. The quantum critical region is characterized by 
T T cross . Here we are in the first regime and the cor- 
relation functions do not factorize. In this regime the 
interplay between quantum and thermal effects is the 
dominant phenomenon affecting the physical behaviour 
of the system. 

Thermal entanglement close to the critical point 
of the quantum XY mod els was recently studied 
in ( Amico and Patariel l200?t ) . In analogy with the zero 
temperature case they demonstrated that the entangle- 
ment sensitivity to thermal and to quantum fluctuations 
obeys universal T ^ 0-scaling laws. The crossover to the 
quantum disordered and renormalized classical regimes in 
the entanglement has been analyzed through the study of 
derivatives of the concurrence d\C and dxC . The ther- 
mal entanglement results to be very rigid when the quan- 
tum critical regime is accessed from the renormalized 
classical and quantum disordered regions of the phase 
diagram; such a 'stiffness' is reflected in a maximum in 
8tC at T ~ T cross . The maximum in the derivatives of 
the concurrence seems a general feature of the entangle- 
ment in the crossover regime. In this respect we mention 
that also the conc urrence between two Kondo impu rity 
spins discussed in (jStauber and GuineaLl200ll2006hh ex- 
periences the largest variation again in the crossover phe- 
nomenon. 

Due to the vanishing of the gap at the quantum crit- 
ical point, in the region T 3> T cross an arbitrarily small 
temperature is immediately effective in the system (see 
Fig. I17p . From analysis of the quantum mutual infor- 
mation (see Eq. [9])) it emerges that the contribution of 
the classical correlations is negligible in the crossover, 
thus providing the indication that such a phenomenon is 
driven solely by the thermal entanglement. It is interest- 
ing to study how the exsistence of the factorizing field h / 
affects the thermal pairwise entanglement (vanishing at 
zero temperature) . It is found that the two tangle still 
vanishes in a region of the h — T plane, fanning out f rom 
just hf, (white region in Fig.^ flAmico et aUBoOrjfr . In 
there entanglement if present, is shared between three or 
more parties. It is further observed that, in contrast to 
the analysis of the ground state, at finite temperature 
we cannot characterize the two separate phases of par- 
allel and antiparallel entanglement by. In fact, the two 
types of entanglement (though well defined also for mixed 
states) can swap by varying T and/or r. The exchange 
between parallel and antiparallel entanglement occurs in 
a non trivial way, that ultimately produces the reentrance 
of T2 seen in FigfTSl 

The common feature in all cases for which the existence 
of entanglement could here be proved is that both high 
temperatures and high values of magnetic field move the 
thermal states away from the region with non-zero en- 
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FIG. 17 The effect of temperature on the anomalies origi- 
nated from the critical divergence of the field-derivative of 
C(R) can be measured by dr[d a C(R)]. The density plot cor- 
responds to 7 = 1 and R = 1. T = T* and T — Tm are 
drawn as dashed and thick lines respectively. Maxima below 
T* are found at T M = (3T cross with (3 ~ 0.290 ± 0.005 and 
they are independent of 7 and R; the crossover behaviour 
is enclosed in between the two flexes of dT[d a C(R)] at T c i 
T c2 ; such values are fixed to: T c i = (0.170 ± 0.005)T cross 
and T c2 = (0.442 ± 0.005)T cross and found to be indepen- 
dent of 7 and R. For T < T cl d T [d a C{R)] ~ 0. Scaling 
properties are inherited in dT[d a C(R)] from d a C(R)] [From 
I Amico and PatanAl2007| |. 




FIG. 18 Contour plot of T2 in the h—T plane, for 7 = 0.3 (i.e. 
hf = 0.9539. . .). Th e white area indicates where ti = 0.[From 
lAmico et ail 120061 ] 



tanglement. This is understandable because high values 
of magnetic field tend to order all spins parallel to the 
field which corresponds to an overall state being a prod- 
uct of the individual spin states. There is upper limit of 
this phenomenon, since an increase in the temperature is 
disruptive for entanglement due to thermal fluctuations. 



C. Thermal entanglement witnesses 

At first sight it may be very surprising that thermo- 
dynamical variables can witness entanglement since the 
only need to obtain them are the system's eigenenergies, 
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and no eigenstate information is required. Since entan- 
glement resides in the fact that the state is inseparable 
(and is not related to the value of its energy) it would 
appear that the partition function is not enough to char- 
acterize entanglement. This logic, although seemingly 
simple, is not entirely correct. The reason is subtle and 
lies in the fact that the whole Hamiltonian is used for 
constructing the partition function, so in a roundabout 
way we do have the information about the states as well. 
Entanglem ent witnesses in spin systems have been con- 
sidere d in (iBrukner and Vedrall . 120041 ; iHide et all l2007t 



Heisenberg model, one has 



iTothl . 1200a IWu et all 120051 ) (see also the next Section 



on multipartite entanglement) 

We will now illustrate how and why we can derive en- 
tanglement witnesses for the partition function. Sup- 
pose that the system is described by an antiferromag- 
netic Heisenberg model. First we have the following 
identity coming from the Hellman-Feynman theorem 



EN I QX c 
i=l\°i D ; 

means that the two point correlations function summed 
over all nearest neighbors can be derived from the parti- 
tion function. This is also a quantity that can be accessed 
experimentally as is usually performed very frequently in 
the solid state experiments. Most importantly, this av- 
erage has in general different values for separable and 
entangled states. It can, therefore, serve as an entangle- 
ment witness as will be seen shortly. The amazing fact 
that will emerge is that, in order to say if a state is en- 
tangled, we really do not need to have the analytically 
calculated form of the eigenstates in order to tell if the 
resulting mixture is entangled. One price to pay for this 
is that we will only be able to derive a sufficient con- 
dition for entanglement that is typical of entanglement 
witnesses. Namely, we will be able to tell if for some 
conditions the resulting thermal state is entangled, but 
we will not be able to say with certainty that the state 
is not entangled if these conditions are violated. 
Using U=(H) and M* = Y,?=i(Sj) we obtain 
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The right-han d of Eq.(|6 2l) is an entanglement witness 
as shown in (Toth, 120051 ): for any separable state, that 
is, for any classical mixture of the products states: p = 
Efc w k pl ® p\ ® ■■• ® Pt, one has 



1 N 1 

^lE« 5 ^+i> + (SfSf +1 ) + (S?S? +1 ))\ < -. 



(63) 



This is also valid for any convex sum of product states 
(separable states). The upper bound was found by using 
the Cauchy-Schwarz inequality and knowing that for any 
state (S x ) 2 + (S y ) 2 + (S z ) 2 < 1/4. It is important to 
note that the same proof can also be applied if one con- 
siders XX model. We now give our thermodynamical 
entanglement witness: if, in the isotropic XXX or XX 



\U + hM z \ > N\J\/A, 



(64) 



then the solid state system is in an entangled state. The 
entanglement witness is physically equivalent to the ex- 
change interaction energy or, equivalently, to the differ- 
ence between the total (internal) energy U and the mag- 
netic energy —hM . From the tracelessness of the Pauli 
operators one can easily see that linrr^oo C/ — > 0. This 
means that the value of the internal energy as given by 
Eq. (|64|) should be defined relatively to the referent value 
of zero energy in the limit of high temperatures. In order 
to complete the proof we need to give an explicit example 
of a state that violates Eqs. (j6"3"|) (or the corresponding 
inequality for the XX model). This implies that that 
Eq. (|64p is indeed an entanglement witness and not just 
a bound that is trivially satisfied by any quantum state. 

As an example of such a state we take the ground state 
of the antiferromagnetic isotropic XXX Heisenberg model 
with zero m agnetic fi e ld. T he energy of this state was 
found to be (|Hulthenl . Il938l) : \E /JN\ = 0.44325 > 1/4. 
Furthermore, due to the symmetry of the XXX Heisen- 
berg Hamiltonian one has E Q /(3NJ) = (SfSf +1 ) = 



(SIS? 



-i/o = 



-.(S?S? +1 )o = -1.773/12 for every i. 



This im- 

plies that 1/N\ Eti({SfSf +1 )o + (5?5? + i)o)| = 0.295 > 
1/4. Therefore, Eq. ([64]) is an entanglement witness 
for the solid state systems described by XXX or XX 
Heisenberg interaction. 

We will now discuss various concrete models of spin 
interaction of which some are exactly solvable and for 
which dependence of internal energy U and magnetiza- 
tion M on temperature T and magnetic field h are known. 
This will help us to determine the parameter regions of 
T and h within which one has entanglement in the solids. 

We first consider XXX Heisenberg model with no 
magnetic field, in this case the magnetization vanishes 
and thermodynamical witness, Eq. ([64|) reduces to \U\ > 
N\J\/4. It was shown that concurrence C(l) is zero at 
any temperature in the ferromagnetic case and that it 
is given by C = ± max [0, \U\/(NJ) - 1/41 in the anti- 
ferromagnetic case (I Wang and ZanardiL l2002h . Thus C 
is nonzero if and only if |Z7|/(iVJ) > 1/4. This shows 
that the thermodynamical entanglement witness can de- 
tect entire bipartite entanglement as measured by con- 
currence. Furthermore, the fact that the value of the 
entanglement witness for the ground state is well above 
the limit of 1/4 suggests that entanglement could ex- 
ist and be detected by the thermodynamical witness at 
nonzero temperatures as well. In the presence of a fi- 
nite magnetic field the low temperature partition func- 
tion is given by Z = e W+M/4(! + e -^/ 2 N/^/2%pJ). 
Using this we obtain \U + h z M\/(NJ) = 1/4 and 
thus no entanglement can be detected in agree ment 
with ijAsoudeh and Karimipourl . 120041 : iPrattl . l2004h . 

The XX Heisenberg model with nonzero magnetic 
field is the most interesting as it is exactly s o lvable , 
the partition function was found in (|Katsural . Il962h . 
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Let us introduce the following dimensionless quanti- 
ties: b — h z /T and K = J/T (note a difference 
of factor 2 in the def initions of J and K with re- 
spect to (|Katsural Il962l )) and the function f(K,b,ui) — 
\J 2K 2 + 2K 2 cos 2ui — AbK cos uj + C 2 for convenience. 



Then the internal energy is given by dKatsural Il962h 



U_ 

N 



T 



f(K, b, u) tanh /(A', b, u)du, (65) 



4tt J a 

and the magnetization by ( Katsural Il962| ) 
M 1 r AK 2 cos 2 uj 



N 



2nJ f(K,b,u) 



tanh f(K,b,w)duj (66) 



both in ferromagnetic and antiferromagnetic case. 

We use Eqs. (|6"5)) and (JBSJ) to determine the parameter 
regions of temperature T and magnetic field h z for which 
entanglement exists in the solid state system. The criti- 
cal values of T and h below which entanglement can be 
detecte d is of the order of J, which can be as high as 10 
Kelvin (|Hammar et all , [l999f ). 

It should be stressed that the analysis based on the 
entanglement witness could be applied to any model 
for which we can successfully obtain the partition func- 
tion. This feature is the main advantage of using ther- 
modynamic witnesses approach to detecting entangle- 
ment. This method for determining entanglement in 
solids within the models of Heisenberg interaction is use- 
ful in the cases where other methods fail due to incom- 
plete knowledge of the system. This is the case when 
only the eigenvalues but not eigenstates of the Hamilto- 
nian are known (which is the most usual case in solid 
state physics) and thus no measure of entanglement can 
be computed. Furthermore, in the cases where we lack 
the complete description of the systems one can approach 
the problem experimentally and determine the value of 
the thermodynamical entanglement witness by perform- 
ing appropriate measurements. It is important to em- 
phasize that any other thermodynamical function of state 
could be a suitable witn ess, such as the magne tic suscep- 
tibility or heat capacity ( Wiesniak et all , |2005| ) (see next 
Section). 

The temperature as well as other thermodynamical 
state variables have been shown to behave also as entan- 



glement witnesses of spatial entanglement ([Anders et all , 
2006). This general feature was explicitly worked out in 
the case of a non-interacting Bosonic gas. It was found 
that entanglement can exist at arbitrarily high temper- 
atures, provided that one can probe smaller and smaller 
regions of space. 

The methods outlined here are not only applicable 
to the models we considered. There are several in- 
teresting questions and possibilities for generalizations 
such as consideration of Hamiltonians with higher spins, 
two- and three-dimensional systems, non-nearest inter- 
actions, anisotropies, other thermodynamical properties 
(e.g. heat capacity, magnetic susceptibility) and so on. 



Similar analysis can be done for continuous thermal en- 
tanglement in a field. It has been shown that for non- 
interacting bosons, entanglement exists when their de 
Broglie therm al wavelength i s smal ler than their average 
separation, a (|Anders et all l2006h . The precise condi- 
tion is kT < h 2 /2ma 2 , where m is the mass of bosonic 
particles. We can now introduce the following correspon- 
dence between spin coupling J and the continuous vari- 
ables bosonic kinetic energy, J = h 2 /2ma 2 . This fur- 
ther implies that we can think of the thermal de Broglie 
wavelength for spins as XdB = ay/ J/T, where a is the 
spin separation. The condition for entanglement that 
the wavelength is larger than the lattice spacing a now 
leads us to the condition that T < J which is exactly the 
result obtained from a more detailed analysis above. 



D. Experimental results 

The question of having macroscopic entanglement is 
not only fascinating in its own right but it also has a fun- 
damental significance as it pushes the realm of quantum 
physics well into the macroscopic world, thus opening the 
possibility to test quantum theory versus alternative the- 
ories well beyond the scales on which theirs predictions 
coincide. It also has important practical implications for 
implementation of quantum information processing. If 
the future quantum computer is supposed to reach the 
stage of wide commercial application, it is likely that it 
should share the same feature as the current (classical) 
information technology and be based on solid state sys- 
tems. It will thus be important to derive the critical 
values of physical parameters (e.g. the high-temperature 
limit) above which one cannot harness quantum entan- 
glement in solids as a resource for quantum information 
processing. 

Recently it was demonstrated experimentally that en- 
tanglement can affect macroscopic properties of solids, 
albeit a t very low (critic al) temperature (below 1 
Kelvin) (iGhosh et all , l2~003h . This extraordinary result 
opens up an exciting possibility that purely quantum cor- 
relations between microscopic constituents of the solid 
may be detected by only a small number of macroscopic 
thermodynamical properties. 

Ghosh et al made measurement on an insulating salt, 
the LiHoF4. At low temperatures the susceptibility de- 
viates from a simple Curie-like law. They find that the 
the temperature dependence is well fitted by a power law 
X ~ T~ a with a = 0.75. The key observation made 
in ( Ghosh et all , [20031 ) is that the experimental data at 
low temperatures cannot be explained by simply resort- 
ing to a classical approximation. By itself this might not 
be enough. It is remarkable however that the authors are 
able to put in close connection the power law divergence 
of the susceptibility with the entanglement present in the 
low-lying excited states. 

Vertesi and Bene studied the magnetic susceptibility of 
NaVs07 and used macroscopic entanglement witnesses 
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FIG. 19 Magnetic susceptibility vs. temperature T 
from simulations of the diluted, dipolar-coupled Ising mag- 
net, compared to experimental data (red triangles). The black 
circles use quantum decimation as well as the correct quan- 
tum mechanical form of susceptibility, utilizing the entangle- 
ment of the low-lying ene rgy doublet with the excited states. 
[From lGhosh et aUl2003l] 
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FIG. 20 The temperature dependence of magnetic suscepti- 
bility of powder CN (triangles) and a single-crystal CN, mea- 
sured at low field parallel (open squares) and perpendicular 
(open circles, crosses, filled circl es) to the monoclin ic b axis. 
The data and the figure are from ijBerger et oZ.l . ll963h . The in- 
tersection point of this curve and the experimental one defines 
the temperature range (left from the intersect ion point) with 
entan glement in CN around ~ 5K. [From iBrukner et all 
120061 ] 



as discussed in the previous Section to estimate the crit- 
ical temperature below which thermal entanglement is 
present. The experimental value of this temperature is 
365K approximately three times higher than the critical 
temperature c orresponding to the van ishing of bipartite 
entanglement (jVertesi and Bend . l2006l) . 

We also mention the very recent experiment 
reporting on macroscopic magnetic measure- 
ments of the pyrob orate MgMnB 9 2C >5 and war- 
wic kite MgTiOBO,, (jRappoport et all 120071 ) (see 
also IContinentinol . 120061 ) 

Earlier experimental data witnessing entanglement 
in bulk properties of solids have been reanalyzed 
in ( Brukner et all 120061 ). They discusses the experi- 
mental results of n eutron scatterin g measurement of CN 
obtained in 2000 ( Xu et all , l2000f ) and show that they 
provided a direct experimental demonstration of macro- 
scopic entanglement in solids. The experimental char- 
acterization of the dynamic spin correlations for next 
neighboring sites enabled them to determine the concur- 
rence and show the existence of entanglement at mod- 
erately high temperatures (as high as 5 Kelvin). In 
the same work they also showed that magnetic suscepti- 
bility at zero magnetic field is a macroscopic thermo- 
dynamical entanglement witness for the class of solid 
states systems that can be modeled by strongly alter- 
natin g spin- 1/2 antiferromagnet chain ( Brukner et all 
200G) .The measured values for magnetic susceptibility of 
CN in 1963 ( Berger et aZ.I . fl963l ) imply presence of entan- 
glement in the same temperature range (below 5 Kelvin). 

An analysis of the experimental r esults of a m a gnetic 
susceptibility measurement of CN ( Berger et all Il963f) 
showing that the values measured at low temperatures 
cannot be explained without entanglement being present 



was performed in ( Brukner et all I2006T ) . This was based 
on the general proof that magnetic susceptibility of any 
strongly alternating antiferromagnetic spin-1/2 chain is a 
macroscopic thermodynamical ent anglement witness. As 
discussed in ( Brukner et all [20061 ) the magnetic suscep- 
tibility for separable states is bounded by the value 



Xsep 



> 



T 6' 



(67) 



The results of their analysis are reported in Fig|2TJl 



VIII. MULTIPARTITE ENTANGLEMENT 

Most of the results reviewed in the previous section 
are for pairwise entanglement. Although much has been 
learned from the study of those quantities, the structure 
of entanglement in many-body systems is much richer 
and it is natural to expect that multipartite entangle- 
ment is present both in the ground state and at finite 
temperatures. Although multipartite entanglement in 
many-body systems is much less studied, some impor- 
tant results have been already obtained. 

A number of groups showed that in certain limits the 
state of a spin chain can resemble that of known multi- 
partite states. For small chains Wang first noted that 
the gr ound s t ate ten ds to have multipartite entangle- 
ment (|Wand . l2002al) . This analysis was further pur- 
sued in Ising a nd Heisenberg rings where the ground 
state has GHZ- 



like (Brufie£ al\ 



Stelmachovic and Buzekl 120041) and W- 



20051 ) correlations, respectively. Multi- 



partite mixed states can be also realized in the case of 
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FIG. 21 One-tangle n and the sum ol squared concurrences 
T2 as a lunction of the applied magnetic field in ^-direction for 
the XY Z model with exchange along y: J y = 0.25 (in unit 
of exhange along z). Inset: contributions to the concurrence 
between j-th neighbors; full symbols stand for G 1 (J), open 
symbols f or C II (j). The dashe d line marks the critical field 
h c . [From lRoscilde et a/.L 12004} 



ferromagnetic rings where the ground state is fully polar- 
ized along the direction of the field and the first excited 
state is a W-state. At finite but low temperatures the 
density matrix is approximately as p — p | • • • 0) (0 • • • | 
+ (1 — p) | W r )(VF | where p is the Boltzman factor and 
| • • ■ 0) is the ferromagnetic ground state. For three 
qubits p has been shown to contain tr ipartite entangle- 
ment in (jBrufi and MacchiavellJ l2003fh 

These examples show that several models naturally 
have multipartite entangled ground states. At the same 
time it also shows that it is important to quantify mul- 
tipartite entanglement in many-body systems. Analysis 
in this direction is reviewed below. 



A. Multipartite entanglement in spin systems 

A first way to estimate multipartite entanglement in 
spin system is provided by the entanglement ratio t 2 /ti 
as the amount of two spin relative to global entanglement. 
For Id — XY Z models it was shown that a small value 
of such a ratio is generic in these syste ms with a cusp 
minim um at the quantum critical point ( Roscilde et all 
2004). This is was shown numerically for the XYZ chain 
in a field (FigHTJ. 

In order to quantify the multipartite entanglement, 
Wei et al. calculated the geome tric measure of entangle- 
ment ( Wei and Goldbarq 120031). see III.El for the trans- 
verse XY chain ()Wei et all l2005h . The calculation can 
be done by a clever use of symmetries: translational in- 
variance and periodic boundary conditions. In this case, 
the set of all possible separable states can be described 
by a global rotation around the y-axis of the fully polar- 
ized state. The maximization is thus reduced to only one 
variable. 

As well as in the case of bipartite entanglement, also 
the multipartite measure of Wei et al. is very sensitive to 
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FIG. 22 The geometric entanglement measure Eq.fT2]) per 
site is plotted versus 7 and h z for the XY model. There 
are three different phases: the ordered oscillatory phase for 
7 2 + h\ < 1 and 7 / 0; the ordered ferromagnetic phase 
between , y 2 + h)z 2 > 1 and h < 1; and the paramagnetic phase 
for h z > 1. There is a sharp rise in the entanglement across 
the quantum phase transition line A" 1 — h z = 1. At 7 = 
lies the XX model, which belongs to a diffe rent universalis 
class than the anisotropic XY model. [From 



ircnt universality 
I Wei et all [20051 ] 



the existence of QPTs. As a paradigmatic example the 
authors analyzed the phase diagram in the anisotropy- 
magnetic field plane. Their results are shown in Figl2"2"l 
As expected there is no divergence in the measure itself 
but in its derivative. The new aspect here is that dif- 
ferently from the concurrence the multipartite entangle- 
ment measure in Eq. (|12p can be expanded as a function 
of n-point correlators. 

The geometric entanglement cannot discriminate be- 
tween different n-particle entanglement classes. A com- 
prehensive classi fication in spin sy stems has been re- 
cently given by (jGuhne et all [2005) via the concept of 
k-producibility fsee lII.E[ ). The systems analyzed in detail 
are the one-dimensional XY and Heisenberg models. Dif- 
ferent types of n-particle quantum correlated states lead 
to distinct bounds for the internal energy (or the ground 
state energy at T = 0). A violation of these bounds then 
implies the presence of certain k-party pro ducable entan- 
glem ent in the system. As pointed out in ( Giihne et all 
2005), aiming at the thermodynamic limit of an infinite 
number of spins, the notion of k-producibility is easier to 
handle than the n-separability (see III. E| ) as its definition 
is independent of the number of sites in the system. 

A systematic approach for deriving energy bounds for 
states without certain for ms of multipartite entang lement 
has been developed in (jGuhne and Tothl . I2006T ). The 
method allowed to investigate also higher dimensional 
and frustrated systems. As an example we report on the 
results for the Heisenberg model. In d dimension, the 
energy bounds per bond for one-, two-, three-, and four- 
producible states are given by 
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FIG. 23 Entanglement in thermal states of the XX-model in 
a transverse magnetic field. The different regions correspond 
to different types of multipartite ent anglement contained in 
the equilibrium thermal state. [From lGuhne and TothL l2006j ] 
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It is striking how relatively close the 2- and 3-producible 
bounds are in all cases. All the bounds given above are 
found to be violated in the ground state. In the previous 
expression the superscripts refer to the dimensionality of 
the model and the subscripts to the k-party entanglement 
for which the bound is obtained. There is a factor 1/4 of 
difference with respect to the original paper because of 
the different notation used in this review. 

Corresponding to the energy scales fixed by the bounds 
there are different temperature scales at which the vari- 
ous n-party entanglement types disappear. As one would 
expect, these crossover temperatures are monotonic in k; 
i-e. Tfc > Tfc+i. The example given in Fig|2"31 clearly 
shows how higher order multipartite entanglement pro- 
gressively disappears as the temperature increases. 



B. Global entanglement 

Multipartite entanglement close to quantum phase 
transitions was quantified by the global-entanglement 
Epi measure of Meyer a nd Wallach ( s ee III. ED in 
(|de Oliveira et all l2006allbl ; ISomma et all . |2004) . To- 
gether with the Meyer- Wallach measure, de Oliveira et 
al also introduced slight extension of it as 
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(68) 



where Pjj+i is the reduced density matrix associated to 
the sites j and j+l. Similarly one can consider also three- 
body reduced density matrices and construct the corre- 
sponding global entanglement measure. According to de 
Oliveira et al. the above hierarchy might provide a com- 
prehensive description of entanglement in many-body 



systems. E xplicit calculation s for th e one-dimensional 
Ising model (|de Oliveira et HI. I2006ah showed that both 



(2) 

E g i and E gl are maximal at the critical point (with 

(2) 

Egi < E gl for any value of the coupling constant) sug- 
gesting that the quantum critical point is characterized 
by the presence of multipartite entanglement. As in the 
case of concurrence the non-analyticity associated to the 
critical point is manifest in the derivative of the global en- 
tanglement m easure. By exten ding an earlier approach 
developed in ( Wu et all l2004h . de Oliveira et al. also 



(n) 

showed how the non-analytic behavior of E gl is related 
to that of the ground state energy. Note that from 
Eq. (f6"8"|) it is possible to define an entanglement length 
proportional to the correlation length £. This differs con- 
siderably from that one defined by the localizable entan- 
glement (see Eg. ((SSI) 1 ); the latter is always bounded from 
below by the correlation length and can even be divergent 
where £ is finite. 

As discussed in ( Facchi et all l2006aH bT) the analysis of 
the average purity might not be sufficient and the analy- 
sis of the distribution of the purity for different partitions 
could give additional information. Rather than measur- 
ing multipartite entanglement in terms of a single num- 
ber, one characterizes it by using a whole function. One 
studies the distribution function of the purity (or other 
measures of entanglement) over all bipartitions of the sys- 
tem. If the distribution is sufficiently regular, its aver- 
age and variance will constitute characteristic features of 
the global entanglement: the average will determine the 
"amount" of global entanglement in the system, while 
the variance will measure how such entanglement is dis- 
tributed. A smaller variance will correspond to a larger 
insensitivity to the choice of the bipartition and, there- 
fore, will be characteristic for different types of multipar- 
tite entanglement. The application of this technique to 
the one-dimensional quantum Ising model in a transverse 
field shows that the distribution function is well behaved 
and its average and second mo ment are good indicator s 
of the quantum phase transition (jCostantini et al. I. l2007lh 
This is in agreement with previously obtained results. At 
the onset of the QPT both the average and the standard 
deviation exhibit a peak that becomes more pronounced 
as the number of qubits is increased. 



C. Generalized entanglement 

A different route for studying the multipartite en- 
tanglement is paved by the general observable based 
entanglement (see Section III. El) . It was first pursued 
by Somma et al for the LMG and the quantum XY- 
model. In the realm of solid state systems an experi- 
mental protocol to measure many-fermion entanglement 
based on this concept has been proposed in ( Kindermannl . 
120061) . A nice connection which emerges from the work 
of Somma et al is the one between the generalized entan- 
glement and the quantum fluctuations of the magnetiza- 
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FIG. 25 The purity P' u 



FIG. 24 Purity relative to the observable S z in th e ground 
state of the LMG model. [From lSomma et q/.l I2004T ] 



tion which are responsibl e for the quantum phase transi- 
tion ( Somma et all [20041 ). Later, Montangero and Viola 



u(N) 



considered the dynamical behavi or of generalized entan- 
glem ent in disordered systems (iMontangero and Violal . 
I200fih . As remarked in (jSomma et aZJ . l2004h . it is impor- 
tant to choose the appropriate subset of observables in 
order to see the critical behavior in the entanglement. 

In the case of the LMG model a natural choice for the 
local observables was to consider the purity relative to 
the spin-iV/2 representation of the angular momentum 

P s = [(S x ) 2 + (S y ) 2 + (S z ) 2 ] . Somma et al. also 

consider the purity relative to the single observable S z : 
4 

P s * = — x \(S Z ) 2 ] . With this last choice the relation be- 
N 

tween the multipartite entanglement and the order pa- 
rameter becomes evident. The result is shown in FiglMl 
Further interesting results emerge from the XY model. 
By choosing the following subset of operators ex- 
pressed in terms of the spinless fermions of the Jordan- 
Wigner transform u(N) = span{ c\cii + c\,Ci,i(c\ci> — 



etc. 



c\ci — 1/2)} it is possible to show that the as- 



sociated purity may be considered as a disorder param- 
eter, i.e. it is non-zero in the symmetric phase only. A 
transparent way to express the purity in this case is to 
relate it to the variance of the number fermion opera- 
tor Nf — J2i c l c i- The result is plotted in Fig|2"5l for 
different values of the anisotropy parameter 7. Two con- 
siderations are in order at this stage. First of all it is clear 
from the previous example that important properties of 
entanglement appear when one moves away from the con- 
ventional picture of partitioning in real space. Secondly 
the purity shows, differently from other measures as the 
concurrence, direct indications of the critical behavior 
(and not only in the derivatives). 

D. Renormalization group for quantum states 

We continue our discus sion on multipartite en tangle- 
ment with a recent work ( Verstraete et all l2006f ) where 



3 u(jv) — (1+7) - of an XF model 
in a transverse field as a function of A for different anisotropies 
7. The constant part has been subtracted in ord er to make the 
conn ection to the disorder parameter. [From ISomma et all 
120041 ] 



the method of renormalization group is applied to quan- 
tum states and not, as traditionally done, to the Hamilto- 
nian. The Renormalization Group is based on a recursive 
transformation which leads to an effective description of 
the low-energy (long distance) physics of a given system 
where all the effects of high energy modes are included 
in a renormalization of certain coupling constants of the 
model Hamiltonian. A study of the flow of these recur- 
sive equations, their fixed point(s) and the behavior of 
the flow close to the fixed points allows to determine 
the critical behavior of the system under consideration. 
Contrary to the conventional renormalization group ap- 
proach, the idea of Verstraete et al. is to analyze how the 
quantum states change under this coarse graining and to 
classify their fixed points. Given a system characterized 
by N sites and a ZJ-dimensional local Hilbert space. The 
steps of this real space renormalization procedure are the 
following. 

(i) Coarse graining, in which clusters of neighboring sites 
are merged into one site of a new lattice. 

(ii) Rescaling of the local variables associated to the new 
sites. 

(iii) Identification of the states which are invariant un- 
der local unitary transformation (long distance behavior 
should not depend on the choice of the local basis). 

(iv) Rescale the distances in order to have again a unit 
lattice constant. 

In this way a representative of a given class of quantum 
states, invariant under local transformations, is trans- 
formed at each step. The irreversibility of the transfor- 
mation is crucially related to the step (iii) of the transfor- 
mation as one loses track of the unitary transformation 
performed (which need to be local over the cluster). 

Verstraete et al carry out explicitly the renormalization 
group transformation for Matrix Product States (MPS), 
see Eq. (|29p . The decimation step in which two neighbor- 
ing sites, say i and are merged together is performed 
here by means of merging the corresponding matrices A Pi 
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means of a singular value decomposition of A it is possible 
to find the representative of the new state (see the step 
(hi)) and therefore to iterate the renormalization group 
map. In the case of D^ps — 2, Verstraete et al. provide 
a complete classification of the fixed points of the trans- 
formation which contains product, GHZ and domain wall 
states. A similar analysis in the case Dmps > 2 and/or 
in higher dimensional systems may turn out to be very 
useful also for the classification of multipartite entangle- 
ment in many-body systems. 



E. Entanglement distribution for Gaussian states 



quantity which is analog to the pairwise tangle has been 
identified as the square of the logarithmic negativity 



It has been observed in Section III.GI that the symme- 
try groups of admissible local transformations of Gaus- 
sian states and qubits are isomorphic. This suggests to 
look for further analogies with the qubit case or to search 
for deviations from it. A striking feature of qubit en- 
tanglement is the monogamy as far as entanglement dis- 
tribution along chains of qubits is concerned (see III.Bjl . 
For continuous variable entanglement such an inequal- 
ity was originally analyzed analytically for fully symmet- 
ric Gaussian states and numerically for randomly chosen 
Gaussian states in ( Adesso and I lluminati. 2005a). The 
general p roof of monogamy has been obtained very re- 
cently in ( Hiroshima et al. 2007). In this Section we use 
the particular case of a symmetric state as a guideline 
to discuss the monogamy for Gaussian states. A more 
detailed discu ssion and the general p r oof ca n be found in 
the review by lAdesso and Illuminati 120071 

A Gaussian state is called fully symmetric if and only 
if its covariance matrix is invariant with respect to per- 
mutations of t he modes. Its covari ance matrix can then 
be written as (jAdesso et all [2004a) 



v„ 



e (3 



e e 



\ 
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where j3 and e are 2x2 matrices. Both can be diag- 
onalized by means of local symplectic transformations 
in phase space, such that $ — diag {b, b} and e = 
diag {ei, £2}- This leads to a highly degenerate symplec- 
tic spectrum where n — 1 doubly degenerate symplectic 
eigenvalues are equal to v — (b — Si)(b — £2) and the re- 
maining eigenvalue is v n — (b+ (n— l)ei) (6 + [n— 1)£2)- 
After a partial transposition, n— 2 symplectic eigenvalues 
v remain unaffected. The negativity is then determined 
solely by which is the s maller one of t he rem aining 
two affected eigenvalues v± ( Adesso et all l20"04al) . This 
is due to the fact that the uncertainty leads to a lower 
bound h 2 /A for the product v + V- 

For a single Gaussian mode with covariance matrix V\ , 
the tangle is given by t\ — 2(1 — det V^/Vdet V\ The 



t 2 (£>__) = [In; 



(70) 



and has been ter med contangle 

in ( Adesso and Illuminatil . l2005al ). The identifica- 
tion of the square negativity as the continuous variable 
tangle is crucial for the demonstration of the monogamy 
inequality. Extensions of these measures to mixed states 
are un derstood to be g iven by the corresponding convex 
roofs ( Uhlmannl . [2000), where the average of the pure 
state measure of entanglement has to be minimized over 
all pure state decompositions of the density matrix in 
consideration. A restriction to decompositions purely 
out of Gaussian states gave rise to the notion of the 
Gaussian entanglement measures (|Wolf et all 120041) 
and the Gaussian entanglement of formation. It clearly 
establishes an upper bound for the entanglement of 
formation. 

Negativity and Gaussian entanglement measures 
have been analyzed with considerable detail in 
( Adesso and Illuminatil . l2005hl ) for two-mode Gaussian 
states with particular focus onto symmetric Gaussian 
states. One important result is that negativity and 
Gaussian measures lead to different ordering of Gaussian 
states with respect to their entanglement. For symmet- 
ric Gaussian states instead, both give the same ordering. 
This result must be handled with care, since it is not 
obvious what precisely the restriction to Gaussian de- 
compositions entails for the ordering of states. Believing 
into the striking similarity to qubit systems, one might 
be tempted to conjecture that for symmetric Gaussian 
states a purely Gaussian optimal decomposition always 
existed. 

There are two extremal classes of two-mode Gaussian 
states considered: for fixed local and global purity, those 
states that maximize the negativity are termed Gaus- 
sian Maxim ally Entangled Mixed States (GMEMS) in- 
troduced by lAdesso et M l2004bL whereas those states 
minimizing the negativity have been termed GLEMS, 
which are states whose covariance matrix has one sym- 
plectic eigenvalue equal to 1 (mixed states with partial 
minimal uncertainty). The symmetric two-mode Gaus- 
sian states have been singled out as those states con- 
tainin g minimal Gaussian entangle ment at fixed negativ- 
ity ( Adesso and Illuminatil . l2005bl ). 

The entanglement sharing inequality for Gaussian 
states assumes the same form as for qubits n,, > 
yijj.j T 2,ij where the indices are the numbers of the var- 
ious distinguishable modes in a multi-mode Gaussian 
state. This inequa lity has been proved f or multi-mode 
Gaussian states by iHiroshima et al. I. I2007L 
In the particular case of symmetric states all ri j = t\ \ 
and T2,ij = T2 12 for all indices i, j. 

As for qubit systems one can define a residual contan- 
gle by the difference r n := T\.i — Yljjii T 2,ij as a quanti- 
fier of the multipartite entanglement present in a Gaus- 
sian state. In particular is the residual contangle an 
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entanglement monotone under Gaussian LOCC. An im- 
portant difference, however, arises as compared to the 
residual tangle for qubits: not even for three modes is 
the residual tangle invariant under permutations of the 
modes and its minimum respect to the common mode i 
has to be taken. Even for symmetric Gaussian states, 
where this anomaly is absent, a promiscuous nature 
of continuous variable entanglement is encountered, in 
sharp co ntrast to the monogamy inherent to qubit entan- 
glem e nt dCoffman et all [20001 : lOsborne and Verstraetd . 
120061 ; lYu and Sons! . 120051 ). The term promiscuous is an 



interpretation of the fact that the maximal residual con- 
tangle T3 in a symmetric Gaussian three-mode state with- 
out pairwise contangle t-i is smaller than the maximum 
possible residual contangle without this restriction. Hav- 
ing in mind the entanglement sharing inequality, this im- 
plies that the local contangle t± is larger, when T2 and T3 
coexist. It is worth noticing at this point that the pecu- 
liarity of Gaussian states is to be completely described by 
two-point correlation functions. Consequently, all types 
of multipartite entanglement are inextricably related in 
that the same type of correlations are responsible for ei- 
ther type of quantum correlation. This is not the case for 
non-Gaussian states and in particular not for qubit sys- 
tems or general higher-dimensional local Hilbert spaces. 



IX. DYNAMICS OF ENTANGLEMENT 

The interest in studying the properties of entanglement 
in many-body system has been recently directed also to 
the understanding of its dynamical behaviour. Entangle- 
ment dynamics has been studied from different perspec- 
tives. In a spirit similar to the study of propagation of 
excitations in condensed matter systems, several works 
analyzed the propagation of entanglement starting from 
a given initial state where the entanglement has been cre- 
ated in a given portion of the many-body system. One 
can imagine for example to initialize a spin chain such 
that all the spins are pointing upwards except two neigh- 
boring spins which are entangled. Due to the exchange 
interaction the entanglement, initially localized on two 
neighboring sites of the chain will spread. This propa- 
gation is ballistic in clean systems. A "sound velocity" 
for the excitations results in a finite speed for the prop- 
agation of excitations. If some weak disorder is present 
one might expect diffusing propagation. Entanglement 
localization and chaotic behaviour can eventually also 
be observed. An alternative approach is to start with 
the ground state of an Hamiltonian Ho and the let the 
Hamiltonian change in time. Most of the attention up to 
now has been devoted to the case of sudden quench, i.e. 
after the preparation the Hamiltonian suddenly changes 
to Hi. Moreover since we are mostly dealing with in- 
teracting systems, entanglement can be also generated 
or it can change its characteristics during the dynamical 
evolution. 

Another important aspect of entanglement dynamics 



is the possibility to generate entangled states with given 
properties by taking advantage of interaction present in 
a many-body systems. This is the natural generalization 
of the case where a Bell state can be obtained by letting 
two qubit interact for fixed time by means of an exchange 
coupling of XX symmetry. In the same spirit one can 
think to generate three-bit entangled GHZ or W states or 
multipartite entangled states by tailoring the appropriate 
exchange couplings in spin networks. 



A. Propagation of entanglement 

1. Pairwise entanglement 

The simplest situation, which we consider first, is the 
propagation of entanglement in the one-dim ensional XX- 
mode l i.e. 7 = and A = in Eq. (|2"3"|) ( Amico et all 
12004 ISubrahmavaml |2004J) . Suppose that the initial 
state of the chain is 



f 



|* ± (i = 0)> = — (of ±o?)|0,. 



.0), (71) 



namely all the spin are in a fully polarized state except 
two, at positions i and j, which are in one of the two Bell 
states |V>±) = 2 _1 / 2 (|01) ± |10)) In this case the problem 
is amenable of a simple analytical solution. The reason 
is that, since the total magnetization is conserved, the 
evolution is confined to the sector where only one spin is 
up. In this sector the dynamics is completely described 
by the states \l) = |0, . . . 0, 1, . . . 0) (I = 1, . . . , N) which 
represents a state of the chain where the Zth spin is pre- 
pared in |1) while all the others N — 1 ones are in |0). At 
later times the state of the chain is to be 



l*±(*)> 



w ±,i 



W) 



(72) 



In the thermodynamic limit, N — > 00, the coefficients 
can be expressed in terms of Bessel functions J n {x) as 
w l lf (t) = ^ {Ji-iiAJt) ± H)i- i J i _ i (4Jt)} . Eq.cza 
together with the coefficients defined previously allows 
to study various kinds of entanglement measures for this 
particular case. 

As a first indication of the entanglement propagation 
we analise the dynamical evolution of the concurrence 
between two sites, located at positions n and m, (initially 
the entangled state is between the sites i and j). The 
concurrence reads 



(73) 



The function C^ m (t), plotted in Fig|2T)l shows for sites 
which are symmetrical respect to the initial position of 
the Bell state \i/j±). The time evolution dictated by the 
XX Hamiltonian, amounts to a propagation of the single 
flipped spin through the chain. The speed of the prop- 
agation is set up by the interaction strengths. The in- 
formation exchange or entanglement propagation over a 
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FIG. 26 Concurrence between sites n = ~x,m = x, symmet- 
rically placed with respect of the state where the singlet was 
initiall y created, from thei r initial position i = — 1 and j = 1. 
[From (|Amico et all 120041 )] . 



distance of d lattice spacings, approximately takes a time 
t ~ Hd/ J. The external field h does not enter in Eq. (|73|) . 
since all of the components of the state are in the same 
sector where one spin has been flipped and therefore it 
contributes only to an overall dynamical phase to the 
state evolution. 

Having in mind the use of spin chains as quantum chan- 
nels, the preparation scheme described above does not 
lead to faithful entanglement transfer. The most nat- 
ural way to perform this task is to prepare the entan- 
gled state and then let only one of the qubit to propa- 
gate along the spin chain (ther eby followin g the protocol 
originally pr oposed bv Bose (iBosd. l2003h or its modi- 
fied version ( Christandl et all |2004[ ) to achieve perfect 
transfer with modulated chain). A detailed analysis in 
thi s direction was recent ly performed by Hartmann et 
al (|Hartmann et a/.l . I2006T ). They considered initial max- 
imally entangled states and used the chain to transfer the 
state of one of the two qubits, found a relation between a 
measurement of the entanglement fidelity at the fidelity 
of state transfer and concluded that there are possibility 
to have perfect entanglement transfer. 

If the chain is initially prepared in |^±(f = 0)), given 
in Eq. (|71[) . the entanglement propagates maintaining its 
original characteristics. This is not the case if, for exam- 
ple, the initial states of the two entangled sites is of the 
type |$±) = (1/V2)(00) ± |11)). These states are super- 
positions of components belonging to different magneti- 
zation sectors. The entanglement propagates with the 
same velocity as before, however, under certain condi- 
tions, it turns out that, the propagating quantum corre- 
lations change their character. After some time a singlet- 
like entangled state propagates even if the initial state 
was not a singlet. A di fferent set of initial states has 
also been considered. In ( Hamieh and Katsnelsonl l2005f ) 
the chain was initialized in a separable state by means 
of a set of projective measurement (in particular they 
considered all the spins aligned in the z-direction except 
one prepared in an eigenstate of S x ). The evolution can 
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FIG. 27 The nearest neighbor logarithmic negativity Em as 
function of the initial transverse field h z and time t for the 
anisotropy 7 = 0.5. At short times t < 1.8 the dynamical 
phase transition is a point of reentrance for the entanglement. 
At later times the state remains separable for values of the 
field above the critical value. For the t = 0, entanglement 
vanishes, as it should, as h z grows. The transverse magne- 
tization of the evolved state does no t show a similar crit ical 
behavior as a function of h z . [From l|Sen(De) et all [2005)]. 



be described using the same approach as outlined above. 
The new ingredient now is the creation of entanglement 
during the dynamics. Hamieh and Katsnelson showed 
that the oscillations of the entanglement wave has the 
same periodicity, but out of phase, with the oscillation of 
the (non-equilibrium) average magnetization. The distri- 
bution of entanglement in the chain has been analyzed 
after launching a si ngle excitation from the central site 
of a XX chain in (Fitzsim ons and Twamlevl . l2005h . It 
was found that the second-order moment of the spatial 
extent of the concurrence grows much faster (with a rate 
increasing as i 5 / 2 ) if some disorder is present in the cen- 
tral portion of the chain (in the ordered case the increase 
goes as t 2 ). 

Additional features emerge in the quantum XY with 
7 7^ 0. In this case the magnetization is not any 
longer a constant of motion (two spins can be flipped 
simult aneously). The ca lculations were done analyti- 
cally (|Amico et all I2004D resorting to the exact calcu- 
lation of the nec essary set of out of equilib rium corre- 
lation functions ( Amico and Osterlohl |2003 ). The most 
notable difference in the two-site entanglement, studied 
is an entanglement production from the vacuum. This 
occurs uniformly along the chain and it superimposes 
to the entanglement wave associated discussed before. 
The velocity of propagation of the entanglement is al- 
most independent on the anisotropy parameter 7. What 
is strongly dependent on 7 is the damping coefficient of 
the entanglement wave. As the anisotropy approaches 
7=1, the Ising point, the wave is strongly damped and 
vanishes approximately after few ~ J -1 . Furthermore 
in the anisotropic case it is possible to observe a non- 
trivial dynamics when the external magnetic field is time- 
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de pendent (Irluang and Kail . 120051 : ISen(De) et ail , |2005|) . 
In ( SenfDel et all 2005 ) the system is initially prepared 
in the ground state; the evolution is analyzed after the 
magnetic field is (suddenly) switched off. Sen(de) et al. 
analyzed the two-site entanglement by means of the log- 
arithmic negativity as shown in Fig[27|and demonstrated 
the existence of a dynamical phase transition, not observ- 
able in the magnetization. As it can be seen in Fig|2"7l 
at a fixed time t the entanglement vanishes (and then 
grows again) for certain values of the initially applied 
magnetic field. For t < 1.8 the critical field is almost 
independent on the time t and is h ~ 0.8. A remarkable 
different situation occurs for t < 1.8, here a dynamical 
phase transition occurs where the entanglement vanishes 
for any value of ft. > h c . For values of the initial field in 
the region near the phase transitions, entanglement be- 
haves non monotonically with respect to temperature of 
the initial equilibrium state. 

The two-site entan glement is non er- 
godic (ISen(De) et all l2004h . The evolution of two-site 
entanglement of a XY chain has studied after a sudden 
change of the external magnetic field. The evolved 
state does not approach its equilibrium value. That is, 
entanglement by itself does not relax to its equilibrium 
value in the absence of further contact with reservoirs. 
Therefore entanglement in such systems cannot be 
described by equilibrium statistical mechanics. 

The entanglement dynamics has been studied, 
to a large extent analytically, also in the LMG 
model (|Latorre et all [2005al: IVidal etaH . l2004tj ) . They 
considered the dynamical evolution of entanglement by 
analyzing both the one-tangle T\(t) and the concurrence 
C(t) for the cases in which the initial state is fully po- 
larized either in the z- or in the x-direction. Because of 
the symmetries of the model and since the initial state 
belongs to the sector with maximal spin, S = N/2, both 
these quantities do not depend on the spins which are se- 
lected. This means that they can be expressed solely in 
terms of the average value of the total spin (S a ) and its 
correlation functions (S a S 13 } (a = x,y,z). An interest- 
ing feature of this model is that its level spacing is finite 
even in the thermodynamic limit. Contrary to the expec- 
tations, however, the Poincare time is dependent on the 
number of spins N and this has important consequences 
on the entanglement dynamics. 

In the case in which the spins are prepared in 
the state |0...0) analytical results can be obtained 
in the limit of zero and very large magnetic field. 
In the limit of zero external field the tangle and 
the concurrence can be expressed as and Ti(t) — 
1 - c^) 2 ^- 1 ) and C(t) = (1/4) (c^" 2 ) - 1 + 

[ c ( t )(JV-2) _ 1]2 + [ 4cW (JV-2) sW (iV-2) ] 2l/2 ); ^ ^ = 

cos(4Ji/7V), s(t) = sin(4Ji/iV), and show perfect anti- 
correlation in time. In the opposite limit of a strong 
applied field it is possible to resort to a semiclassi- 
cal a pproximati o n of the Heisenberg equation of mo- 
tion (|Law et all l200lf) . S x = 2hS y , S y = -2hS x + 
(2J/N)[S Z , S w ]+, S y = -(2J/N)[S Z , S y ]+ by noting that, 



for N ^> 1, the z-component of the magnetization has 
negligible fluctuations (S z (t) ~ S = N/2). There- 
fore the set of equations governing the dynamics of the 
system becomes linear and can be easily solved. The 
concurrence as a function of time for positive values 
of the external field evolves as C(t) = 1 — Ch(t) 2 — 
( J 2 /4hl)s h (t) 2 with c h {t) 2 = cos 2 (^h(h- J)t) and 
Sh(t) 2 = sm 2 (y / h(h — J)t). Dynamics of two-site entan- 
glement was disc ussed also in the context of spin-boson 
like models. In (jCiancio and ZanardiL l2006f) the nega- 
tivity for a two-modes Jaynes-Cummings was analyzed 
with particular emphasis on the entanglement between 
the two bosonic modes mediated by the qubit. The re- 
laxation dynamics of the entanglement, quantified trough 
the concurrence, in the spin-boson model was discussed 
in (|Lim et aLl . l2006h 

We finally mention a study where it was noticed that 
the entanglement encoded in the states caused by the 
splitting of the degeneracy determined by the transverse 
field in the quantum XY model is not preserved by an 
adiabatic perturbation. Separable states can become en- 
tangl ed after the field is varied adiabatically, and vicev- 
ersa (|Cao et aLl . l2006fh 



2. Dynamics of the block entropy 

The dynam ical behaviour of the block en tropy was first 
considered in (jCalabrese and Cardvl . [20051 ) in the general 
framework of conformal field theory and via an exact so- 
lution of the quantum Ising model. Later Diir et al con- 
sidered Ising models with long range interaction and De 
Chiara et al (|De Chiara et all I2005T ) performed numeri- 
cal simulations of the XX Z chain. 

Calabrese and Cardy showed that a quench of the sys- 
tem from a non critical to a critical point leads the block 
entropy first to increase linearly in time and then to satu- 
rate. For periodic boundary conditions and given a block 
of dimensions I, the time at which the entropy saturates 
is given by t* — £/(2v) where v is the spin wave velocity. 



S e (t) 




t < t* 
t > t* 



(74) 



Thus there is an arbitrary large entanglement entropy in 
the asymptotic state, contrarily to the ground-state case, 
where the entropy diverges only at critical point. An ex- 
ample of the time dependence of the block entropy for the 
Ising and XX Z models is shown in Figl28l Calabrese and 
Cardy proposed a simple interpretation of this behaviour 
in terms of quasiparticles excitations emitted from the 
initial state at t = and freely propagating with veloc- 
ity v. The argument goes as follows. The initial state 
has a very high energy relative to the ground state of the 
Hamiltonian which governs the subsequent time evolu- 
tion, and therefore acts as a source of quasiparticle exci- 
tations. Particles emitted from different points are inco- 
herent, but pairs of particles moving to the left or right 
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FIG. 28 Evolution of the entropy for different types 
of quenches in the c ase o f the Ising (left panel) [from 
(ICalabrese and Cardvl 120051')] and Heisenberg (right panel) 
[from (|De Chiara et all l2005r )]models. The block are of 60 
and 6 sites in the Ising and the Heisenberg cases respectively. 
In the Ising case the quench is obtained by changing the ex- 
ternal magnetic field from ho to hi = 1. In the Heisenberg 
model the anisotropy parameters is used instead, Ao = 1.5 
while Ai = 0.0,0.2,0.4,0.6,0.8. The time is expressed in 
units of the spin wave velocity. The inset of the left panel 
shows the behaviour of the initial slope as a function of the 
final value of the anisotropy 



from a given point are highly entangled (see Fig 
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FIG. 29 Schematic representation of the dynamics of block 
entropy. Entangled particles are emitted from the region 
A, they will contribute to the block entropy as long as 
one of the two part i cles ends in the region B. [from 
(jCalabrese and Cardvl l2005h ] 

idea is that a point x a in the region A will be entangled 
with a point Xb in the region B if a pair of entangled 
particles emitted at an earlier time arrive simultaneously 
at xa and xb- In this picture the block entropy is sim- 
ply proportional to the length of the interval where this 
can happen. The saturation is reached when the most 
distant quasiparticles (which started in the middle of the 
block) have already got entangled. In the presence of 
disorder the dynamical behaviour is strikingly different. 
This picture applies in the more general context of the 
study of dynamical corr elation functions after a que nch 
as recently discussed in ( Calabrese and Cardvl . 120061 ). 

Possible evidence of localization in the blo ck entropy 
has been discussed in ( De Chiara et all |2005[ ) . The case 
of Ising models in zero magnetic field with long range in- 



teraction the dyna mics was analyze d in the case of factor- 
ized initial states (|Diir et all 120051) . The general picture 
is confirmed also in this case although in the short time 
limit additional oscillations appear probably due to the 
various different phases related to the different couplings. 

We finally mention that the entropy in the case 
of a finite time quenc h has been considered in 
ICherng and Levitovl . 120061 and by ICincio et ail 120071 
Very recently rigorous bounds for the time evolution 
of the block entropy were obta ined by Osborne and 
Eis ter (lEisert and Osborne] . 120061 ) and by Verstraete et 
al (jBravvi et aLl . l2006f ). 



3. Chaos & dynamics of entanglement 

The evolution of entanglement is very sensitive to 
the different properties of the spectrum in the tran- 
sition to chaos as in the case of a quantum com- 
puter with imperfections. A small inaccuracy in the 
coupling constants induces errors in gating or a un- 
wanted time evolution if the Hamiltonian cannot be 
switched exactly to zero. If the imperfection strength 
increases, new phenomena occur and above a certain 
threshold the core of the computer can even "mel t " due 
to the setting in of chaotic behavior ( Benenti et aZ.I . [200Tt 
iGeorgeot and Shepelvanskvl l2000f) . In addition to the 
understanding of the behaviour of the fidelity as an in- 
dicator t o measure the sta bility of the quantum mem- 
ory (see iGorin et all 120061 for a review), a more com- 
plete characterization has recently included the behav- 
ior of entanglement on approaching the transition to 
the ch aos either by considering a dynamics of a disor- 
dered (IMontangero et aE . 20031 : iMontangero and Violal . 
120061: ISen(De) et all 120061) or time-dependent Ising 
model ( Lakshmin aravan and Subrahmavaml 120051 see 



2007 for a recent 



iProsenl 

namical complexity analysed 
model) an d by studying the dynamics of a quan- 
tum map (iBandvopadhvav and Lakshminaravan . 2002 ; 



review on the dy- 
on the kicked Ising 



Bettelli and Shepelvanskvl. 120031: iGhose and Sanders 

2004l)~ 



2004t iMiller and Sarkarl . Il999t iRossini et all I2004D . In 
partic ular the disordered Ising model has been pro- 
posed ( Georgeot and Shepelvansk v. 2000) to describe the 
hardware of a quantum computer, in which system im- 
perfections generate unwanted inter-qubit couplings and 
energy fluctuations. Three different regimes appear de- 
pending on the variance of the fluctuations of the ex- 
change couplings. At a first critical value 6 J c the systems 
crosses from the perturbative to the chaotic regime, while 
at the second point<5 Je the system goes into an ergodic 
state characterized by a Gaussian local density of states. 
All these dynamical regimes h ave been detected in the 
dynamics of entanglement in ( Montangero et all , 120031 : 
IMontangero and Violal 120061 ). 

We finally me n tion the interesting connec- 
tions ( Boness et all [2006) found for the properties 



of the entanglement in a Heisenberg chain with pulsed 



53 



magnetic field with the localization behaviour of the 
quantum kicked rotator. 



B. Generation of entanglement 

The Hamiltonians discussed so far can be also used to 
generate particular entangled states. 

The simplest case is to consider the AX-model in 
the sector in which only one spin is flipped. The fact 
that the Hilbert space is spanned by the basis set \l)) 
hints to the fact that W-states can be generated. For 
short chains, N = 3,4, generalized W-states of the type 
\W) = ^j[|1000) - i|0100) - |0010) - i|0001>] appear 
at discrete times is the initial state of the chain is fully 
polarized ( Wand . [20011 ). This simple scheme cannot be 
extended easily to an arbitrary number of qubits. The 
reason is related to the fact that for generic terms the 
various frequencies appearing in the dynamical evolution 
of the state are not commensurate. 

An interesting example of generation of entangled 
state s is that of cluster states ( Briegel and Raussendorj . 
2001) which has been shown to be essential for the one- 
way quantum computation (iRaussendorf and Briegel 
12001b iRaussendorf et all l2003f ) fsee iHein et all 12009 for 
a review). Remarkably, they can be generated by certain 
spin Hamiltonians. In fact, the dynamical evolution of 
an Ising model in zero magnetic held, H — . Jijafaj, 
is equivalent to a series of conditional phase shifts. In the 
case of a regular lattice this Hamiltonian generates clus- 
ter states which are readily generalized to graph states 
for an Ising model defined on a simple graph. It is easy 
to realize that the evolution operator at time Jt = it/ 4 
can be written asC/(7r/4) = fl i . i(l + ofcr|) If the initial 
state is a product state with all the spin pointing along 
the x-direction, than at time Jt — tt/4 the state is a graph 
state characterized by a maximal bipartite entanglement 
and by maximal entanglement persistence. An example 
of cluster state generate by an Ising chain with nearest- 
neighbour couplings, for N — 4, is given by \*£) c i = 
|0000) + |1100) + |0011) - |1111) (for N = 2,3 cluster 
states coincide with Bell and GHZ states respectively). 
The Ising model is not the only case when graph states 
can b e created. Borhani and Loss ( Borhani and Loss) . 
2005) showed how to gener ate them using He isenberg in- 
teraction while Clark et al ( Cl ark et aD,[2005) considered 
the XX-mode\ Hamiltonian. 

An appropriate tailoring of the initial state or the 
spin graph can lead to the production of prope rly tai- 
lored entangled states ( Koniorczvk et all , 120051 ). In a 
XX-mode\ in a star network it is possible to control 
the concurrence between two spins by varying the ini- 
tial state of the central spin only. Such a spin me- 
diates the int eraction between t he ou ter ones as also 
discussed in ( Hutton and Bosel . 120041 ). The pairwise 
entanglement can be maximized by choosing all the 
outer spin down and the central spin up. These dy- 
namically generated states saturate the CKW inequality 



Eq.® (|Coffman et all |2000D and hence have the maxi- 
mal possible two-site entanglement. Koniorczyk et al an- 
alyzed also the concurrence of assistance, i.e. the max- 
imum amount of entanglement which could be concen- 
trated on two qubits by means of assisted measurements 
on the rest of the system. Depending on the initial sys- 
tem, the behaviour of the concurrence and the concur- 
rence of assistance can be quite different. 



C. Extraction of entanglement 

The entanglement naturally contained in a many-body 
state can in principle be extracted and therefore used for 
quantum information processing. This means that such 
entanglement can be transferred to a pair of particles 
and subsequently used, in principle, for quantum com- 
putation or to test the violation of the Bell's inequali- 
ties. De Chiara et al proposed a scheme of entanglement 
swapping by means of scattering of a pair of particles 
with a spin chain or an optical lattice ( De Chiara et all , 
200G). To this end one sends simultaneously a pair of 



probe particles toward the entangled spin chain in such 
a way that each probe interacts with a different spin. 
The entanglement between the probes has been extracted 
from the spin chain and cannot exist without entangle- 
ment in the chain. This is a genuine non local pro- 
cess between the two probes like in the case of entan- 
glement swapping. In practice the scattering interaction 
between probes and spins in the chain must be capa- 
ble of (at least partially) swap their state. This is the 
case of very common interactions like the Heisenberg or 
the XY ones. The most natural way to extract entan- 
glement from entangled electron spins in solids would 
be to scatter pairs of neutrons off the solid. Another 
possibility would be to realize Hamiltonians of entangled 
spin chains or ladders that can be realized using trapped 



cold atoms dDuan et al. 



120041 ; iGarcia-Ripoll etal 



20031: 



Garcia-Ripoll and Cirac, 



20031) and, as probe particles 
one can use marker qubits ( Calarco et all 120041 ) . 



D. Time evolution of the entanglement in Gaussian states 

The dy namics of Gaussian states has hrst been dis- 
cussed in ( Plenio et all |2004[ ) , where essentially two as- 
pects have been highlighted: the creation of entangle- 
ment from an initially disentangled state and the prop- 
agation of an entangled state on top of a disentangled 
background, both induced by Hamiltonian dynamics. 

The initially entangled two-mode state is characterized 
by the squeezing parameter r, which appears in the co- 
variance matrix as V^ a ^ a = coshr for all phase space 
variables of the zeroth and the first oscillator mode, 
whereas V qo _ qi — V Po _ Pl — sinhr; in absence of further off- 
diagonal elements, all other diagonal elements are equal 
to 1. Two different types of nearest neighbor couplings of 
the oscillator have been considered: ideal springs obey- 
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ing Hooke's law, and its Rotating Wave Approximation 
(RWA). In this approximation both the ki netic and po- 
tenti al energy terms assume the same form (jPlenio et all 
120041) . The RWA conserves the number of bosonic exci- 
tations in the system, whereas the ideal spring coupling 
does not. The initial entanglement is encoded in a zeroth 
oscillator, itself decoupled from the chain of oscillators, 
and one oscillator within the chain as described above. In 
this case, the Hamiltonian for the harmonic lattice (with 
an appropriate choice of the matrix U, see Eq. (|34| ) reads 



H= g (q +Po + 



N \ 

J2(9l+Pt + Kl(q k+ l-qk) 2 )\ (75) 
fe=l / 



After discarding the terms da and a) a\ in RWA, the _ 
can be written (up to a constant) as ( Plenio et all [200 



Hi 



RWA 



fc=l 



(76) 

For both ([75]) and (|76|) the time evolution of the position 
and momentum operator can be evaluated analytically. 

The initially entangled state of two oscillators - a de- 
coupled and one harmonic oscillator within a periodic 
chain - is released at time t = into the background 
of initially disentangled oscillators, all being prepared in 
their ground state. In regard to the pairwise entangle- 
ment we can observe that the n-th and the zeroth os- 
cillators will become entangled after a finite time t c (n), 
which is essentially given by the velocity of sound of the 
underlying model Hamiltonian. After this entanglement 
wave arrives, the n-th oscillators' entanglement exhibits 
damped oscillations with the characteristic time scale of 
the model. For the model ([75)1 , the velocity of sound has 
been determined as v — K\j\J\ + 2K\. Within the RWA 
this velocity appears enhanced: wrwa = K\. The height 
En f of the first local maximum in the logarithmic nega- 
tivity has been used to define the transmission efficiency 
T e ff for entanglement in the chain as T e s = En f / 'En j, 
where Eni is the logarithmic negativity initially prepared 
between the oscillators zero and one. 

In both models En / is observed to saturate when 
cranking up r and hence the initial entanglement. As 
expected, the saturation value is a decreasing function 
of the distance n. The main difference is that for the 
generic coupling the saturation value also depends on the 
coupling strength K\ and is the smaller, the stronger is 
the coupling. Within the RWA it is independent of the 
coupling strength. Given the initial logarithmic nega- 
tivity Eni = — 2r/ln2, the transmission efficiency be- 
haves differently. For the generic model it is zero for zero 
squeezing, has a maximum at a finite squeezing r max 
and anishes for r —> oo. In sharp contrast, in RWA 
T e g is a monotonically decreasing function of r. Inter- 
estingly enough, for the model in RWA, the oscillator 
frequencies and coupling strengths can be optimized as 
to have perfect transmission from one end of an open 



chain of M sites to the other end. This is a chieved by 
virtue of Uj J+ i = Uj+ij = K\^Jn{M — n), Ujj = 1, 
u>„ = 1 — K\ y/nE — K\ y/(n— l)(h + 1) where we defined 
n = M — n. The same turns out to be impossible for 
the generic model and M > 2. The only possibility is to 
choose the couplings and frequencies as in the RWA case 
in the limit of K\ — > (where the RWA is exact). In this 
way, the transmission efficiency can be pushed arbitrar- 
ily close to perfect transmission but with a transmission 
speed tending to zero. The transmission of quantum in- 
formation has been found to be robust against noise in 
Ki for both models (jPlenio et aLl . l2004h . 

Another effect occurring in the generic model is entan- 
glement creation from a disentangled state, which is not 
an eigenstate. This can be realized by an abrupt change 
of the coupling strength. As in the spin case (see Sec- 
tion HX.Ap no entanglement creation can be generated 
within the RWA. In an open chain the oscillators at the 
end points become entangled after half the time a sig- 
nal needs to travel through all the chain. This indicates 
that this initial pairwise entanglement is mediated from 
the center of the chain. Actually this is the first signal 
that possibly can arrive when essentially pairwise entan- 
glement is created or contributes to the eventual pairwise 
entanglement of the boundary oscillators. Raising a fi- 
nite temperature for a thermal state, the arrival time is 
slightly enhanced and the signal is subject to damping. 
However, this effect turned out to be much more sensitive 
to noise in the coupling strength than the propagation of 
quantum information (Ple nio et aZ.I . [200^ ) . 



X. CONCLUSIONS AND OUTLOOK 

During the last years it became evident that quan- 
tum information offers powerful instruments to grasp the 
properties of complex many-body systems. This is the 
reason why this area of research is undergoing an im- 
pressive expansion. In this review we touched only one 
particular aspect of this activity, namely the properties 
of entanglement in many-body systems. 

As mentioned in the introduction, there are several 
important aspects which, for space limitations, were ig- 
nored. In this respect we should certainly mention the 
increasing interest towards the optimization of numeri- 
cal simulations of quantum systems. There were already 
a number of breakthroughs (see also the introduction) 
that, for example, lead to the development of the time- 
dependent DMRG. The design of variational methods to 
study the ground state and finite temperature proper- 
ties of many-body Hamiltonians has been exploited in 
numerous interesting works. We already mentioned in 
the introduction the large body of activity on quantum 
state transfer in spin chains. Here the knowledge of the 
low-lying excitations of the chain (spin- waves) has helped 
in finding new quantum protocols. More is expected to 
come in the next future. 

It is very tempting, although very difficult, to try to 
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envisage the possible evolution of this line of research. 
The study of the topological entanglement entropy is a 
very important aspect that will be investigated in the 
coming future. Adiabatic quantum computation and 
the one-way quantum computation will benefit from the 
study of entanglement in complex systems. The study 
of the topological entanglement may also have remark- 
able spin off to understand many puzzling phenomena in 
condensed matter physics amoung which high T c super- 
condu ctors and heavy fermions are o f paramount impor- 
tance (|Coleman and Schofieldl . 120051) . Furthermore will 
the analysis of new measures for multipartite entangle- 
ment provide with additional insight necessary for under- 
standing the role of entanglement in such complex phe- 
nomena, which might also reveal valuable information, 
e.g. on how to simulate these systems efficiently on a 
computer. 

Many interesting results have been already obtained, 
but the overlap between quantum information and many- 
body physics is not yet fully unveilled. There is a number 
of open questions which provide a fertile ground for a field 
of lively exciting research. 
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